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SUMMARY 


Calculation procedures and computer codes for the thermal design of 
multiple-tube potassium condensers and boilers are presented. 

The condenser code provides for using either sodium, lithium or NaK as 
the shell-side cooling fluid. The main output from the condenser code is 
the condenser tube length required to satisfy specified design-point conditions. 
Additional output includes the local values of vapor quality, potassium tem- 
perature and pressure, and thermal parameters such as the local heat flux and 
heat transfer coefficients along the condenser tubes. 

The boiler code uses either liquid sodium or lithium as the shell-side 
heating fluid. The main output from the boiler code is the boiler tube length 
required to satisfy specified design-point conditions. The boiler code is 
based on using overall average values for the potassium boiling heat transfer 
coefficient from 0% to 100% quality. 

Both the condenser code and the boiler code are programmed in FORTRAN IV, 
and the source languages are compatible with either the IBM-7094 or the 
G.E. - 625/635 machines. 



I INTRODUCTION 


Calculation procedures and digital computer codes, developed for 
preliminary thermal design studies of multiple-tube condensers and once- 
through boilers applicable to Rankine cycle space power systems using 
potassium as the working fluid, are presented in this report. This work 
was done under NASA Contract NAS 3-2528 as a supporting effort to the 
potassium multiple-tube test condenser design study described in Reference 1. 

The condenser computer code calculates the condenser tube length required 
to satisfy specified design-point conditions. In addition, local values of the 
thermal parameters, such as potassium temperature and pressure, local heat 
flux and heat transfer coefficients are computed as functions of distance along 
the condenser tubes. The local potassium condensing heat transfer coefficient 
results reported in Reference 5 are used as the basis for the condenser code. 

The boiler computer code employs average boiling potassium heat transfer 
coefficients, averaged over the entire boiling region from 0% quality to 100% 
quality, to calculate the boiler tube length required to satisfy specified 
design-point conditions. The average potassium boiling heat transfer coefficient 
results reported in Reference 11 are used as the basis for the boiler code. 

The basic thermal and fluid dynamic equations used in the codes, the 
formats of the code inputs and outputs, and illustrative test cases are given in 
following Sections of this report. Listings of the computer codes from which 
duplicate program cards can be made are given in the Appendices. 
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II POTASSIUM CONDENSER THERMAL DESIGN PROGRAM 


The condenser thermal design code is applicable to potassium condensers 
using uniform tubes without inserts in either counterflow or cocurrent flow 
with saturated potassium vapor at the inlet and subcooled liquid potassium at the 
outlet* The code is based on the single-tube condensing potassium heat transfer 
results reported in Reference 5. 

Design Relations 

The system considered is a counterflow or cocurrent flow condenser with 
saturated potassium vapor at the inlet and subcooled liquid potassium at the 
outlet. Distance L is measured in the direction of potassium flow and heat 
transferred from the potassium is taken as positive. Since vapor quality decreases 
in the direction of potassium flow, the differential rate of heat transfer is 
obtained as follows, where changes in vapor kinetic energy and liquid sensible 
heat are ignored. 


Expressing an energy balance between the potassium condensing fluid 
(secondary fluid) and the condenser cooling fluid (primary fluid) for a length- 
increment dL gives after re-arrangement. 



where 


w h^ 

r = 5 fg 

1 tt d. IT U. 
i t i 


(la) 


(lb) 


Distance along the tubes is positive in the direction of potassium flow and 
heat transfer from the potassium is taken as positive. 


Differentiating Equation (la) results in the following equation involving 
the temperature derivatives dT/dx, which will be used in the subsequent analysis, 
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dT 
s 

dx 


dT p 

dx” 


(lc) 


It remains to determine expressions for 
The shell- side coolant temperature 


the -t— terms in Equation (lc). 
dx 

gradient can be written as: 


dT p 

dx" 


+ 


w 

s 



( 2 ) 


where the + sign pertains to counterflow, and the - sign indicates parallel 
flow. 


The secondary temperature gradient may be written as: 


dT dT dP 

s s s dL 

dx ” dP dL dx 
s 


(3) 


where: 

dT 


— - = Slope of saturated vapor temperature - pressure curve 
dP 


dP 


dL 


= Total two-phase pressure gradient 


Ignoring elevation effects, the total two-phase pressure gradient is 
composed of the frictional and momentum contributions and can be expressed as 


dP s _ f G 2 0 _ d_ (1. dx 

dL = ' 2 P f g c D eq ” g c dx A dL 


(4) 


where: 


7T ( “ ) 

dx A 


i_ 

P 


1 + x (K-l) 



r 

P 1 



1 

g 

t 

1 

K 

' p fJ 

T 


a-,) ^ + f 


(K-l) 


(4a) 
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In Equation (4) , the first term on the right hand side is the friction 
pressure gradient and the second term is the momentum or acceleration pressure 
gradient. Note that the friction pressure gradient term is negative whereas 
the momentum pressure gradient is positive (since dx/dL is negative). The net 
effect is that the pressure gradient may be either positive or negative depending 
on the relative magnitude of these terms. 


Substituting Equation (4) into Equation (3), 


dT 

s 

dx 


f G 0 


2 P* g. 


D 

eq 


dT 

s 

dP 

s 


dL 

dx 



dT 

s 

dP 

s 


Then, substituting Equations (2) and (5) into Equation (lc) gives 


(5) 




- f G 


2 P g 
f *c 



dT 

s 

dP 

s 


dL 

dx 


^2 ^ n dT 

5L_ — (i) - + c 

g dx V dP ± 2 

& c A s 


( 6 ) 


where: 


w 

s 



(+) = Parallel flow 
(-) = Counterflow 


Equation (6) is the differential equation for the condensing length as 
a function of gravity. The additional equations to be solved along with 
Equation (6) are Equation (5), and the following differential equation: 


dq 

dx 


w h „ 
s fg 


(7) 


These latter two relations yield the local secondary saturation 
temperature and rate of heat transfer as functions of local quality. 
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Equations (5) to (7) may be considered a set of linear differential 

equations with non-constant coefficients. Due to the complexity of the 

system^ numerical methods are employed for the integrations. In order to 

reduce the system to a set of equations more amenable to computer solution^ 

dL 

it is appropriate to make the variable transformation J = — . The system 
is then reduced to first order in the dependent variables h, J, T g and q, and 
becomes 


dL 

dx 


= J 


(8a) 


= A(x) J + B(x) 
dx 


(8b) 


dT 

c 

dx 


- f G 0 

2 p . g D 
f c eq 


dT _2 . _ dT 

j _ 2_ d_ ,1. * £ 

dP g dx ' dP 

s c As 


(8c) 


dq 

dx 


w h „ 
s fg 


(8d) 


where: 

A(x) = 


(T s - Tp) 


2 A dT 

f G 0 s 

2 p. g D dP 

f & c eq s 


(8e) 


B(x) = 


(T g - V 


2! <L 

g c dx 


, dT 
/In S 

(_) dP~ ± C 2 
s 


$ 


(8f) 


The initial conditions at x = x. are: 

in 


“C- 


J = 


< T S - Tp) 

in 


(9a) 


L = 0 


(9b) 


- 8 - 



s. 

m 


(9c) 


T = T 
s 


q = 0 


(9d) 


Equations 8 are integrated numerically in the condenser code using the 
5-point Lanczos scheme (Reference 2). This technique was chosen for its 
inherent stability and capacity for yielding minimum round-off error. The 
first four steps are taken using the Runge-Kutt a-Gi 11 method (Reference 3) to 
generate 5 points for the more accurate and faster Lanczos equations. Step 
size changes (halving and doubling) are made on the basis of pre-set error 
criteria. 


The assumed shell-side geometry consists of a tube bundle enclosed by 
a cylindrical shell. The tubes are assumed to be arranged in a triangular 
array. The inner periphery of the shell consists of a cylindrical geometry 
with triangular lands positioned in the spaces between the tubes in the outer 
row. With this arrangement y the net flow area of the shell may be divided 
into hexagonal channels surrounding each tube in the array. This geometry 
is illustrated for a seven-tube case by sketch (a). 



Typical Shell-Side Flow Geometry 

Sketch-a -9- 


For a geometry of this type, the following equation yields the shell-side 
(primary) fluid flow area (A ) . 


tt D 


ot 


Sh 


N, 




TT D 


Ot 


( 10 ) 


The equivalent diameter to be utilized in the shell-side Reynolds 
number may be evaluated as 


D 

eq 


2 

Sh 

Wetted Perimeter 


4 A 
WP 


Sh 


(Ha) 


The total wetted perimeter WP includes wetted surfaces due to the tubes, 
the triangular lands, and the shell circumference not occupied by lands, as 
given by Equation (lib), 


WP = TTnd/ 3N p 
t ot os 


(lib) 


The shell-side heat transfer coefficient is calculated by the method of 
Dwyer and Tu (Reference 4). The equation recommended in Reference 4 for 
parallel flow through rod bundles is 


N = a + e (f N ) 
Nu Pe 


Y 


( 12 ) 


where 


P p 2 

a = 0.93 + 10.81 (^ ) - 2.01 (^ ) 

ot ot 


3 = 0.0252 (-^ ) 
ot 


(12a) 


(12b) 


Y = 0.8 


(12c) 


The quantity ? represents the ratio of the eddy diffusivity for heat transfer 
to the eddy diffusivity for momentum transfer. From Reference 13, 


f = 1 - 


1.82 


N (S /V) * 

Pr m max 


(13) 
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where: 


= f <V 

max 


(14) 


The condensing heat transfer relations are based on the forced convection 
potassium condensing coefficient data reported in Reference 5. The local 
condensing heat transfer resistance may be represented as the sum of a vapor 
phase and a liquid film resistance as follows: 



where: 

h = k„/& (15a) 

f f 

6 = liquid film thickness 


The vapor phase heat transfer coefficient is given as: 

dP P_ 

dT 2T 

where T is in °R 


M S c 

h v “ °c h fg < '2 tt RT ) 


(16) 


It is shown in Reference 5 that the condensation coefficient a for potassium 

c 

vapor is approximately 0.2 over the temperature range from 1100°F to 1400°F. 

The local film thickness 6 for use in Equation (15a) is estimated using 
the Dukler model (Reference 6). In the condenser program, the liquid film 
coefficient is implicitly represented in Nusselt's condensing ratio, which is a 
function of the local film Reynolds number and the local shear exerted by the 
vapor on the liquid film interface. The relevant defining equations are: 
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1/3 





h f N* 2 

t *) = -i (: — ) 

V k f g c 


( 17 ) 


where: 


_ 4(l-x)w 
Re T “ tt D VI 

Lt 


The dimensionless vapor shear, T *, is estimated from the local two-phase 
frictional pressure gradient as: 


T* - 
v 


P ,l, ( ! i s 

8 «c 2 9 t e c °eq 


P £ Cf-) 


D. 1/3 

>(“£> <* 5 ) 


(18) 


Overall heat transfer coefficients are utilized in both the condensing 
region and the subcooling region. The overall heat transfer coefficient based 
on the inside heat transfer area is defined as: 


U. = 


,1 1 _i . 

( h + h D h 5 

s w ot P 


where: 


1 _ P it 0 D ot 

h "2k D.. 

w w it 


(19) 


( 20 ) 


The wall thermal conductivity is assumed to be a linear function of temperature, 
for which the mean wall conductivity is given as: 



or 


(a + bT) dT 


k 

w 


a + 


£ (T + 
2 wo 


T .) 
wi 


( 21 ) 
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where: 


T = Outer wall temperature, °F 
wo 


T . = Inner wall temperature, °F 
wi 

The overall heat transfer coefficient is applied on a local basis in the 
condensing region. 


In the subcooled region, the liquid potassium heat transfer coefficient 
is calculated using: 


N = 3.66 + .0055 N_ (22) 

Nu Pe 

Equation (22) is based on liquid heat transfer coefficient data reported in 
Reference 5. The tube length in the subcooled region is calculated from: 


_ sub 

J sub “ tt j) n U AT T 
it t iL Lm sub 


(23) 


where: AT , = Logarithmic mean temperature difference 

Lm sub , . . , , 

in the subcooled region 

The shell-side fluid pressure drop is composed of the friction losses, 
the radial cross-flow losses, and the entrance and exit expans ion- contract ion 
losses. Only the friction and radial pressure drops are accounted for in the 
condenser design program. It is assumed that the expansion-contraction losses will 
be calculated separately after the inlet-exit geometry is established. 

The frictional pressure losses are given as: 


AP 


f L G 
_ tot s 

fr-shell 2 p„ g D 

f c eq 


where 


D G 

f = 0.316/ (— jj — — ) 


0.25 


(24) 


(24a) 
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The analysis of Reference 7 is utilized to compute the radial pressure losses 
in the crossflow regions at the primary fluid inlet and exit. The resulting 
expressions for the radial pressure loss in the shell are: 


AP 


radial 


2 P g 
K f g c 


Oo 25 + 


0. 118 


<f - D 1 * 08 

o 


D 

(— ) 
^ li ; 


- * 16 


n 


y (g .) 

L-> max i 


1.84 


( 25 ) 


i=l 


The maximum mass velocity in any tube row is given as a function of the total 
flow rate and the shell geometry. 

3(i-l) 2 

(3n 2 - 3n + 1) 

max i ” 6(i-l) (p-D )L 

ot c 



where: n = the total number of tube rows 

i = the summation index on the tube rows 

L = the length along the tubes over which crossflow occurs 
c 


The homogeneous model K = 1 is used to calculate the two-phase friction 
pressure drop multiplier. The equation for the local two-phase friction 
multiplier is (Reference 8): 


0 = 


1 + X 


{> 


1 + X 


iW - j 

«w - 1 


(27) 


Thermodynamic properties of the saturated potassium liquid and vapor are 
taken from Reference 9, Transport properties of the potassium liquid and 
vapor are taken from Reference 10, and are compiled directly into the program 
as tabular functions of temperature. 


The shell-side cooling fluid may be NaK, sodium or lithium. The properties 
of these fluids are taken from Reference 10 and are programmed directly into the 
code as tabular functions of temperature. 
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Computer Code 

The condenser design code is programmed in FORTRAN-IV. The source programs 
are compatible with both the IBM-7094 and the GE-625/635 machines. The subsequent 
discussion, however, is limited to programs run on the IBM-7094. 

The code is composed of a main program and a series of shorter subprograms. 
Common data storage is utilized to facilitate data transmission between the main 
program and the subroutines. The majority of the numerical calculations are per- 
formed in the subroutines. The main program is primarily utilized to call the 
calculation routines, read the input, and transmit the output data to the 
peripheral I/O media. The individual programs and their respective primary 
functions are summarized in Table I. A detailed listing of the main code and 
the subroutines is given in Appendix A. 

Data are read into the program by means of the standard "namelist" routine. 
The input consists primarily of the shell and tube geometry, the total power, 
the primary flow rate and inlet temperature, the secondary inlet temperature 
and quality, and the level of subcooling at the condenser exit. Additional 
input includes the starting increment in quality, the flow direction of the 
primary and secondary fluids, and the parameters pertinent to the per step 
error of the integration routine. 

The output produced by the program consists of the total tube length and the 
shell and tube parameters which provide a measure of the design point performance. 
The printout of local conditions along the condenser is optional. If desired, 
the local values of quality, length, static pressure, primary and secondary 
temperatures, heat transfer coefficients, secondary film Reynolds number, heat 
flux and cumulative heat transferred are printed. In this case, the interval 
between local print-points may be specified as 1, 2, 4 or 5. The first page of 
output contains the tube length and the relevant design point performance 
information. The local conditions are printed consecutively on the subsequent 
format pages. 
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In the numerical integration of a system of ordinary differential 
equations by a conventional predictor-corrector technique, the truncation 
error in the i-th dependent variables refers to the per-step difference 
between the predicted and corrected values of the particular variable being 
integrated. For the Lanczos scheme, which is a pseudo predictor-corrector 
method (Reference 2), this type of error is given as: 


E. = 0.11 

i 



(28) 


where the superscripts c and P designate corrected and predicted values, 
respectively, of the y - variable being calculated. These errors may be 
treated independently by requiring the individual errors to be less than a 
set of corresponding maximum tolerable error limits. However, in the condenser 
design program, it has been found to be more convenient to use a weighted error 
for the entire system of differential equations as recommended in Reference 3. 


E 


tot 


0 . 11 



W. 

l 



(29) 


The weighting factors (W.) and an allowable maximum error limit (E ) are 

i max 

read into the program at object time. Within the program, an additional minimum 
error bound is established as: 


E . = E / 40 (30) 

mm max 

For a given incremental step in quality, E^ q ^ is tested against E^^. If 

E, , > E , the step size is halved and the step is repeated. If E, , is still 
tot max 7 tot 

greater than E , the step size is again halved, and the integration is restarted 
max 

at this location. Five such restarts are allowed for each case being calculated. 

If E^ , < E it is then tested against E . . Then, if E, , ^ E . , the integra- 
tot max min tot min' 

tion is continued with the current step size. If, on the other hand, E^^ < E min 7 
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the current step size is doubled, and the next step is taken with Ax^ ew = 2 Ax q ^. 
Doubling is restricted to every fourth step to prevent too rapid a growth of the 
step size in regions of small truncation error. 

In addition to truncation errors, round-off errors can accumulate as the 
solution proceeds along the length of the condenser. Both the Runge-Kutta Gill 
starting routine and the Lanczos technique are formulated to minimize the 
round-off errors. 

When differential equations are approximated by a difference expression, 
as in the present case, extraneous solutions are normally carried along with 
the desired solution. This situation occurs because the approximating difference 
equations are of higher order than the original differential equation. An 
integration procedure may become unstable when errors cause these unwanted 
solutions to grow exponentially. The ultimate result is that the true solution 
is swamped by the extraneous solutions and is completely lost. The stability 
criterion that must be satisfied for the Lanczos method is: 


dy i 

^ ( dT> < 

I . 


0.3, 


i 


1 , 4 


(31) 


In the condenser design program, a stability check is made after each quality 
step. If the conditions of Equation (31) are not satisfied, the step size 
in quality is halved and the step is repeated. If the system is still unstable, 
the step size is again halved, and the integration is restarted using the current 
values of the dependent variables. As in the case of the truncation errors, 
five restarts are allowed for each case. 


Data are put into the program by means of the standard FORTRAN-IV "namelist" 
routine. Input data are loaded into the machine on punched cards. The card 
input format is given below. 
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CARD INPUT FORM 


Column 2 
$ INPUT 


TITLE = 48H 

(48 Hollerith 

Characters) , 


NT = 

y 

ITP = , NCF = 

, NQ = 

MINC = , ITWMAX = 

, JHP 

DT = 

y 

XW = , DCB - 

= y PODT : 

= WP = , GS = 

y 

POWER = 


, TSATIN = 

, TPIN = 

, TSUBO = , HP = 

~y 

HC0N = _ 


* 1 

li 

> 

= , BK = 

, RADPL = , 


SLIPEX = 


, DELTAY = 

PRINT = 

, TOLP = , TOLW = 

y 

ERLIMT = 


, XIN = 

SIGMAC = 

, WGHT = 

$ 


The input variables and relevant units are defined in Table II. Variables 
beginning with the characters I, J, N, M must be entered as integer variables; 
viz ; (NT = 50). The remainder of the input consists of real variables. These 
data may be entered using either an F or an E type format; viz, (POWER = 300), 
or (POWER = 3. E + 2). Change cases may be run by simply entering new values 
for the variables being changed. These new input values are stacked behind 
the initial case input cards. 

Some special instructions are necessary for some of the variables in the 
input list; as given in the following table. Where pertinent, values found to 
be satisfactory for the cases run to date are indicated. 


Symbol 




Ins 

ITP 

Sati 

sf actory 

value r= 

20 

NCF 

1) 

Enter NCF 

- 0 if 

count erf 1. ow 


2) 

Enter NCF 

= 1 if 

parallel flow 

NQ 

Sati 

s factory 

value - 

20 

ITWMAX 

S at i 

s factory 

value - 

30 
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Symbol Inst ruction 

JHP 1) Enter JHP = 0 if primary heat transfer coefficient is to be 

calculated by program 

2) Enter JHP = 1 if primary heat transfer coefficient is input 
DCB Enter DCB = 0 for tube with no centerbody 

GS Enter GS = 0 if number of tubes is specified. This value must be 

entered for every case . 

HP Enter HP = 0 if primary heat transfer coefficient is computed by 

the program 

HCON Enter HCON = 0 if condensing heat transfer coefficient is computed 

by program 

PRINT 1) Enter PRINT = 0 if local condition printout is not desired 
2) Enter PRINT = 1 if local condition printout is desired 
TOLP Satisfactory value = .02 (2 percent accuracy) 

DELTAY Satisfactory value = - .001 

ERLIMT Satisfactory value = .0005 (.05 percent accuracy) for DELTAY = - .001 

WGHT Satisfactory values are WGHT = .2, .5, . 1, .2 for DELTAY = - .001 


The complete input deck for an IBM-7094 calculation consists of the 
requisite IBSYS monitor control cards, the binary decks (main and subprograms), 
and the input data. The following deck set-up pertains specifically to the 
IBM-7094 computer installed at General Electric Company in Evendale. 


Yellow Sequence Card 


DECK SET-UP (IBM- 7094) 


Pink ID Card 


Col. 1 


I STOP 
13L 


$J0B 
$REWIND 
$EXECUTE 
$ IB JOB 

(Binary Programs) 


$DATA 

(Input Data Cards) 

gl* - — (Punch in Col. 1) 

$ IBSYS 

$ENDFILE 

$ENDFILE 

$ENDFILE 

$STOP 


Col. 7 

IBSYS 

IBSYS 


IBJOB 
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Col. 16 


SYSPP1 

GO, MAP, DLOGIC 


SYSPP1 

SYSPP1 

SYSPP1 



Certain types of program errors will cause the execution of a particular 
case to be aborted. If these situations occur, appropriate error comments 
will be printed before the case is terminated. In most instances, these comments 
will be accompanied by the program results up to the point where the error was 
encountered. Error printouts will occur for the following situations: 

(1) If the exit temperature from the condensing region is less than the specified 
subcooled outlet temperature, the outlet temperature is lowered, and the case is 
re-run with the adjusted subcooled outlet temperature. If the same condition 
prevails at the end of this pass, an error comment is printed and the case is 
terminated. 

(2) The primary temperature profile is established by nodal iteration of the 
coolant heat balance relation. If any of the iterations fail to converge, an 
appropriate error comment is printed. The node at which non- convergence occurred 
is included in the printout. 

(3) If the mean wall temperature iteration fails to converge, the case will be 
terminated. The error comment in this situation includes the local quality at 
which the iteration failed to converge. 

(4) Under certain conditions the condensing fluid temperature may become lower 
than the primary coolant temperature. This temperature cross indicates a physically 
impossible situation in the condenser. If a cross occurs, an error comment to 

this effect is printed, along with the computed values of local quality, length 
and temperatures. 

(5) The program will be terminated if the total number of restarts exceeds 5. 

In this instance, the final variable errors and the local values of quality, length, 
and temperatures will be printed before an exit is made from the program. 


- 20 - 



Sample Case 

Results of some sample calculations are given to illustrate the use of 
the program and to show typical output produced by the program. The specified 
design conditions for these sample calculations are: 


Tube material 

Number of tubes 

Tube OD, inches 

Wall thickness, inches 

Tube, p/D , 

7 ot 

Power, KW 
Secondary Fluid 

Potassium inlet temperature, °F 

Potassium inlet quality 

Potassium exit temperature, °F 

Primary coolant 

NaK inlet temperature, °F 

NaK flow rate, lb /sec 
m 

Initial quality increment Ax 
Maximum allowable error 
Average wall temperature tolerance 
Length of shell-side radial pressu 
Condensation coefficient 


Stainless steel-316 
19 

0.6875 

0.035 

1.3 

225 

Potassium 

1300 

1.0 

1191 

NaK (counterflow) 
1190 
10.2 
-.001 

.0005 (or .05%) 
.0005 

drop region 0 

0.2 


The printed output for this sample calculation is shown in Figures 1 and 

2 . 


The programs used in the condenser code are listed in Table I, the nomenclature 
for the condenser code input is given in Table II and the nomenclature for the 
printed output is given in Table III (Appendix A). 
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CONDENSER OESIGN (GENERAL PARAMETERS) 
TITLE* POWER-225KW, TSAT-1300, 10 APPROACH 


GEOMETRY 


NT 


X 

19 

COND LENGTH 

= 2. 

6834 FEET 

TUBE ID 


* 0 • 

6175 INCHES 

SUBC LENGTH 

* 0. 

0079 FEET 

TUBE OD 


* 0 • 

6875 INCHES 

TOTL LENGTH 

= 2. 

6914 FEET 

MALL XW 


* 0 . 

0350 INCHES 

TUBE P/D 


* 

1.30 

TUBE FLOW 

AREA* 0.03951 SQ FT 

SHELL ID 


* 4 

.091 INCHES 

CENTER BODY 

o 

o 

N 

O 

• 

INCHES 

SHELL FLOW AREA* 0. 

0423 SQ FT 




THERMAL PARAMETERS 




POWER KW 
225.00 


QTOT f BTU/SEC 
213.25 

QCOND»BTU/SEC 

207.93 

QSUB 

, BTU/SEC 
5.32 


PRIMARY 



SECONDARY 


TPIN 

X 

1190.00 

DEG F 

TSIN 

3 

1300.00 

DEG F 

TPOUT 

X 

1290. OC 

DEG F 

TSATCE 

3 

1304.33 

DEG F 

WP 

3 

10.1360 

LBM/SEC 

TSOUT 

3 

1191.00- 

DEG F 

GP 

3 

239.653 

L8M/FT2SEC 

WS 

3 

0.24386 

LBM/SEC 

VP 

3 

5.403 

FT/SEC 

GS 

3 

6.172 

LBM/FT2SEC 

REP 

3 

83635.6 


VMAX 

- 

319.283 

FT/SEC 





HSL 


1332.85 

B/FT2HROEG 





REL 

= 

3299.8 


DPFRIC 

= 

0.191 

PSIA 

PEL 

- 

11.42 


OPRAD 

3 

0. 

PSIA 

UIL 

= 

842.69 

B/F T2HROEG 





DELTAP 

3 

-0.221 

PSIA 

DPTOT 

3 

0.19 1 

PSIA 

DPFRIC 

= 

0.203 

PSIA 





DPMOM 

3 

-0.425 

PSIA 
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CQNDENS C-R DESIGN (LOCAL CONDITIONS) l 

. 7ITLE= POWER-2 25 KW » TSAT-1300, 10 APPROACH 


X 

z 

TSC 

PLQCAL 

TPC 

HSC 

HP 

UI 

0 

00 A 

REZ 


FT. 

DEG F 

PSI A 

OEG F 

HR UNITS 

HR UNITS 

HR UNITS 

ft/SEC 

8/ F T 2 HR 


1.00000 

0. 

1300 .00 

8.822 

1290.00 

28903.3 

5476.4 

2331.8 

0. 

23326.8 

0. 

0.999C0 

0.0X04 

1299 .99 

8.821 

1289.90 

32173.3 

5476.4 

2351.1 

0.21 

23710.0 

3.4 

0.99800 

0.0207 

1299 .97 

8.821 

1289.80 

32170.9 

5476.5 

2351.1 

0.42 

23916.7 

6.8 

0.99700 

0.0310 

1299.96 

8.820 

1289.70 

32168.6 

5476.5 

2351.1 

0.62 

24115.8 

10.2 

0.99600 

0.0411 

1299 .95 

8.819 

1289.61 

32166.3 

5476.5 

2351.1 

0.83 

243L5.3 

13.6 

0.99500 

0.0512 

1299.94 

8.819 

1289.51 

32164.0 

5476.5 

2351.1 

1.04 

24515.3 

17.0 

0.99400 

0.0612 

1299 .93 

8.818 

1289.41 

32161.7 

5476.5 

2351.1 

1.25 

24715.8 

20.4 

0.99300 

0.0711 

1299 .91 

8.817 

1289.32 

32159.4 

5476.6 

2351.1 

1.46 

24916.6 

23.8 

0.99200 

0.0909 

1299 .90 

8.817 

1289.22 

32157.2 

5476.6 

2351.1 

1.66 

25117.8 

27.2 

0.99100 

0.0907 

1299 .89 

8.816 

1289.12 

32154. 9 

5476.6 

2351.1 

1.87 

25319.4 

30.6 

0.99000 

0. 1004 

1299 .88 

8.816 

1289.02 

32152. 7 

5476.6 

2351.0 

2.08 

25521.5 

34.0 

0.98800 

0.1195 

1299 .86 

8.815 

1208.83 

32148.4 

5476.7 

2351.0 

2.50 

25926.8 

40.8 

0.98600 

0. 1384 

1299.84 

8.813 

1288.63 

32144. 1 

5476. 7 

2351.0 

2.91 

26333.6 

47.6 

0.98400 

0.1570 

1299 .81 

8.012 

1288.44 

32139.8 

5476.0 

2351.0 

3.3 3 

26741 .8 

54.4 

0.98200 

0. 1753 

1299 .79 

8.011 

1208.24 

32135.7 

5476.3 

2351. C 

3.74 

27151.5 

61.2 

0.98000 

0. 1933 

1299.77 

8.810 

1288.05 

32131.6 

5476.0 

2351.0 

4. 16 

2 7562.6 

68.0 

0.97600 

0.2286 

1299 .74 

8 .809 

1287.66 

32123.6 

5476.9 

2 150.9 

4 . 99 

20389.0 

81.6 

0.97200 

0.2629 

1299.70 

8.807 

1287.27 

32115.8 

5477.0 

2350.9 

5.82 

29220.4 

95.3 

0.96800 

0 .2963 

1299 .67 

8.005 

1286.88 

32004. 7 

5477.1 

2350.3 

6.66 

30049.6 

108.9 

0.96400 

0.3287 

1299 .64 

8.803 

1286.49 

31847.0 

5477.2 

2349.5 

7.49 

30879.8 

122.5 

0.96000 

0.3604 

1299 .61 

8 . 802 

1286.10 

31698.8 

5477 . 3 

2 348 . 7 

8.32 

31714,5 

136.1 

0.952C0 

0.4212 

1299 .55 

8.799 

1235.33 

31425.0 

5477.5 

2347.2 

9,98 

3 3396.2 

163.3 

0.94400 

0.4791 

1299 .51 

8. 797 

1284.55 

31175.4 

5477.6 

2345.8 

11.65 

35093.0 

190.5 

0.93600 

0.5343 

1299 .47 

8.795 

1283.77 

30930.3 

5477.8 

2344.5 

13.31 

36802.2 

217.7 

0.92800 

0 .5871 

1299.43 

8.793 

1282.99 

30693. 1 

5478.0 

2343.1 

14.97 

3 8 5 2 2 . 8 

244.9 

0.92000 

0.6376 

1299 .40 

8. 792 

1282.21 

30469. 2 

5478.2 

2341.8 

16.64 

40254.2 

272. 1 

0.90400 

0.7324 

1299 .35 

8.789 

1280.65 

30038.3 

5478.5 

2339.3 

19.97 

4374^,4 

326.5 

0.88800 

0.8201 

1299 .32 

8.788 

1279. 10 

29626.2 

5478.9 

2336.8 

23.29 

47266 . 1 

380.9 

0.87200 

0.9016 

1299 .31 

8.787 

1277.54 

29241.9 

5479.2 

2334.5 

26.62 

50814.9 

435.4 

0.85600 

0.9776 

1299 .30 

8.787 

1275.98 

28880. 1 

5479.5 

2332.2 

29.95 

5 4 3 d 7 . 4 

489.8 

0.84000 

1 . 0490 

1299 .31 

8.787 

1274.42 

2 8537.1 

5479.9 

2330.C 

3 3.28 

5 7979. 1 

544.2 

0.92400 

1.1162 

1299.32 

8.708 

1272.86 

2 8209.9 

5480.2 

2327.8 

3 6 . 6 i 

6 1 5 8 H . I 

598.6 

0. 90800 

1.1796 

1299 .35 

8.789 

1271.31 

27877.5 

5480.6 

2 32 5.6 

3 9.9 3 

65207.4 

6 5 3,0 

0.79200 

1.2397 

1299 .38 

8.790 

1269.75 

27557.2 

5480 . 9 

2323.4 

4 3.26 

68038.3 

7 C 7 * 5 

0. 776C0 

1.2960 

1299 .41 

8.792 

1268.19 

2 7248.0 

5481.3 

2321,3 

4 6.5 9 

7 2 4 7 9.6 

76 1 . 9 

0. 74400 

1.4031 

1299.51 

8.797 

1265.07 

26645.7 

5482.0 

2316.9 

5 3.25 

79 781 .4 

8 70. B 

0.71200 

1.5001 

1299.62 

8.803 

1261.96 

26067.4 

5482.7 

2312.6 

5 9.90 

07106.1 

9 79. 7 

0.680 00 

1 .5398 

1299.75 

8.809 

1258.84 

25491 . 7 

5403.4 

23C8.0 

t> 9 . 6 

° 4432.6 

1 0 8 8 . / 

0.648)0 

1.6725 

1299 .90 

8.817 

1255.72 

24931 . 8 

54B4.1 

2 303.5 

7 3.21 

101764.0 

119 7. o 

0.61600 

1 . 7500 

1300 .06 

8.825 

1252.60 

24388.4 

5484.8 

2298.0 

7 9.8/ 

10 9 09^ . 1 

1 30 >.6 

0.58400 

1 .8225 

1300 .23 

H. 833 

12 4 9.48 

2 3858.3 

5485.5 

2294.2 

86.5 2 

1 16412. 1 

1413.7 

0.55200 

1 .8906 

1 300 .40 

8.842 

1246. 36 

2 3 338.5 

5486.2 

2289.4 

93. Id 

123719.1 

1524. 7 

0. 52000 

1 .95 50 

1 300 .59 

8.052 

124 3.24 

22826.0 

5486.9 

2284.4 

99.0 3 

131007.0 

16 3 3.8 

0.48800 

2.0159 

1 300 .79 

8.862 

1240.12 

22320. 3 

5487.6 

2279.4 

106.49 

130271.8 

1742.9 

0.45600 

2.0739 

1300 .99 

8.872 

1237.00 

21818. 3 

5488.3 

22 74. 1 

113.14 

146508.4 

1852.0 

0.42400 

2 . 1290 

1 301 .19 

8.883 

12 33.88 

21316.7 

5489.0 

2268.7 

119.79 

152 7 09, 7 

1961.2 

0.39200 

2.1817 

1301.41 

8 . 893 

12 30.76 

2G812.4 

5489. 7 

2263.0 

126.46 

159868, 7 

2070. 3 

0.36000 

2.2321 

1301 .63 

8.904 

122 7.64 

20299.6 

5490.4 

2256.9 

133.10 

166974.9 

217 9,6 

0.29600 

2 . 3268 

1302 .07 

8.92 7 

1221.40 

19228.9 

5491 . 8 

2243.2 

146.41 

180976.8 

2398. 1 

0 . 2 32 GO 

2.4143 

1302 .54 

8.951 

1215.15 

18053.2 

5493.2 

2226.5 

159. 71 

194669.0 

2616.8 

0.16800 

2 .4958 

1303 .02 

8.976 

12C8.90 

16751.2 

5494.6 

2205.5 

173.02 

207571.8 

2 8 3 5.5 

0 . 10400 

2.5718 

1 303 .51 

9.001 

1202.66 

16277.4 

5496.0 

2197 . 3 

106.32 

221605.5 

3064 , 5 

0.04000 

2.6433 

1304 .01 

9.027 

1196.41 

15767.6 

5497.4 

2188.0 

199.62 

2 35433 , 9 

3? 7j, 5 

0. 

2.6834 

1 304.33 

9.043 

1192.50 

15113.9 

5498 . 3 

2175,1 

207.93 

24323 J. 8 

3 410, 5 

0. 

C . 

0 . 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 , 


" 

-I'' 

j ) ;• 1 1 1 

>• > r ! v • f 
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Ill POTASSIUM BOILER THERMAL DESIGN PROGRAM 


The boiler thermal design code is applicable to once-through potassium 
boilers using uniform tubes with or without helical inserts in counterflow 
with subcooled liquid potassium at the inlet and saturated or superheated 
potassium vapor at the outlet. The code is based on the single-tube boiling 
potassium heat transfer and pressure drop results reported in Reference 11. 

Design Relations 

Determination of the required boiler tube length for a fixed tube geometry, 
power level, and fluid conditions involves a simultaneous solution of the 
applicable heat transfer and pressure drop relationships. 

The total boiler length may be divided into a liquid region, a boiling region, 
and a superheat region. Distance L along the tubes is positive in the direction 
of potassium flow, and heat transfer to the potassium is taken as positive. The 
respective channel lengths may be related to the total heat transfer rate, the 
regional overall heat transfer coefficients, and the logarithmic temperature 
difference in each region as follows: 


sub 


sub 


U K N ^D. 
sub t 1 


(T - T ) 

V PI sr 


n 


(T - T ) 

v Pout Sin 


(T 


PI 


T ) 

sr 


Pout 


Sin 


(32) 




U N tt D. 
B t i 

< T P2 - 

T S2> - 

(T - T ) 

v pi sr 



in 

<T P2 

- T“ ) 
S2' 




(T 

- T ) 

- 



( 33 ) 
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(34) 


sup 


U N n D. 
sup t 1 


sup 

(T - T ) 
Pin Soul: 


(T 


P2 


T ) 
S2' 


JL 


(T - T ) 
v Pin Sout' 


(T 


P2 


T ) 
S2' 


where the subscripts on the temperatures are keyed to the following sketch. 



T 


Pin 


T 

Sout 


Fluid Temperature Distributions 
Sketch-b 


The remaining heat transfer relations are supplied by the following heat 
balance relations: 


q sub W s C Ps (T sl T Sin^ Cp P (T P1 ^out^ 


(35) 


q B = W s h fg = W P CP P (T P2 " T P1> 


(36) 
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“.up - w s £p , (T Sout ' T s2> - w p 0p P (T Pin • W 


(37) 


^total ^sub + + ^ 


sup 


(38) 


The various heat capacities in the above equations are evaluated at the average 
fluid temperatures in the respective regions. 


There are seven equations and nine unknowns; viz: L , L_> L , q . q„ » 

SUD D Sup 15 SUp 

T } T y T y and T . The remaining two equations are the pressure drop 
P2 Pi S2 Si 

relations for the boiling and the superheat regions. 


The pressure drop in both the boiling and the superheat regions are composed 
of frictional and momentum contributions. 


AP_ = AP„ + AP 


B 


fr 


B 


“b 


AP = AP,* + AP 

sup f r m 

sup sup 


or 


^B = P sl " P s2 


= * H Sh V 0 + V ( I _ i_ ) 

2P g D ° 1 g c P g2 P fl 
f c eq 


(39) 


2 2 

f~L ■ G G " . 

. _ _ _ H supH H H ,1 1 

sup s2 sout 2 p g D g c P gout P g2 

g C eq 


) 


(40) 


The pressures P , and P „ are related to the corresponding temperatures T , and 
sl s2 Si 


sl s2 

T 0 by the saturation properties of the boiling fluid. 
s2 


“ 27 - 



“sup * W S Cp s <T Sout - T s2> - W P Cp P (T Pln ' T P2 ) 


(37) 


q total q sub + ^B + q sup 


(38) 


The various heat capacities in the above equations are evaluated at the average 
fluid temperatures in the respective regions. 


There are seven equations and nine unknowns; viz: L 9 L. } L , q , q„, j 

sub B sup B sup 

T , T , T and T . The remaining two equations are the pressure drop 
P2' PI s2 si 

relations for the boiling and the superheat regions. 


The pressure drop in both the boiling and the superheat regions are composed 
of frictional and momentum contributions. 


AP = AP 

B 


fr. 


+ AP 


B 


"b 


AP = AP. + AP 
sup fr m 

sup sup 


or 


= P sl - P s2 


f H Sh °H 
2P_ f g c D eq 


0 + — S- (i - — ) 

o-i g c ( P g2 P. 


fl 


(39) 


f H L supH °H °H -1 1 

AP — P - P - I— 4 . — 1— ( - — 

sup - s2 sout - 2 - g c p t p 2 

g C eq 


) 


(40) 


The pressures P , and P - are related to the corresponding temperatures T , and 
sl s2 si 

T by the saturation properties of the boiling fluid, 
sz 
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may be either lithium or sodium. The secondary boiling fluid is potassium. 


The thermodynamic and transport properties of lithium and sodium are taken 
from Reference 10 and are programmed directly into the code as tabular functions 
of temperature. The thermodynamic properties of potassium are taken from Reference 9 
and are included as a subroutine in the code. 

The average potassium heat transfer coefficients used in the code for the 
boiling region are supplied as an input. Experimental boiling potassium heat 
transfer coefficients obtained from Reference 11 are used. 

For tubes without inserts, the superheated vapor heat transfer coefficient is 
calculated from the conventional Colburn relation (Reference 12): 

N = 0.023 (N )°' 8 (N_ ) 1/3 (43) 

Nu Re Pr 

In the case of tubes with helical inserts, the coefficient is calculated from 
the empirical relation developed in Reference 11: 

N Nu ' °' 39 <V °' 56 <H Pr >1/3 <44) 

The length mass velocity and equivalent diameter employed in Equation (44) 
are the helical values as given by Equations (41) and (42) . 


The local two-phase pressure drop multiplier is based on the homogeneous 
model (K = 1) as follows (Reference 8): 


0 = 


[<p f /p g ) - 1] 

( i * * [<w - *]}” 


( 45 ) 


Where n = 0.25 (Reference 8). 
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For a uniform axial heat flux distribution in the boiler, quality is linear 
with length, which approximates the actual distribution. Equation (45) can 
be integrated from x = 0 to x = 1, assuming linear variation in quality, to 
yield an average integrated two-phase multiplier for the boiling region: 


d + b) 1_n i + § 


0—1 T 1 (l-n) T) ( 1-n) T 2 


1 1 11 
(n-l)(l+ri) n_1 " (n-2) ( 1-ri) n-2 + (n-2) n-1 


where 


P f 

? = (p 1) and T| = (j-j 1) 

g g 


Computer Code 

The boiler design code is programmed in FORTRAN-IV. The source programs 
are compatible with both the IBM-7094 and GE-625/635 machines. The following 
discussion, however, is limited to programs run on the IBM-7094. 


The code consists of a main program and several smaller subprograms. The 
majority of the numerical calculations are performed in the main program, along 
with the relevant input and output operations. The subroutines are utilized 
primarily to perform repetitive calculations. The main program and the subroutines 
utilized are listed in the following table: 


Program Primary Function 

PBDES (Main Program) 

RDP Computes radial pressure drop 

FRIC Computes helical flow friction factors 


DICHOA 


Iteration routine 

LININT 


Linear interpolation 

THERMA 


NRL Thermodynamic Properties 

(1) 

DUBIN 


(2) 

CURV 



( 46 ) 
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A complete listing of these programs is given in Appendix B. 


The requisite input to the program consists of the shell and tube geometry, 
the total power level, the primary flow rate and inlet temperature, and the 
secondary inlet and exit temperatures. The secondary exit temperature is fixed 
by specifying the exit pressure and the degree of superheat above the saturation 
temperature corresponding to this pressure. 

The output produced by the program consists primarily of the required tube 
length to transfer the specified power at the given temperature levels of the 
primary and secondary fluids. Additional output includes the design point 
performance as indicated by the tube and shell-side pressure drops, the potassium 
temperature at boiling initiation and the radial acceleration produced by the 
helical insert . 

Data are put into the program by means of the standard FORTRAN-IV "namelist” 
routine. Input data are loaded into the machine on punched cards. The card 
input format is given below. 


CARD INPUT FORM 


Column 2 
$ INPUT 

TITLE = 48H 


(48 Hollerith Characters) 

NT = 

, ITP = 

ITS = , ITSUP = 

, JHP = 

POWER = 

, TRINS = 


, TSIN = , DTSUP = , 

PSOUT = 

, DT = 

y 

XW = , PODT = 

, PODI = , DCB 

HB = 

, HSUP = , 

HP 

= , HPMARG = 

, GS = 

WP = 

, XINST = 

y 

RADPL = , TOLP = 

, TOLS = , 

TOLSUP = 

li 

Oh 


CSLOPE = , CS = 

, XN= , 

TW = 

> > 

— 

— up to 20 values in 

Table , 

XKW = 

y y ”” 

— 

up to 20 values in 

Table .$ 
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The input variables and relevant units are identified in Table IV. 
Variables beginning with the characters I, J, N must be entered as integer 
variables; viz, (NT = 50). The remainder of the input consists of real 
variables. These data may be entered using either an F or an E type format. 
Change cases may be run by entering new values for the variables being changed. 

Some special instructions are necessary for some of the variables in the 
input list. These instructions are given in the following table. Where 
pertinent, values found to be satisfactory for the cases run to date are 
indicated. 


Symbol 


Inst ruction 

I TP 

Satisfactory 

value = 20 

ITS 

Satisfactory 

value = 50 

ITSUP 

Satisfactory 

value = 50 


JHP 

1) Enter JHP = 0 if primary heat transfer coefficient 
is to be calculated by program 


2) Enter JHP = 1 if primary heat transfer coefficient 
is input 

DCB 

Enter DCB = 0 for tubes with no insert 

HSUP 

Enter HSUP = 0 if superheated vapor heat transfer 
coefficient is calculated by program 

HP 

Enter HP = 0 if primary heat transfer coefficient 
is calculated by program 

HPMARG 

Set HPMARG = 1 for most cases 

GS 

Enter GS = 0 if number of tubes is specified. 
This value must be entered for every case. 

TOLP 

Satisfactory value = .25°F 

TOLS 

Satisfactory value = .005 

TOLSUP 

Satisfactory value = .003 

C SLOPE 

Set CSLOPE = 1 

CS 

Set CS = 1 
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The deck set-up is the same as for the condenser thermal design code 
(Section II) . 

Execution of the program will be terminated if any of the following types 
of errors occur: 

(1) The primary exit temperature is established by iteration of the primary 
fluid heat balance. If this iteration fails to converge, an error comment is 
printed and the program calls EXIT. 

(2) If the number of iterations for either the boiling length or the superheat 
length are exceeded, subroutine DICHOA will call DUMP. The difference between 
the iterated variables will be printed. Execution of following cases will be 
suppressed. 


Sample Case 

Results of some sample calculations are given to illustrate the use of the 
program and to show typical output produced by the program. The specified design 
conditions for these sample calculations are: 


Tube material 

Number of tubes 

Tube OD, inches 

Tube p /D 

s ot 

Helical insert p/D 
Power, KW 

Primary heating fluid 
Secondary fluid 

Potassium inlet temperature, °F 

Net superheat of potassium outlet, °F 

Potassium exit pressure, psia 

Sodium inlet temperature, °F 

Sodium flow rate, lb /sec 
7 m 

Potassium boiling heat transfer coefficient, 
Btu/ft-hr-°F 

Helical insert vane thickness, inches 

Length of shell-side radial pressure drop 
region 


Haynes-25 alloy 
7 

0.75 

1.3 

2 , 4 and 6 
300 

Sodium 

Potassium 

1200 

50 

56.7 ( 1700° F sat.) 
1850 

9.03 
4000 

0.0625 

0 
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2 

The potassium boiling heat transfer coefficient of 4,000 Btu/hr-ft -°F, 
average for the entire boiling region from 0% to 100% quality, is typical 
of the average boiling potassium heat transfer coefficient data reported in 
Reference 11 for tubes with helical inserts. Further discussion of the use 
of the Reference 11 data for boiler thermal design is given in Reference 1. 

The printed outputs for these sample calculations are shown in Figures 3a, 
3b and 3c. 

The nomenclature for the boiler code input is given in Table IV and the 
nomenclature for the printed output is given in Table V (Appendix B) . 
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POWER BOILER THERMAL/HYDRAULIC DESIGN 
TITLE* 7-TUBE BOILER. CASE 1-B 


GEOMETRY 


NT 

= 

7 


TUBE P/D 

1.30000 

TUBE 

ID 

= 

0.6900 

INCHES 

TU8E FLOW AREA = 

0.018177 FT 2 

TUBE 

OD 

~ 

0.7500 

INCHES 

INSERT P/D 

4.00 

WALL 

XW 

= 

0.0300 

INCHES 

SHELL ID = 

2.709 INCHES 

TUBE 

LTOT 

= 

7. 0434 

FEET 

SHELL FLOW AREA* 

0.0185 FT 2 

TUBE 

LB 

= 

4.6559 

FEET 



TUBE 

LSC 

S 

0.6078 

FEET 




THERMAL PARAMETERS 


POWER 


QTOT 

QB 



QSC 

Q/A BOILING 

KW 


BTU/SEC 

BTU/SEC 


BTU/SEC 

BTU/FT2 HR 

300. 


284. 34 

250.10 



30.56 

152933.29 


PRIMARY 



SECONDARY 


TP IN 

- 

1050.00 

DEG F 

TSIN 

— 

1200.00 

DEG F 

TPOUT 

= 

1750.18 

DEG F 

TSAT 

= 

1717.10 

DEG F 

WP 

= 

9.080 

LBM/SEC 

TSDUT 

= 

1705. 10 

DEG F 

GP 


489.645 

LBM/FT 2SEC 

WS 


0.315 

XffRTSTC 

VP 


10.845 

FT/SEC 

GMAX 

= 

25.500 

LBM/FT2 SEC 

REP 


180115.0 


vMax 

- 

143.26 

FT/SEC 

PEP 

iift‘ 


614.10 


ACCL 

= 

713808.62 

FT/SEC2 

HP 


13483.46 

B/M 2HRDEG 

“ACCL/G 

as 

22 lb7 . 9T 

"GELT - 

HWALL 

= 

6819.52 

B/FT2HR0EG 

HSL 

= 

849. 10 

B/FT2HR DEG 

DELTAP 

= 

1.791 

PSI A 

HB 

- 


BYFT2HR DEG " 

DPR I N 

= 

0 . 

PSI A 

REL 

= 

12797.3 


DPROUT 


0 . 

PSI A 

PEL 

= 

51.93 


... ... . 




UIL 

- 

715.25 

B/FT2HR DEG 





U I B 

= 

2 1 5 1 . r4 

B/FT2HR DEG 

. _ 




DEL TAP 

= 

2.643 

PSI A 





DPFRIC 

~=~ 

1.3 71 

PSI A - - 





DPMOM 

= 

1.2 72 

PSI A 


SUPERHEAT REGION 


Q SUP 

~ ft ti I 7r r r — 

HSUPS 

TSUPS 

TUBE L SUP 

DPSUP 

DPSFRIC 

DPSMOM 

dTU/S EC 

HR UN I I S 

UE(TT~ 

FT 

ptt ” 

~~P ST 

P57 — 

3.68 

57.37 

1750.00 

1.7798 

1.098 

1.027 

0.0709 


Pig upo ya . 


T/pical Computer Printout of Overall Data From The 
Boiler Thermal Design Code (Case 1-B, Insert p/D - i| ) 
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* 

* 


POWER i 

BOILER 

THERMAL/HYDRAULIC DESIGN 




T I TLE= 

• 

7-TUBE BOILER. 

CASE 1-C 





GEUMETRY 


NT = 

7 


TUBE 

P/D 

1.30000 



TUBE 

ID = 

0.6900 

INCHES 

TUBE 

FLOW AREA = 

0.018177 

FT 

2 

TUBE 

OD = 

0.7500 

INCHES 

INSERT P/D 

6. 00 



WALL 

XW = 

0.0300 

INCHES 

SHELL 

ID 

2.709 

INCHES 

TUBE 

LT0T = 

7.0105 

FEET 

SHELL 

FLOW AREA = 

0.0185 

FT 

2 

TUBE 

LB = 

4.4770 

FEET 






TUBE 

LSC = 

0.5858 

FEET 







THERMAL PARAMETERS 


POWER 

0 TO T 

QB 

osc 

Q/A BOILING 

KW 

BTU/SEC 

BTU/SEC 

BTU/SEC 

BTU/FT2 HR 

300. 

284. 34 

250.26 

30.30 

159146.26 



PRIMARY 



SECONDARY 



TP IN 

— 

1850.00 

DEG F 

TSIN 

— 

1200.00 

DEG F 


T POUT 

= 

1750.18 

DEG F 

TSAT 

- 

1712.86 

DEG F 


WP 


9.030 

LBM/SEC 

TSOUT 

= 

1703.95 

DEG F 


GP 

= 

489.645 

LBM/FT2SEC 

WS 

= 

0.315 

LBM/SEC 


VP 

* 

10. 8<t5 

FT/SEC 

GMAX 

= 

22.637 

LBM/FT2 

SEC 

REP 

= 

180115.0 


VM AX 


95.88 

FT/SEC 


PEP 

= 

614.10 


ACCL 

= 

319781.24 

FT/SEC2 


HP 

= 

13483.46 

B/FT2HRDEG 

ACCL/G 

= 

9931.09 

“GEFS 


HWALL 

= 

6816.14 

B/FT2HRDEG 

HSL 

= 

849. 40 

B/FT2HR 

DEG 

DELTAP 

= 

1.783 

PS I A 

Hb 

= 

4000700 

B/FT2HR 

OK 

DPR IN 

= 

0. 

PSI A 

REL 

= 

12785.0 



QPRQUT 

= 

0. 

PSI A 

PEL 

= 

5 1.80 







UIL 

= 

715. 50 

B/FT2HR 

DEG 





UfR 

' = 

2130.80 

B/TT2HR 

DEG 





DELTAP 

- 

1.947 

PSI A 






DPFRIC 

= 

0.940 

PSTA 






dpmom 

= 

1.006 

PSIA 





SUPERHEAT REG 

ION 





Q SUP HSUPS TSUPS TUBE LSUP 

dpsup 


DPSFRIC 

DPSMOM 


BTU/bhC HR 

UNITS DEG 

F FT 

FS1 


PSI 

PSI 


3.77 53 

.66 1 75 C . 

00 1.9478 

0.850 


0. 798 

0.0519 


Figu re 


Typical Computer Printout of Overall Data From The 
Boiler Thermal Design Code (Case 1-0, Insert p/D = e) 
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POWER 9GILFR THERMAL / HYDRAULIC DESIGN 
TITLE* 7-TUBE BOILER. CASE 1-F 


GEOMETRY 

MT * 7 TUBE P/D = 1.30000 

TUBE IC * 0.6900 INCHES TUBE FLOW AREA = 0.C18177 FT 2 

Tuee oc = 0.7500 inches insert p/d * 2.00 

WALL XW * 0.03^0 INCHES SHELL ID = 2.709 INCHES 

TUBE LTCT* 9.6543 FEET SHELL FLOW ARE A = 0.0185 FT 2 

TUBE LB * 7 .4808 FEET 

TUBE L SC * 0.9115 FEET 


THERMAL PARAMETERS 

PCWER QTOT 30 OSC C/A BCILING 

KW BTU/SEC BTU/SEC BTU/SFC BTU/FT2 HR 

300. 2B4.34 248.67 32.46 94638.28 


PR IMARY 


T P IN 

= 

1850.00 

DEG F 

TPCUT 


1750.18 

DEG F 

WP 

= 

9.0R0 

L8M/SEC 

GP 

= 

439.645 

LBM/FT2SEC 

VP 

= 

10.845 

FT/SEC 

REP 


130115.0 


PEP 

= 

614.10 


HP 

- 

13483.46 

13/ FT 2HRDEC 

HWALL 

= 

6 q 3 7 . 90 

B/FT2HRDEG 

CELT AP 

= 

2.455 

PS I A 

DPRIN 

= 

0 . 

PSIA 

DPROUT 

= 

0 . 

PSIA 


SECONDARY 


TSIN 

= 

120C. 00 

DEG F 


TSAT 


1 749.03 

DEG F 


TSOUT 


1 71C.89 

OEG F 


WS 

— 

0.315 

LBM/SEC 


GMAX 

= 

37.343 

LBM/FT2 

SFC 

VMA X 

= 

28C.84 

FT/SFC 


ACC L 


743367.4 7 

FT/SEC2 


ACCL/G 

= 

85197.75 

GEES 


HSL 

s 

847.59 

B/FT2HR 

DEG 

HB 

= 

40CC.00 

R/F T2HR 

OEG 

REL 


12884.2 



PEL 


52.85 



U IL 

= 

71 3.64 

B/F T2HR 

DEG 

•JIB 


2152.97 

B/FT2HR 

DEG 

DELTAP 


8.866 

PSIA 


OPFR IC 

= 

6.192 

PSIA 


DPMOM 

- 

2.673 

PSIA 



SUPERHEAT REGION 


QSUP HS'JPS TSUPS TUBE L SOP OPSUP DPSFRIC DPSMOM 

BTU/SEC HR UNITS DEG E FT PSI PSI PSI 

3.21 71.08 1750.00 1.2619 2.362 2.155 0.2062 


Figure J>c . Typical Computer Printout of Overall Data From The 

Boiler Thermal Design Code (Case 1-F, Insert p/D = 2) 
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IV CONCLUDING REMARKS 


The analytical procedures and digital computer codes presented in this 
report are useful for thermal design studies of multiple-tube condensers and 
boilers applicable to Rankine cycle space power systems using potassium as the 
working fluid. The relatively short machine times of 3 to 7 seconds for each 
code enable rapid computation of several design-point cases to be done .with 
a minimum of cost and effort for parametric studies. 

Currently available experimental heat transfer results for condensing 
potassium (Reference 5) and once-through boiling of potassium in tubes with 
and without helical inserts (Reference 11) are the bases for the computer codes 
and associated analytical procedures. When more experimental data and improved 
analytical methods become available, the computer codes given in this report 
can be used as a guide for the preparation of more extensive programs. 
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APPENDIX A 


POTASSIUM CONDENSER THERMAL DESIGN CODE 
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TABLE I 


Program 

1) MAIN (CONDES) 

2) CONSEC 

3) DERIV 

4) TPRIMC 

5) CONDH 

6) SHELL 

7) PMULT 

8) GLANCE 

9) INTERP 

10) LININT 

11) THERMA 

a) DUBIN 

b) CURV 


PROGRAMS USED IN THE CONDENSER DESIGN CODE 


Primary Function 

Input and Output 

Calls derivative and integration subroutines. 
Stores computed data 

Computes derivatives of dependent variables 

Determines primary coolant temperature profile 

Computes the local condensing heat transfer 
coefficient 

Computes primary coolant heat transfer 
coefficient and shell side pressure drops 

Computes two-phase pressure drop multiplier 

Integration routine ( RKG + Lanczos) 

Performs double linear interpolation 

Performs linear interpolation 

NRL Thermodynamic properties 


TABLE II 


Symbol 
TITLE 
NT 
I TP 

NCF 

NQ 

MINC 

ITWMAX 

JHP 

OT 

XW 

DCB 

PODT 

WP 

GS 

POWER 

TSATIN 

TPIN 

TSUBO 

HP 

HCON 

*Symbols 


NOMENCLATURE FOR CONDENSER DESIGN CODE INPUT* 


Quantity 

Identification block, 48 characters 
Number of tubes 

Number of nodal iterations for primary coolant 
temperature profile 

Flow direction indicator 

0 = count erf low 

1 = parallel flow 

Number of nodes for primary coolant 
temperature profile 

Print interval indicator. May be 1, 2, 4 or 5 

Maximum number of iterations for average 
wall temperature 

Indicator for primary fluid heat transfer 
coef f ici ent 

0 = calculated by program 

1 = input 

Tube outside diameter 

Tube wall thickness 

Center body diameter 

Tube pitch to diameter ratio 

Primary coolant flow rate 

Secondary fluid mass velocity 

Total power 

Secondary inlet temperature 
Primary inlet 

Secondary subcooled outlet temperature 
Primary fluid heat transfer coefficient 
Secondary (condensing) heat transfer coefficient 


Units 


dimensionless 


inches 

inches 

inches 

dimensionless 

lb /sec 
m 2 

lb /ft -sec 
m 

KW 

°F 

°F 

°F 

Btu/ft 2 -hr-°F 

Btu/ft 2 -hr-°F 


are given in literal order of appearance on input cards 
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TABLE II (Cont 1 d) 


Symbol 

PF 


AK 

BK 

RADPL 

SLIPEX 


DELTAY 

PRINT 


TOLP 

TOLW 

ERLIMT 

XIN 

SIGMAC 

WGHT 


Quantity 

Primary fluid indicator 
-1 = NaK 

0 = Lithium 

1 = Sodium 


Intercept of linear equation for wall thermal 
conductivity = f(T) 

\ Slope of linear equation for wall thermal 
conductivity = f(T) 

Length of radial pressure drop region 


Reciprocal of slip ratio 


i 


exponent; viz, 2. = (^“) 

g 


Starting increment in quality 

Print indicator for local properties 

0 = Do not print 

1 = Print 

Tolerance on primary temperature iterations 
Tolerance on average wall temperature 


Maximum allowable per-step itegration 
error limit 

Secondary inlet vapor quality 

Condensation coefficient 0 O’ - 1 

c 

factors for integration errors, 
= 1 


Weighting 


4 

E 

i=l 




Units 


Btu-ft/ft^-hr-°F 

Btu-ft/ft 2 -hr-°F 

inches 

dimensionless 


dimensionless 


°F 

Percent 

100 

Percent 

100 

dimensionless 

dimensionless 

dimensionless 
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TABLE III 


NOMENCLATURE FOR CONDENSER DESIGN CODE PRINTOUT 


Symbol 

Quantity 

Units 

DELTAP 

Net condensing pressure drop 

psia 


DPFRIC 

Primary frictional pressure drop or condensing 
friction pressure drop 

psia 


DPMOM 

Condensing momentum pressure drop 

psia 

o 


GP 

Primary mass velocity 

lb /ft - 
m 2 

■sec 

GS 

Secondary mass velocity 

lb /ft - 
m 2 

■sec 

HP 

Local primary heat transfer coefficient 

Btu/ft - 

0 

1 

u 

A 

HSC 

Local condensing heat transfer coefficient 

Btu/ft 2 - 

S3* 

1 

o 

*1 

HSL 

Secondary heat transfer coefficient in 
liquid region 

Btu/ft 2 - 

■hr-°F 

NT 

Number of tubes 



PEL 

Secondary Peclet number in liquid region 

dimensionless 

POWER 

Total power 

KW 


PLOCAL 

Local secondary static pressure 

psia 


Q 

Cumulative heat transfer rate 

Btu/sec 


QCOND 

Heat transfer rate in condensing region 

Btu/sec 


QSUB 

Heat transfer rate in subcooled region 

Btu/sec 


QTOT 

Total heat transfer rate 

Btu/sec 


QOA 

Local heat flux 

Btu/ft 2 - 

•hr 

REL 

Secondary Reynolds number in liquid region 

dimensionless 

REP 

Primary Reynolds number 

dimensionless 

REZ 

Local secondary film Reynolds number 

dimensionless 

TPC 

Local primary fluid temperature 

°F 


TPIN 

Primary fluid inlet temperature 

°F 


TPOUT 

Primary fluid exit temperature 

°F 


TSC 

Local secondary fluid saturation temperature 

°F 


TSIN 

Secondary fluid inlet temperature 

°F 
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TABLE III (Cont'd) 


Symbol 

Quantity 

Units 

TSOUT 

Secondary fluid subcooled exit temperature 

°F 

UI 

Local overall heat transfer coefficient 

Btu/ft 2 -hr-°F 

UIL 

Overall heat transfer coefficient in liquid 
region 

Btu/ft 2 -hr-°F 

VMAX 

Inlet secondary vapor velocity 

fps 

VP 

Primary fluid velocity 

fps 

WP 

Primary flow rate 

lb /sec 
m 

ws 

Secondary flow rate 

lb /sec 
m 

X 

Local secondary 

dimensionless 

xw 

Wall thickness 

in. 

z 

Length 

ft . 



CONDENSER C O DE LISTING 

$JBF*TC CONDES FULIST.SDD 
C 

‘ DIMENSION T I TLE ( 8 ) 

C 

DIMENSION V I SL { ?0 ) »VISG(20).TT (20) *TCKL(20) 

DIMENSION WGHT ( 4 ) ,T PCM 100) ,QPRIM( 10 0 ) , X ( 1 00 0 ) * 

1TSCI 1000) ,Q( 1000) *HSC( 1000 ) , TPC ( 1 00 0 ) , U I ( 10 00 ) »QOA( 1000) ♦ 
2Z ( 1 000 ) 

DIMENSION HP ( 1000 ) .PLOCAL ( 1000 ) .REZ ( 1000 ) 

COMMON SHFA.SHID.GP , VP ♦ REP ♦ D ELPP . DP R AD ♦ DPCE X , DPTOTP 
COMMON NT » I TP » I TWMAX ,NCF ,NQ ♦ DT .XW.DIT.DLM.PODT » 
1WP.WS.GS.TSATIN.TPIN.P0WER.QT0T .HP * HCON ,TOLP * TOLW, 

2 PF » X I N » AK »BK .EXN.TPQ.QPRIM.X .TSC.HSC.TPC.UI .QOA. JL . 

3WGHT ,ERL IMT .DELTAX ,Q»Z 
COMMON SLIPEX 
COMMON P L OC AL ♦ R E Z * DEOV 
COMMON DCB.ZTOT 
COMMON SIGMAC 
COMMON JHP 
C 

C PROGRAM IS APPLICABLE TO A COUNTERFLOW OR 

C CO-FLOW CONDENSER WITH ENTERING QUALITY ; 

C XIN. SECONDARY FLUID EXITS SATURATED OR 

C SURCOOLED. 

C 

NAMELIST/INPUT/TITLE»NT»ITP*NCF,NQ»MINC» 

1ITWMAX»DT.XW»RADPL»DCB»P0DT»WP»GS*P0WER» 

2TSAT IN »TP I N * TSUBO ,HP»HCON»PF»XlN»TOLP»TOLW» 

3EXN * AK »RK » ERL I MT ♦ WGHT • PR I NT * DELTA Y ♦ SL I PEX » S I GMAC . JHP 

1 READ ( 5 » INPUT ) 

R I T S = 0 • 

DFLTAX = DELT AY 

RDP=RADPL 

ISURsO 

2 QT0T=.9A78*P0WER 
D I T = DT-2 • *X W 
DLM=2.*XW/AL0G ( DT/DI T ) 

C 

C COMPUTE WS FROM AN ENTHALPY BALANCE. THEN 

C COMPUTF GS OR NT FROM TOTAL POWER 

C AND SECONDARY FLOWRATE 

C 

2 TSATIR=T SATIN +4 59.69 
TSU*R0=TSUB0+459 . 69 

CALL THERMA ( PU »TSAT I R.PSIN.HLIN .HGI N.SL7 *SV7» 
1VL7.VGIN.CPG7.CPL7.EX.AV.-1. »0. .HS7 .SS7.VS7 .CPSH7 . EXS7 ♦ 
2AVS7) 

CALL THERMA ( PU *T SUB RO.P SO. HLOUT . HGOUT .SL7.SV7 *VL7 ♦ 
1VG7.CPG7.CPL7. EX »AV*-1« .0. »HS7 »SS7 * VS7 .CPSH7 ♦ EXS7 » A VS7 ) 
WS=QTOT / ( X I N^HGI N+ ( 1 . -X I N ) *HL I N-HLOUT ) 

I FIGS )4, 4 .5 

4 AS=3.1416*(DIT*#2-DCB*#2)/4.*FL0AT(NT)/144. 

GS-WS/AS 
GO TO 6 

6 ATUBE=3. 1416*( DI T **2-DCR^#2 ) /4. /144. 
xnt=ws/gs/atube 
NT= I F I X ( XNT ) 

6 CONTINUE 

WFT=3.1416*( DCB+DIT ) / 1 2 . *FLOAT ( NT ) 

DEQV=4.*AS/WET*12. 

C 
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COMPUTE THE PRIMARY TEMPERATURE PROFILE. 

CALL SL I TET ( 2 » JSEN ) 

CALL TPRIMC 

CALL SLITET ( 2, JSEN2 ) 

I F ( JSEN2-1 ) 1*1 ,666 
666 CONTINUE 

I F ( JHP >7,7,8 

8 SHFA=B ITS 
PODT^BITS 

sh i d=b its 

GP=B ITS 
VP=B ITS 
REP-BITS 
DPRAD=B ITS 
DELPP = B I T S 

9 CONTINUE 
7 CONTINUE 

ENTER SUBROUTINE CONSEC TO INTEGRATE THE 
DEOS IN THE CONDENSING REGION. 

SAVEX=DELTAX 

CALL SLITET ( 1* JSEN ) 

CALL CONSEC ( JL ) 

CALL SLITET ( 1* JSEN1 ) 

I F ( JSFN1-1 ) 1*1*11 
11 continue 

TEST TO SEE IF THE EXIT SATURATION 
TEMPERATURE IS LOWER THAN THE SUBCOOLED 
EXIT TEMPERATURE. 

T SATCE=T SC ( JL ) 

DTSUBC=TSATCE TSUBO 
10 IFfDTSUBC >15,16,16 

15 I F ( ISUB) 17, 17, 18 

18 WRITE (6,500) (TITLE(K) *K = 1 * 8 ) 

GO TO 1 
17 I SUB=1 

DELT AX = SAVEX 
GS=0 . 

T SUBO=TSUBO+DT SUBC-1 0 • 

GO TO 2 

16 CONTINUE 


COMPUTE QSUR AND THE LENGTH IN THE 
SUBCOOLED REGION. DETERMINE THE MEAN WALL 
TEMPERATURE BY SUCCESSIVE APPROXIMATION. 

T CWM 1=0. 

I TW= 1 

I F ( NCF ) 20 ,20,21 

20 TPP=TP I N 
TPOUT =TPQ ( 1 ) 

GO TO 22 

21 NZQ=NQ+1 
TPP=TPO ( NZQ ) 

TPOUT = TPP 
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2? CONTINUE 

TWO= ( TPP+TPC ( JL ) )/2« 

TWl = (TSC(JL)+T SUBO ) / 2 » 

TW2 = TW I +4 59 .69 

CALL THERMA ( PU»TW2 »PS7»HL7*HG7 *SL7* SV7 »VL7 » 
1VG7*CPG7,CPLS*EX*AV*-1.*0.»HS7*SS7*VS7»CPSH7»EXS7*AVS7) 

CALL LININT (TT »VlSL*TWl *VlSLS*OUT *2 0) 

CALL LININT (TT*TCKL*TW I *TCKLS*OUT *2 0) 

PRL=CPLS*VI SLS/TCKLS 
REL=DEQV/ 12 .*GS*3600 . /VISLS 
PECL=PRL*REL 
ZNUSC=3.66+.0055*PECL 
HTSC=ZNUSC* 1 2 • *T CKL S /DEQV 
QSC=QTOT-Q( JL) 

I F ( JHP ) 3 1 *31,32 

31 CALL SHELL ( 0 »TWO » TWO »RDP *HPSC » 1 ) 

GO TO 33 

32 HPSC=HP(1) 

33 CONTINUE 
T 1 =TW I 
T2=TWO 

23 TCWM2=AK+BK* ( TWO+TWl ) 

DELT A= ( T C WM2-T CWM 1 ) / T CWM2 

I F ( ABS (DELTA ) TOLW ) 100 , 100 , 24 

24 RESW=XW*DIT/12./TCWM2/DLM 

U I L= 1 • / ( 1 • /HTSC+RESW+D I T /HPSC/DT ) 

DTl-TSC( JL)-TPC( JL ) 

DT2 = TSUBO— TPP 

DTLMSC=(DT1-DT2) /ALOG ( DTI /DT2 ) 

ASC=OSC/UIL/DTLMSC*3600. 

DTI F=OSC* 3600* /HTSC/ASC 
DTOF=QSC*3600. /HPSC/ASC*DT /D IT 
T W I =T 1 -DT I F 
TWO=T2+DTOF 
T CWM1 =T CWM2 

I F ( ITW-ITWMAX) 25*25*26 
26 WRITE (6,1000) 

CALL EXIT 

25 I TW= I TW+ 1 
GO TO 23 


100 CONTINUE 

ZSC=ASC/3.1416*12./D IT /FLOAT (NT ) 

ZCON=Z ( JL ) 

ZTOT=ZSC+ZCON 

COMPUTE TOTAL CONDENSING PRESSURE DROP 
T69=TSATCE+459.69 

CALL THERMA ( PU * T 6 9 , PSOUT ,HL7 *HG7*SL7»SV7 ,VLOUT » 
1VG7,CPG7 ,CPL7* EX * AV »-l . *0. ,HS7 * S S 7 * VS7 , CPSH7 , EXS7 * AVS7 ) 
DELP- ( PS I N-PSOUT ) 

PMOM=GS**2/32. 2/ 144. * ( VLOUT-VG I N ) 

PFRIC=DELP-PMOM 
VMAX=GS*VG I N 
QCOND=Q ( JL ) 

DO 156 K= 1 * JL 
T SCRAN=TSC ( K ) 459.69 

156 CALL THERMA (PU*TSCRAN » P LOCAL (K) ,HL8»HG8»SL8»SV8»VL8» 

1VG8*CPG8»CPL8,EX*AV,-1. ,0 . , HS8 , SS8 » VS8 * CPSH8 » EXS8 , AVS8 ) 
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SET UP THE OUTPUT. 

TPIN1=TP I N 
I F ( PODT )1 09 ,109,111 

111 CALL SHELL ( 1 ,TPIN1 »TPOUT »RDP »HP1 ♦ 1 ) 

DPTOTP=DELPP+DPRAD 
GO TO 131 
109 DPT0TP=B ITS 
131 CONTINUE 

WRITE THE FIRST PAGE OF OUTPUT 
WRITE (6.5999) (TITLE(K) ,K = 1 ,8 ) .NT.ZCON.DIT ♦ 

1ZSC.DT .ZTOT .XW .PODT ♦ AS ♦ SH I D ♦ DCB ♦ SHF A .POWER » OTOT ♦ 

2QCOND ,QSC .TPIN.TSATlN.TPOUT.TSATCE.WP.TSUBO 
WRITE (6.6001) GP.WS.VP.GS.REP.VMAX.HTSC.REL.DELPP. 
1PECL .DPRAD.UIL .DELP ♦ DPTOTP .PFR I C ♦ PMOM 

IF(PRINT) 700,700.601 

601 K = JL+1 
DO 602 J = K » 1000 
X( J)=RITS 
T SC ( J ) =B I TS 
Q ( J ) =R I TS 
TPC ( J ) =B I TS 
U I ( J ) =R I TS 
HSC ( J ) =B I TS 
HP ( J ) =B I TS 
PLOCAL ( J ) =B I TS 
REZ ( J ) =B I TS 
QOA ( J ) =B I TS 

602 Z ( J ) =B I TS 
600 CONTINUE 

MINC=INTERVAL BETWEEN PRINT POINTS. 

= 1 ♦ 2 ♦ 4 ♦ 5 
MSTART= 1 
NPAGES=20/MINC 
I PAGE=0 
MF INAL=0 

DO 610 L= 1 » NPAGES 
I PAGE= I PAGE+1 
MFINAL=MF INAL MINC*50 
C 

WRITE (6,6000) I PAGE , ( T I TLE ( K ) , K* 1 , 8 ) , ( X ( K ) .Z ( K ) , 
1TSC(K) .PLOCAL (K ) , TPC ( K ) ,HSC( K) » HP ( K ) .U I ( K ) ♦ Q ( K ) ♦ 

2QOA(<) , REZ ( K ) ,K=MSTART»MF INAL ,M INC ) 

MSTART=MSTART 50*MINC 
I F( X (MSTART )-B ITS >610.700 ,610 
610 CONTINUE 
C 

700 CONTINUE 

DATA TT/ 6 00. ,700. ,800. ,900. » 1000. .1050. ,1100., 1150.. 

11200. . 1250. .1300. .1350. .1400 ..1450. .15 00. .1550. . 1600. » 

21650. . 1700. .1800./ 

DATA VISL/.58, .54 , . 496 ♦ . 449 ♦ .408 , . 3 93 , . 380 , .368* .356, 

1.346. . 336. .327.. 319. .312..305.. 298. .292. .287. .281. .271. 
DATA TCKL/25.3 ,24. 2, 23. 1,22. 1,21. 0,20. 5, 20. 0,19. 5,18. 9, 

118. 4, 17. 9, 17. 3, 16. 8, 16. 3, 15. 8, 15. 2, 14. 7, 14.1 ,13«6, 12. 6/ 
5777 GO TO 1 
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500 FORMAT GENERATOR 
RESTORE 
T I TLE=-X 
SPACE 2 

SUBCOOLED TOUT SPECIFIED TOO HIGH. 

END OF FORMAT 
1000 FORMAT GENERATOR 
RESTORE 

WALL TEMPERATURE IN SUBCOOLED REGION 
NON CONVERGENT. 

END OF FORMAT 
5999 FORMAT GENERATOR 
RESTORE 

CONDENSER DESIGN (GENERAL PARAMETERS) 


SPACE 
T ITLE=- 
SPACE 

X 




-A 





GEOMETRY 



SPACE 






-F4 FEET 

NT 

= 


-I 

COND LENGTH = 


SPACE 
TUBE ID 

rr 


-F4 INCHES 

SUBC LENGTH 


-F4 FEET 

SPACF 
TUBE OD 

«- 


-F4 INCHES 

TOTL LENGTH 


— F4 FEET 

SPACF 
WALL XW 



-F4 INCHES 

TUBE P/D 


-F2 

SPACE 







TUBE FLOW AREA= 


-F5 SQ FT 

SHELL ID 


-F3 INCHES 

SPACE 







CENTER BODY 0D= 


-F4 INCHES 

SHELL FLOW AREA= 


— F4 SQ FT 

SPACE ? 










THERMAL PARAMETERS 



SPACE- 
POWER KW 


OTOT 

* BTU/SEC 

QCOND, BTU/SEC 

OSUB 

♦BTU/SEC 

-F 2 
SPACE 



-F2 

-F2 


-F2 


PRIMARY 


SECONDARY 


SPACE 
TP I N 
SPACE 

= 

-F 2 

DEG F 

TSIN 

-F2 

DEG F 

TPOUT 

SPACE 


-F2 

DEG F 

TSATCE = 

-F2 

DEG F 

WP 

= 

-F4 

LBM/SEC 

TSOUT 

-F2 

DEG F 

END OF 

FORMAT 






FORMAT 

GENERATOR 





SPACE 






LBM/SEC 

GP 

= 

-F 3 

LBM/FT2SEC 

WS 

-F5 

SPACE 

VP 

SPACE 

= 

-F3 

FT/SEC 

GS 

-F3 

LBM/FT2SEC 


REP 


-FI 


VMAX - 

-F 3 

FT/SEC 

SPACE 




HSL 

-F2 

B/FT 2HRDEG 

SPACE 




REL 

-FI 


SPACE 

DPFRIC 

SPACF 

= 

-F3 

PSIA 

PEL 

-F2 


DPRAD 

= 

-F 3 

PSIA 

UIL 

-F2 

B/FT2HRDEG 

SPACF 




DELTAP = 

-F3 

PSIA 
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SPACF 

DPTOT = -F3 PSIA DPFRIC = -F3 PSIA 

SPACE ^ , 

DPMOM = -F 3 PSIA 

END OF FORMAT 
6000 FORMAT GENERATOR 
RESTORE 
SPACE 

CONDENSER DESIGN (LOCAL CONDITIONS) -I 


SPACE 

TITLE=-X 

SPACE 

X Z T SC PLOCAL 

X Q QOA REZ 

FT. DEG F PSIA 

X B/SEC B/FT2-HR 

-F5 -F4 -F2 -F 3 

X -F2 -FI -FI 

REPEAT 1 
END OF FORMAT 


-A 

TPC HSC HP UI 

DEG F HR UNITS HR UNITS HR UNITS 
-F2 -FI -FI -FI 


END 
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SIBFTC CONSFC FULIST.SDD 

SUBROUTINE CONSEC(JK) 

C SUBROUTINE INTEGRATES THE CONDENSER DEQ FOR 

C A SPECIFIED INLET VAPOR SATURATION TEMPERATURE. 

C THE FIRST 4 STEPS ARE RKG TO GENERATE 5 

C POINTS FOR THE LANCZOS SCHEME. 

DIMENSION VISL(20) » V I SG ( 20 ) , TT ( 20 ) ♦ TCKL ( 20 ) 

DIMENSION BITS(1 i 

DIMENSION WGHT (4) »TPQ( 100 ) *OPRIM( lOO).X(lOOO). 

1TSC( 1000) .0(1000) .HSC( 1000 ) .TPC ( 1 00 0 ) , U I ( 1 0 00 ) .QOA( 1000) ♦ 

2Z ( 1000 ) 

DIMENSION HP (1000) , PLOCAL ( 1000 ) .REZ(lOOO) 

COMMON SHFA.SHID.GP . VP ♦ REP , DELPP . DP RAD . DPCEX ♦ DPTOTP 
COMMON NT. I TP* ITWMAX.NCF .NQ.DT .XW.DIT .DLM.PODT ♦ 

1 WP , WS .GS .TSATIN.TPIN, POWER » QTOT .HP , HCON . TOL P » TOLW ♦ 

2PF.XIN.AK .BK.EXN .TPQ.QPRIM.X .TSC.HSC.TPC.UI .QOA. JL ♦ 

3 WGHT .ERLIMT .DELTAX.Q.Z 
COMMON SLIPEX 
COMMON PLOCAL.REZ.DEQV 
COMMON DCB.ZTOT 
COMMON SIGMAC 
COMMON JHP 

DIMENSION ZZ ( 3 ) .TSCC(3) .OCC(3) .XX (3 ) 

DIMENSION V A L ( 4 ) ,DVAL(4) .REFER (4) » E RR ( 4 ) 

DIMENSION DSTB ( 4 ) 

ERRMIN=ERLIMT/40. 

START=0 • 

6000 M0DE=4 

SAVE = DELT AX 
J = 1 

X( J)=XIN 
TSC ( J ) =TSAT I N 
Q( J) =0. 

Z(J)=0. 

1 K = 0 

TSC1 =TSC ( J ) 

X 1 =X ( J ) 

Q1=Q( J) 

5 CALL DERIV( Xl.Ql , TSC1 .VARJ.TPC1 .DVAL.HSC1 »U 1 1 .FLUXI .HPl ,REZ1 ) 

STARTING VALUES OF THE DERIVATIVES ARE NOW 
AVAILABLE FOR USE. 

I F ( MODE— 4 )80l,7,801 
7 VAL ( 1 ) =0 • 

VAL ( 2 ) =DVAL ( 1) 

VAL ( 3 ) =TSC1 
VAL ( 4 ) =Q1 
801 I F ( K ) 80 * 80 » 2 
80 HSC ( J ) =HSC1 
TPC ( J ) =TPC1 
UI (J)=UI1 
QOA ( J ) =FLUX I 
HP ( J ) =HP 1 
REZ ( J ) =REZ1 
K = 1 

2 CALL GLANCE (M0DE.X1.DELTAX.DVAL.VAL.REFER) 

I F ( MODF ) 3 .4 . 3 

S M0DE=-4 
VAR J=VAL ( 2 ) 

TSC1 =VAL ( 3 ) 


- 56 - 



nnrv r* n n r> nnn 


Q1 =VAL ( 4 ) 

GO TO 5 

C 

4 J=J 1 

2 ( J ) =VAL ( 1 ) 

VAR J=VAL ( 2 ) 

TSC ( J ) =VAL ( 3 ) 

Q( J)=VAL (4) 

X{ J)=X1 

IF( J-5)6,777,777 
6 M0DE=-4 
GO TO 1 
C 

777 I XA=4 

10 TSC1 =TSC ( J j 
X1=X ( J ) 

Q1=0( J) 

CALL DERIV(X1*Q1 , TSC 1 , VAR J , TPC 1 ,DVAL ,HSC1 ,U 1 1 *FLUXI ,HP1»REZ1 ) 
HSC( J)*HSC1 
TPC ( J ) =TPC1 
UI ( J )*UI1 
HP ( J j =HP 1 
REZ ( J) =REZ1 
QOA ( J ) =FLUX I 
IF(TSC1-TPC1 )6 09»609,21 
21 CONTINUE 

8 CALL GLANCE(M0DE»X1»DELTAX»DVAL,VAL»REFER) 

I F ( X 1 ) 340 ,340*222 

CHECK THE STABILITY OF THE SYSTEM 

22? XSTAB=X 1 

OSTAB=VAL (4) 

TSTAB=VAL ( 3 ) 

VSTAB=VAL (2 ) 

CALL DERIV( XST AB , QST AB , TSTAB » VSTAB * TPSTAB,DSTB, 
1HSTAB,UISTAB,FLXTAB»HPSTAB*REZTB) 

DO 273 K1 = 1 ,4 

IF(DELTAX*DSTB(K1 )/VAL(Kl ) +. 3 ) 30 » 27 3 ,2 73 
273 CONTINUE 

COMPUTE THE ERROR AND DOUBLE OR HALVE 
THE QUALITY INTERVAL ACCORDINGLY. 

11 DO 20 K= 1 ,4 

20 ERR(K) = .11*ABS( ( VAL ( K ) -REFER ( K ) )/VAL(K ) ) 

ERTOT =WGHT ( 1 )*ERR ( 1 ) +WGHT ( 2 ) *ERR ( 2 ) +WGHT ( 3 ) *ERR ( 3 ) 

1+WGHT ( 4 ) *ERR (4 ) 

I F ( ERTOT- ERL IMT)41, 30,30 
41 I F ( ERTOT-ERRMI N)40,40,50 

STEP SIZE IS HALVED IF ALLOWABLE. 

30 IF(MODE+l )32 ,197 ,33 
33 MODE=-l 
GO TO 8 


32 MODE=0 
I XA= I XA+1 
GO TO 50 


C 
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C STEP SIZE IS DOUBLED IF ALLOWABLE 

C 

40 I F ( I XA-4 ) 51 * 49 * 5 1 

49 MODE=-2 
I XA = 0 

GO TO 50 
51 M0DE=0 
I XA= IXA+1 

50 J=J 1 - 

Z( J)*VAL( 1 ) 

VARJ=VAL ( 2 ) 

TSC ( J ) =VAL ( 3 ) 

Q ( J ) =VAL ( 4 ) 

XTJ) =X1 
GO TO 10 
C 

C ... 

340 JK= J+l 

Z( JK)=VAL (1) 

VAR J=VAL ( 2 ) 

TSC ( JK ) =VAL ( 3 ) 

Q ( JK ) =VAL ( 4 ) 

X( JK)=X1 
C 

C TEST TO SEE IF THE FINAL QUALITY IS 

C ZERO OR NEGATIVE. IF 0, CONTINUE. IF 

C NEGATIVE. THE VALUES AT X=0,J=JL= ARE 

C DETERMINED BY INITIATING A RKG 

C INTEGRATION USING VALUES AT JL-1 WITH 

C A DELTAX OF X(JK-l). 

C 

I F ( X ( JK ) ) 105,106.106 
105 J=JK-1 

DELT AX =-X ( J ) 

M0DE=4 

555 TSC1=TSC ( J ) 

X1=X ( J ) 

Q1=Q ( J ) 

556 CALL DERI V( X1.Q1 ,TSC1 .VARJ.TPC1 .DVAL.HSC1 ,U II . FLUX I .HP 1 ♦ REZ 1 ) 
IF(MODE**4)557, 558.557 

558 HSC ( J ) =HSC1 
TPC ( J ) =TPC1 
UI( J)=UI1 
HP ( J ) =HPl 
REZ ( J ) =REZ1 
QOA ( J ) =FLUX I 
VAL ( 1 ) =Z ( J ) 

VAL ( 2 ) =DVAL ( 1 ) 

VAL ( 3 ) =TSC1 
VAL ( 4 ) =01 

557 CALL GLANCEIMODE, XI, DELTAX, DVAL, VAL. REFER) 

I F ( MODE ) 5 60 , 56 1 .560 

560 M0DE=-4 
VARJ=VAL ( 2 ) 

TSC1 =VAL ( 3 ) 

Q1 =VAL ( 4 ) 

GO TO 556 

C 

561 JK= J+l 

Z( JK)=VAL(1 ) 

VAR J=VAL ( 2 ) 
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TSC( JK)=VAL (3) 

Q ( JK ) =VAL ( 4 ) 

X( JK)=0. 

106 CALL DERIV( X(JK) » Q ( JK ) » TSC ( JK ) * VARJ »TPC (JK) »DVAL »HSC ( JK ) * 
1UI (JK) *QOA( JK) » HP ( JK ) ,REZ ( JK ) ) 

111 RETURN 


RESTART INTEGRATION. ALLOWABLE NUMBER OF 
REST ARTS= 5 • 

197 IF(START-5. ) 198,198,199 

198 ST ART = ST ART+1 • 

DELTAX=DELTAX/2. 

M0DE=4 

K = 0 

T SCI =TSC ( J ) 

X1=X ( J) 

Q1=Q(J) 

CALL DERI V( X1»Q1 » TSC 1 > VARJ » TPC 1 »DVA L »HSC1 »U 1 1 » FLUX I ♦ HP 1 , REZ 1 ) 
VAL ( 1 ) =Z ( J ) 

VAL ( 2 ) = DVAL ( 1) 

VAL ( 3 ) =TSC1 
VAL ( 4 ) =Q1 
GO TO 801 


199 WRITE (6,1002) ERL I MT » ERTOT , ( ERR ( K ) > K= 1 , 4 ) ♦ 

1 ( X (K ) *Z ( K ) » TSC ( K ) ,TPC(K)»K«1»J) 

CALL SLITE(l) 

GO TO 111 

609 WRITE (6,1001) J , ( X ( K ) , Z ( K ) , TSC ( K ) * TPC ( K ) , K = 1 , J ) 
CALL SL I T E ( 1 ) 

GO TO 111 


1001 FORMAT GENERATOR 
RESTORE 

J= -I 

SPACE 

TEMPERATURE CROSS HAS OCCURRED. 

SPACE 

X Z TSC TPC 

FT. DEG F DEG F 

-F5 -F5 -F2 -F2 

REPEAT 1 
END OF FORMAT 

1002 FORMAT GENERATOR 
RESTORE 

NUMBER OF RESTARTS EXCEEDED 
CHECK ERROR LIMITS. 

SPACE 

ERLIMT- -F5 ERTOT= -F5 ERR= -F5 -F5 -F5 -F5 

SPACE 

X Z TSC TPC 

FT. DEG F DEG F 

-F5 -F5 -F2 -F2 

REPEAT 1 
END OF FORMAT 
END 
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SIRFTC DERIV FuL 1ST *SDD 

SUBROUTINE DER I V ( X 1 * OX » TSC1 ♦ VAR J ♦ TP Cl » DVAR ♦HSCl » U 1 1 ♦ 
1FLUXI ,HP1 »REZ1 ) 

SUBROUTINE RETURNS VALUES OF THE DERIVATIVES AND 
THE LOCAL CONDITIONS TO SUBROUTINE SECOND. 

DIMENSION VISL(20) ,VISG(20) ♦ TT ( 20 ) * TCKL ( 20 ) 

DIMENSION BITS(1) 

DIMENSION WGHT < A ) »TPQ( 100 ) *QPRIM( 100)*X(1000)» 

1TSC( 1000) »Q ( 10 00 ) »HSC( 1000 ) »TPC(1000) »UI (1000 ) »QOA( 1000) » 
2Z ( 1000 ) 

DIMENSION HP (1000 ) »PLOCAL( 1000 ) »REZ ( 1000) 

COMMON SHFA.SHID»GP»VP»REP»DELPP»DPRAD»DPCEX,DPTOTP 
COMMON NT ♦IfP»lTWMAX»NCF*NQ»DT ♦ XW »D I T »DLM»PODT » 

1 WP » WS *GS »TSATIN»TPIN» POWER » QTOT » HP ♦ HCON ♦ TOL P * TOLW ♦ 
2PF.XIN.AK *BK »EXN * TPQ »QPR IM » X ,TSC ♦ HSC • TPC »U I »00A* JL» 
3WGHT.ERLIMT *DELTAX»Q,Z 
COMMON SLIPEX 
COMMON PLOCAL. REZ ,DEQV 
COMMON DCB»ZTOT 
COMMON SIGMAC 
COMMON JHP 

DIMENSION TEMP (20 ) , CPL I ( 20 ) * CPNA ( 20 ) ,CPNAK( 20 ) 

DIMENSION DVAR (4) 

NTQ=NQ+1 

CALL LIN I NT (QPRIM .TPQ.QX »TPC1 »OUT *NTQ) 

TP=TPC1 

COMPUTE DERIVATIVES 
S I G=S I GMAC 
WSS^WS 

TEM=TSCl+459.69 

CALL THERMA ( PU»TEM,PVAP1 *HL1 >HG 1 » SL 1 ♦ SV 1 » VL 1 ♦ VG 1 ♦ 

1CPG1 ,CPL1 * EX » ALPHA * -1. »0. *HS»SS*VSH *CPSH*EXS»AVS ) 

HFG1 =HG1-HL 1 
RHOL 1 =1 ♦ / VL 1 
16 IF(PF)22923»24 

22 CALL LININT (TEMP»CPNAK*TP»CPAV »OUT * 20 ) 

GO TO 25 

23 CALL LININT (TEMP»CPLI »TP»CPAV»OUT » 2 0 ) 

GO TO 25 

24 CALL LININT (TEMP»CPNA»TP»CPAV»OUT ♦ 2 0 ) 

25 RH0G1=1./VG1 

CALL LININT ( TT *VlSL *TSC1 »VISL1 »OUT*20 ) 

CALL LININT ( TT » V I SG » TSC 1 » V I SGI »OUT»20) 

PH11=PMULT(X1» RHOL1 » RHOG1 »VlSLl »VISG1»EXN) 

PVAP1=PVAP1 /14.696 

DPDT=l./( PV API *14. 696*144. *( 18717 .2 2/TEM**2-. 53299/TEM ) ) 

26 DTDP1=DPDT 
RELOC=DEQV*GS*300./VlSLl 
FRIC=.316/RELOC**.25 

GROUP=-l 2 .*PHIl*FRIC*GS**2/2./R HO Ll/32. 2/DEQV 
I F ( HCON ) 2 »2 ♦ 3 

2 CALL CONDH(TSCltWSS»Xl»SIG»HSCl*REVV*GROUP) 

REZ1 =REVV 
GO TO 53 

3 CONTINUE 
HSCl =HCON 

53 CONTINUE 
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I F { JHP ) 5 1 *51*52 

51 CONTINUE 
KC = 0 
RDP=0. 

CALL SHELL (0.TPC1*TPC1*RDP.HPP.KC) 

HP1 =HPP 
GO TO 57 

52 HPl =HP ( 1 ) 

57 CONTINUE 

DETERMINE LOCAL O/A AND LOCAL UI BY SUCCESSIVE 
APPROXIMATION 

TCWM1=0. 

ITW=1 
TWO=TPCl 
TWI =TSC1 

4 TCWM2=AK+BK* ( TWO+TWl) 

DELTA=(TCWM2-T CWM 1 ) /TCWM2 

I F ( ABS (DELTA ) TOLW)10»10»5 

5 RESW=XW*DIT/12./TCWM2/DLM 
UI1=1./(1./HSC1+RESW+DIT/HP1/DT ) 

FLUXI=UI1*( TSC1-TPC1 ) 

FLUXO=FLUXI*(DIT/DT ) 

TCWM1=TCWM2 

TWI=TSC1-FLUXI/HSC1 

TWO=TPCl+FLUXO/HPl 

I F ( ITW-ITWMAX) 6*6*7 
7 WRITE (6*1000) XI 
CALL EXIT 

6 I TW= I TW+1 
GO TO 4 

10 CONTINUE 

C1=WS*HFG1/3.1416/DIT*12./UI 1 /FLOAT (NT )*3600. 
IF(SLIPEX ) 6 1 *61*62 
61 SL I P= 1 • 

GO TO 63 
6? CONTINUE 

SL I P= ( RHOL1 /RHOG 1 ) ** ( 1 . /SLI PEX ) 

63 CONTINUE 

27 DVHAT= ( (l.+Xl*(SLIP-l») )*(1«/SLI P-RHOG 1 /RHOL1 ) + 
1 ( ( 1 .-XI ) *RH0G1 /R HOL 1+X1/SLIP )*( SLIP-1. ) )/RHOGl 

C2=WS*HFG 1 /WP/CPAV 

TEST TO SEE IF CONDENSER IS COUNTERFLOW OR 
COFLOW. 

I F ( NCF ) 29 *29 *2 8 

29 C2= C2 
GO TO 30 

28 CONTINUE 

30 A0FX=C1*GR0UP*DTDP1/(TSC1-TPC1)**2 
GROUP2=Cl/( TSC1-TPC1 )**2 
BOFX=GROUP2*( DVHAT*DTDP1+C2 ) 

I F ( X 1-X I N ) 40 ,41,40 
41 DLDX=-C1/ (TSC1-TPC1 ) 

GO TO 42 
40 DLDX=VARJ 
4? DDLDX=AOFX*DLDX+BOFX 

DTSDX=GROUP*DTDPl*DLDX-GS**2/32 .2*DVHAT*DTDP1 
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DQDX=-WS*HFG1 
DVAR ( 1 ) =DLDX 
DVAR ( 2 ) =DDLDX 
DVAR ( 3 ) =DTSDX 
DVAR ( 4 ) =DQDX 
?000 RETURN 

DATA TT/6 00. ,700. >800. ♦ 900. » 1000. *1 050. ♦1100. *1150. » 
11200.»1250.*1300.»1350.»1400.»1450.*1500.*1550.»1600** 
21650. ♦1700. *1800./ 

DATA VlSL/.58».54».496».449».408».393».380*.368».356» 
1.346*.336*.327*.319*.312»*305».298» .292 *.287 *.281 *.271/ 
DATA V I SG/ . 0349 ♦ .0365 » *0379 » .0395 * • 041 1 * .0419 » .0427 » 

1. 0435*. 0442 *.0449». 0457*. 046 4*. 0471*. 0478*. 04 85». 0493* 
2. 0499*. 0506, .0 512 *.0524/ 

DATA TCKL/2 5 .3 *24. 2*23. 1*22. 1*21.0*20. 5*20. 0,19.5*18. 9* 
118. 4*17. 9*17. 3*16. 8*16. 3*15. 8*15. 2*14. 7*14. 1*13. 6*12. 6/ 
DATA TEMP/400. *500 ., 550 . *600 .» 650 .» 700. » 

1750. *800. *850. *900. *950. *1000. *1050. *1100. *1150. * 1200. * 

21300. . 1400. .1500. .1600./ 

DATA CPNAK/.2166 , .2 138 , .2127 ♦ .2118* .21 10 *.2 103 , 

1.2097*. 20 93 *.2090*. 2088*. 2088*. 2089*. 2090*. 2093*. 2096* 
2. 210 *.2 11 *.2 122*. 2 139 *.2 159/ 

DATA CPNA/. 320*. 316*. 314*. 312*. 310*. 308, 

1.306. . 304 ..3032*. 3022*. 30 14*. 30 08 , .3004, .3*. 3* 

2. 3 *.3006* .3018 ,.3 03 8*. 3064/ 

DATA CPU/. 992 *. 992 *. 992 , .992*. 992*. 992 , .992, 
1.992».992».992*.992».992*.992*.992*.992».992».992» 

2. 992*. 992*. 992/ 

1000 FORMAT GENERATOR 
RESTORE 

WALL TEMPERATURE NON CONVERGENT. 

X= -F 5 

END OF FORMAT 
END 
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BFTC SUBR0l1t INE^CONDH? TLOC »WS *X1 » SI GMAC » HCON » RES »G ROUP ) 

SUBROUTINE RETURNS LOCAL CONDENSING HEAT 
TRANSFER COEFFICIENT BASED ON THE SUM 
OF A VAPOR PHASE AND A LIQUID FILM 
RESISTANCE 

DIMENSION V I SL ( 2 0 ) » V I SG ( 20 ) * TT (20) » TCKL ( 20 ) 

DIMENSION BITS(I) 

DIMENSION WGHT ( 4 ) »TPQ ( 100 ) »QPR I M ( 1 0 0 ) » X ( 100 0 ) ♦ 

1TSC( 1000) ,Q( 1000 ) ,HSC( 1000) , TPC ( 1 00 0 ) , U I ( 1 0 00 ) *QOA( 1000) » 
2Z ( 1000) 

DIMENSION HP ( 1000 ) » P LOCAL (1000) >REZ(1000) 

COMMON SHFA»SHID»GP»VP»REP*DELPP»DPRAD*DPCEX»DPTOTP 
COMMON NT,ITP» I T WMAX » NCF »NQ » DT » XW » D I T » DLM » P OD T » 

1 WP ♦ WS »GS »TSATIN»TPIN» POWER »QTOT ,HP » HCON * TOLP ♦ TOLW * 
2PF,XIN,AK,BK,EXN»TPQ,QPRIM*X,TSC,HSC,TPC»UI »qoa » jl * 

3WGHT »ERLIMT »DELTAX»Q*Z 
COMMON SLIPEX 
COMMON PLOCAL » RFZ »DEQV 
COMMON DCB »ZTOT 
COMMON SI GMAC 
COMMON JHP 

DIMENSION REF ( 10 ) »T AU ( 6 ) ,CONRL (10*6) 

DIMENSION T ( 20 ) ,THCONK(20) ,VISK(20) 

1 TRANK=TLOC+459.69 

CALL THERMAIPU. TRANK, PVAP1,HL1.HG1,SL1.SV1,VL1,VG1, 

1CPG1 »CPL1 » EX , ALPHA »-l • »0. ,HS »SS »VSH »CPSH »EXS * AVS ) 

2 HFG1 =HG1-HL 1 
RHOL 1 = 1 • / VL 1 
RHOG1 =1 • / VG1 

CALL LININT ( T » THCONK » TLOC »THC0N1 » OUT ,2 0) 

CALL LININT ( T , VI SK , TLOC , VI SKI ,OUT ,2 0) 

3 PRES S =144 »*PVAPl 

4 DPDT=PRESS*( 1871 7.22 /TRANK**2-. 5329 9/TRANK) 

5 GROVP=39. 1*32. 2/2. /3. 1416/15 45. /TRANK 
GR0UP1=SQRT (GROVP ) 

HV=3600.*SIGMAC*HFG1*GROUP1* (DPDT-P RES S/2./ TRANK) 

COMPUTE LIQUID FILM RESISTANCE FROM THE 
DUKLER FILM THICKNESS 

IF ( XI— XIN ) 6 , 7 , 6 

7 RES =100. 

REST = 0. 

TVST AR = 500. 

RATI 0=32 .2*3600. **2 *RHOL 1 **2 /VISK1**2 
GO TO 801 

6 DPDL=GROUP 

FTERM=(RH0G1-DPDL ) /RHOL1 
RAD=D I T / 2 4. 

RAT 1 0=32 . 2*360 0.**2*RHOL1**2/VI SKI* *2 
TVSTAR=FTERM*RAD/2.*RATIO**( 1./3. ) 

I F ( TVSTAR-500. >20,20,30 
30 TVST AR = 500. 

20 CONTINUE 

TVSL=ALOG10(TVSTAR) 

GAMMA=WS* ( 1 • -X 1 ) *43 2 00. / 3. 14 16 /D IT/ FLOAT (NT ) 

8 RES=4.*GAMMA/VISK1 
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REST = RES 

I F ( RES-100. ) 35*21 *21 
35 RES=100. 

21 CONTINUE 
801 CONTINUE 

REZ1 =ALOG 1 0 ( RES ) 

CALL INTERP( REF , T AU » CONRL ♦ REZ 1 *TVSL * CONRU ♦ 1 0 ♦ 6 ♦ ERR ) 
C0NRL1=EXP( 2.303*CONRU) 

HL=THC0N1*RATI0** ( 1./3. )*C0NRL1 
HC0N = 1 • / ( l./HV+l ./HL ) 

RES=REST 
500 RETURN 

DATA BITS/0-377777777777/ 

DATA T/600.* 700. *800. *900. *1000. *10 50. *1100. *1150.* 

11200. *1250. *1300. *1350. ,1400 ., 1450., 1500. *1 550. *1600. » 

21650. . 1700. .1800./ 

DATA VI SK/. 58*. 54*. 496* .449* .40 8*. 393*. 3 80 *.368*. 356* 
1.346*.336*.327,.319,.312*.305».29 8, .292 ♦ .2 87, .281 ,.271/ 
DATA THCONK/25.3 *24.2 *23. 1 *22. 1 ,21. 0,20. 5*20. 0,19. 5*18. 9* 
118.4,17.9,17.3,16.8,16.3,15.8,15.2,14.7,14.1,13*6,12.6/ 
DATA REF/ 2. ,2.30103 ,2.47712,2.77815 ,2.90309,3.0,3.30103, 
13.47712,3.8451 ,4.0/ 

DATA TAU/-10., 1.30103,1.69897,2 .00, 2.30103 ,2.69897/ 

DATA CONRL /-. 6903 7, -.72 81 6,-. 882 73, -.91 722 ,-.96257, 

1- 1.00044,-1.13077,-1.20761,-1.38722 , -1 . 47496 ,-. 2 0066 , 

2- . 33724 ,-.42 59 7, -.59346 ,-.67 16 2 ,-.73283 ,-.92082,-1.03621, 

3- 1.28819,-1 .39794, .00432,-. 14267, -.2403 3 ,-.42022 ,-.49757 , 

4- . 56067, -.76 700, -.8 8606, -1.1 5490,-1 .27572,. 14922 ,-.00437, 

5- .10237, -.28 400, -.36653, -.43180 ,-.64782 ,-.76955,-1.05061 , 

6- 1.16749, .25285 , . 1 30 33 , . 0453 2 , - . 1 30 77 , - . 207 61 ,-.27165, 

7- .48812 *-.55284,-. 92 082 ,-1.05061 » . 5 3 148 , . 3 8 20 2 , . 2 8 330 , 

8.09691 .. 01284, -.05552 ,-.27984,-. 41454, -.71670 ,-.84771/ 

END 
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I F ( JK ) =0 COMPUTES HPi 

I F ( JK=1 COMPUTES FRICTION AND RADIAL 

PRESSURE DROP 

DIMENSION REG ( 5 ) » EPSMAX ( 5 ) 

DIMENSION TEMP (20) ,CPNAK ( 20 ) »CPLl (20) ,CPNA ( 20 ) * 

1RHONAK (20)»RHOLI(20) , RHONA (20) » VISNAK ( 20 ) ♦VISLI (20) « 

2V I SNA ( 20 ) ,T CNAK ( 20 ) , TCL I (20) »TCNA(20) 

DIMENSION VISL(20)»VISG(20)»TT(20)*TCKL(20) 

DIMENSION B I TS ( 1 ) 

DIMENSION WGHT ( 4 ) »TPQ( 100) »QPRIM( 100) »X(1000) » 

1 TSC ( 1000 ) »Q ( 1000 ) ♦ HSC ( 1000 ) ♦ TPC ( 1000 ) *U I ( 1000 ) *Q0A( 1000) ♦ 

2Z( 1000) 

DIMENSION HP(1000) , PLOCAL ( 1000 ) ,REZ(1000) 

COMMON SHFA»SHID»GP*VP,REP,DELPP*DPRAD»DPCEX»DPTOTP 
COMMON NT » I TP » ITWMAX »NCF »NQ * DT * XW »D I T *DLM»PODT , 
lWP,WS»GS»TSATIN»TPIN»POWER»QtOT»HP* HCON »TOLP»TOLW» 
2PF»XIN,AK»BK»EXN»TPQ,QPRIM»X»TSC*HSC»TPC»UI »QOA * JL ♦ 

3WGHT ,ERLIMT ,DELTAX,Q»Z 
COMMON SLIPEX 
COMMON PLOCAL, REZ »DEQV 
COMMON DCB »ZTOT 
COMMON SIGMAC 
COMMON JHP 


INITIALLY COMPUTE PROPERTIES 


1 IF(PF)2»3»4 

2 CALL LINI NT ( TEMP »CPNAK *TEMP1 »CP1 »OUT»20 ) 

CALL LININT ( TEMP ,CPNAK * TEMP 2 >CP2*OUT»20) 

CALL LININT ( TE MP ♦ RHONAK » TEMP 1 » RHO 1 * OUT » 20 ) 

CALL LININT ( TEMP * RHONAK , TEMP2 » RH02 * OUT » 20 ) 

CALL LININT ( TEMP » V I SNAK » TEMP 1 • V I S 1 » OUT » 20 ) 

CALL LININT ( TEMP , V I SNAK * TEMP2 » V I S2 * OUT • 20 ) 

CALL LININT ( TE MP « TCNAK » TEMP 1 » TC 1 » OU T » 2 0 ) 

CALL LININT ( TEMP » TCNAK ♦ TEMP2 ♦ TC2 » OUT *2 0 > 

GO TO 5 

3 CALL LININT (TEMP»CPLI »TEMP1»CP1 »OUT *20) 

CALL LININT ( TEMP »C PL I »TEMP2 »CP2 »OUT *20) 

CALL LININT ( TEMP ♦RHOLI » TEMPI *RH01 »OUT»20) 

CALL LININT ( TEMP, RHOLI ,TEMP2 ,RH02 ,0UT,20) 

CALL LININT (TEMP, VISLI , TEMPI »VIS1 *OUT ,20) 

CALL LININT ( TEMP, VISLI *TEMP2 ,V I S2 ,OUT ,20 ) 

CALL LININT ( TEMP , TCL I , TEMPI , TC 1 , OUT , 20 ) 

CALL LIN I NT (TEMP, TCL I , T EMP2 , TC2 ,OUT ,20) 

GO TO 5 

4 CALL LININT ( TEMP ,CPNA, TEMPI ,VIS1 , OUT, 20 ) 

CALL LININT ( TE MP , CPN A , TEMP2 , V I S2 » OU T , 2 0 ) 

CALL LININT ( TEMP , RHONA , TEMP 1 , RH01 ,OUT,20) 

CALL LINI NT ( TEMP , RHONA , TEMP 2 ,RH02 ,OUT ,20 ) 

CALL LININT ( TEMP , V I SNA , TEMP 1 *V I Si »OUT *20 ) 

CALL LININT ( TEMP ,V I SNA , TEMP 2 »V I S2 »OUT ,20) 

CALL LININT ( TEMP , TCN A , TEMPI , TC 1 ,OUT ,20 ) 

CALL LININT(TEMP,TCNA,TEMP2*TC2»OUT ,20) 

5 CONTINUE 
IF(KC)6»6»7 

6 SHID=DT*PODT*SQRT(2.*SQRT(3. ) / 3 . 1 41 6*FL0AT ( NT ) ) 
SHFA=FLOAT(NT)*3.1416*DT**2/4.*(2.*SQRT (3. ) / 
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1 3 « 1416*P0DT **2-1 • )/144. 

NT1 = 1 
1=1 

4111 NTl=NTl+6*( 1-1 ) 

I F ( NT1-NT )4112*4113»4114 
411? 1=1 1 

GO TO 4111 
4114 1=1 1 
411"? CONTINUE 

P I TCH=PODt*DT 

RATlOsPITCH/SQRTfSHl D**2-P I TCH#*2 ) 
ANG*ATAN( RATIO ) 

NOUTR=6*( I-i ) 

XNROWS= FLOAT ( NOUTR ) 

XNTl=FLOAT(NT) 

WETPER=3.14l6* (SHID+XNT1*DT ) +XNROWS* ( 2 . * 

IP ITCH /SORT ( 3 • ) -SH I D*ANG ) 
DEQSH=4.*SHFA*12./WETPER 
GP=WP/SHF A 
KC = 1 

7 I F ( JK ) 8 ♦ 8 »8 0 

COMPUTE LOCAL PRIMARY HEAT TRANSFER 
COEFFICIENT. 

8 RE1=DEQSH*3600.*GP/VIS1 j 

PR1=CP1*V IS1/TC1 , 

PEC1 =RE1*PR 1 

REPL=ALOG 10 ( RE 1 ) ; 

CALL LININT (REG»EPSMAX»REPL» EPS ♦ OUT »5 ) 
EPSP=EXP(2.303*EPS) / 

PSIBAR=1.-1 .82/PRl/EPSP**1.4 
ALPHA = • 93 + 1 0 • 8 l*PODT-2 » 01*PODT**2 
BETA=.0252*P0DT**.273 
GAMMA=. 8 

I F (PS I BAR ) 666 > 666*667 

666 ZNUP=ALPHA 
GO TO 668 

667 CONTINUE 

ZNUP=ALPHA+RETA*( PSIBAR*PEC1 )**GAMMA 

668 CONTINUE 
HP1 =ZNUP*TC 1 /DEQSH 

9 RETURN 

COMPUTE FRICTION PRESSURE DROP BASED 
ON AVERAGE SHELL TEMPERATURE. 

80 TPAV=(TEMP1+TEMP2 ) /2. 

I F ( PF ) 81 ♦ 82 » 83 

81 CALL LININT(TEMP»RHONAK»TPAV*RHOAV»OUT»20) 
CALL LININT (TEMP*VISNAK»TPAV ♦ V I SAV » OUT »20) 
GO TO 85 

8? CALL LININT(TEMP,RHOLI.TPAV.RHOAV.OUT,20) 
CALL LININT ( TEMP » VI SL I »TPAV»VISAV »OUT »20 ) 

- GO TO 85 

83 CALL LININT (TEMP»RHONA»TPAV*RHOAV»OUT *20) 
CALL LININT ( TEMP » V I SNA » TPAV » VI SAV »0UT»20) 
85 CONTINUE 

REP=DEQSH*3600.*GP/VISAV 
VP=GP/RHOAV 
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FRIC=.316/REP**.25 

DELPP=FR I C*ZTOT*GP**2/2 . /RHOAV/32 .2 /DEOSH 
1/144. 

IF(RADPL) 500 .500 ♦ 501 
COMPUTE RADIAL PRESSURE DROP 

501 F=0 • 

XNT=FLOAT (NT) 

DO 502 J = 2 » I 
NROW= J-l 

F1=FL0AT (NROWl/XNT 
F=F Fl**1.84 

502 CONTINUE 

TERM 1=. 25 /( PODT-1. ) **1.84 
TERM2=.118/(PODT-l. )**2.92 
GAM = ( TERM1+TERM2 ) *F 

TERMI=(VIS1*RADPL/12 ./WP/3600. )**.16 
TERMO=( VIS2*RADPL/12 ./WP/3600. ) **.l 6 
TERM3=(WP/RADPL*144./DT)**2 
C 

DPRAD=(TERMI /RH01+TERM0/RH02 )*. 56 2/32. 2/1 44. 

1*GAM*TERM3 
GO TO 510 
500 DPRAD=0. 

510 CONTINUE 

DPTOTP=DELPP+DPRAD 
1067 RETURN 
C 

DATA REG/3. *4. »5. *6. »7./ 

DATA EPSMAX/ .477 ♦ 1 .2455 ♦ 2.072 »2 . 9 77 5.3.929/ 

DATA TEMP/400. >500. *550. *600. >6 50. ♦700.»750.»8 00. ♦ 

1850. >900. >950. >1000. >1050. >1100. >11 50. >1200. >1300. > 

21400. . 1500. >1600./ 

C 

DATA CPNAK/. 2166 > .2 138 >.2126 >.2118 > . 2 1 1 0 > . 2 103 > . 20975 > 

I . 2093 >.20 90, .2 08 8 >.2088 > .2089 >.2090 >.2 0925 > .2 096 >.2100 > 

2. 21 09 >.2122 >.2138 >.2158/ 

C 

DATA CPL 1/1.0115 >1.0099 >1.00 91 >1.00 83 >1.007 5 >1.0067 >1.0059, 

II. 0051 >1.0043, 1.0035 >1.0027, 1.0019, 1.0011 >1 .0003,0.9995, 

20. 9987.0. 9971 >0.9955 >0.9939,0. 9 92 3/ 

C 

DATA CPNA/. 3 25 2 > .3175>.3142>.3114>. 3090>.3070>.3055> 

1 .3042 >.30 32 > .3022 > . 301 5 > . 300 8 > . 3004 > . 3 000 > . 2999 > .3000 > 

2 .3006 >.3018 > .3038 > .3064/ 

C 

DATA RHONAK / 51.42 ,50. 58,50.16,49. 74, 49.32,48.90, 

148. 48,48. 06,47. 64,47. 22,46. 80,46. 38 >45. 96,45. 54,45.12, 
244.70,43.84,43.00 >42.16,41.32/ 

C 

DATA RHOL 1/3 1.27 >30. 98 >30. 85 >30.70 >30. 56 >30.42 >3 0.28, 

130. 15, 30. 01, 29. 87, 29. 74 >29. 60,29. 48 >29. 35,29.23, 
229.10,28.85,28.60,28.36,28.12/ 

DATA RHONA/56. 48 , 55 . 50 ,55. 01 ,54.54, 54.07 >53 .60 >53. 14, 

152. 65, 52. 21, 51. 78, 51. 35, 50. 92, 50. 50 >50. 11 >49. 71, 49. 33, 
248.59,47.86,47.17,46.47/ 

C 

DATA VI SNAK/. 8 95, .780, .730, .68 8, .64 8, .612, .580, .551, 

1 .523 >.500 ,.476 >.4 54, .435 ,.418, .402, .388 ,.362, .342 , 
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2. 323*. 308/ 

C 

DATA VI SL 1/1.37*1. 22*1. 16*1. 10 *1*055*1. 01 **973** 935* 

1 #900 * #866 * • 8 AO ♦ • 812 * #786 * #765 * • 744 * #725 * #686* #652 ♦ 

2. 620*. 592/ 

C 

DATA VI SNA/. 860*. 805 *.777* .751 » .725 *.699 *.6 74 » .648 ♦ 

1.624 *.600 ,.576 *.5 56 *.53 3 *.513 *.493* .47 6 *.444* .416, 

2.392 ,.372/ 

C 

DATA TCNAK/14. 34*14.82*14.93*15.01* 15.067 *15.115*15.146* 
115. 167 *15. 173* 15. 160, 15. 135*15. 096* 15. 040 *14. 985 *14. 92 5* 
2 14. 86 *14. 732 *14. 5 87 *14. 432 *14. 266/ 

C 

DATA TCLI/26.62»26.775»26.855»26. 928*27.005 *27.08, 

127. 16. 27. 24, 27. 32, 27. 395, 27. 47 *27. 55, 27. 63 *27. 705* 
227.78*27.86,28.01 *28.16*28.32*28.48/ 

C 

DATA TCNA /48 .90*46. 25*45.04*43. 97 ,42.95*42.02 *41.20* 

140. 42 *39. 63 *38. 95, 38. 25, 37. 63*37. 06, 36. 53*36. 01 *35. 53* 
234.60,33.72*32.94*32.23/ 

END 
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$IBFTC TPRIMC FuLlST.SDD 
SUBROUTINE TPRIMC 

SUBROUTINE EVALUATES PRIMARY TEMPERATURE 
PROFILE AS A FUNCTION OF Q TRANSFERRED. 

DIMENSION VISL(20)*VISG(20)»TT(20)»TCKL(20) 

DIMENSION BITS(l) 

DIMENSION WGHT (4 ) *TPQ( 100 ) »QPRIM( 100 ) »X ( 1000 ) ♦ 
lTSC(lOOO) *Q C 1000 ) »HSC ( 1000 ) . TPC ( 100 0 ) *U I ( 10 00 ) »Q0A(1000) ♦ 
2Z ( 1000 ) 

DIMENSION HP ( 1 000 ) »P LOCAL (1000) »REZ ( 1000 ) 

COMMON shfa,shid.gp.vp»rep,delpp,dprad.dpcex,dptotp 
COMMON NT * I TP * ITWMAX »NCF »NQ > DT > XW »D I T .DLM.PODT ♦ 
1WP.WS.GS.TSATIN.TPIN » POWER »QTOT .HP ♦ HCON • TOLP ♦ TOLW » 
2PF.XIN.AK.BK.EXN.TPQ.QPRIM.X.TSC.HSC.TPC.UI »QOA » JL * 

3WGHT .ERLIMT .DELTAX.Q.Z 
COMMON SLIPEX 
COMMON PLOCAL.REZ.DEOV 
COMMON DCB » Z TOT 
COMMON SIGMAC 
COMMON JHP 

DIMENSION TEMP (20)»CPLI(20)»CPNA(20) »CPNAK ( 20 ) 

DELTAO=QTOT /FLOAT ( NQ ) 

TEST TO SEE IF PARALLEL OR COUNTERFLOW. 

IF(NCF)100»100»200 
100 K5=NQ+1 

TPQ( K5 ) =TP I N 
QPR I M ( K5 ) =QTOT 

50 TP=TPQ(K5) 

I T 1 = 1 

6 IF(PF)2»3»4 

2 CALL LININT (TEMP»CPNAK»TP»CPAV»OUT .20) 

GO TO 5 

S CALL LININT (TEMP.CPLI »TP»CPAV»OUT *20) 

GO TO 5 

4 CALL LININT ( TEMP ,CPNA»TP »CPAV»OUT *20) 

5 TP1=TPQ(K5)+DELTAQ/WP/CPAV 
TZ= ( TP 1 + T PQ ( K5 ) ) / 2 • 

I F ( ABS ( T Z-T P ) TOLP ) 7 *7 * 8 

8 IFC IT1-ITP) 10.11.11 
11 WRITE (6.1000) K5 
CALL SL I T E ( 2 ) 

GO TO 500 
10 I T1 = I Tl + 1 
TP = TZ 
GO TO 6 

7 K5=K5-1 
TPQ(K5)=TP1 

QPR I M ( K5 ) =QPRIM(K5+1 l-DELTAQ 
I F ( K5-1 ) 5 1 » 5 1 ♦ 50 

51 QPR I M ( 1 ) =0. 

GO TO 500 

PARALLEL FLOW 

200 K6=NQ+1 
K5 = l 

TPQ( K5 )=TPIN 
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QPRIM(K5) =0. 

60 TP=TPQ(K) 

I Tl = l 

16 IF(PF >22*23*24 

22 CALL LININT ( TEMP ♦ CPN AK , TP »CP AV ,0UT » 20 ) 
GO TO 25 

23 CALL LININT (TEMP»CPLI »TP»CPAV»OUT *20) 
GO TO 25 

24 CALL LININT (TEMP»CPNA*TP»CPAV*OUT *20) 

25 TP1=TPQ(K5)+DELTAQ/WP/CPAV 
TZ= ( TP1 + TPQ ( K5 ) )/2. 

I F ( ABS ( TZ-T P ) TOLP >27*27*28 
28 I F ( IT1-ITP)30»31 »3l 
31 WRITE (6*1000) K5 
CALL SL I T E ( 2 ) 

GO TO 500 
30 IT1=IT1+1 
TP=TZ 
GO TO 16 
27 K5=K5+1 

TP0(K5)=TP1 

QPR I M { K5 ) =OPRI M ( K5-1 ) +DELTAQ 
I F ( K5-K6 )60»61*61 

61 QPR I M ( K6 ) =QTOT 
GO TO 500 


DATA TEMP/400. >500. *550. ,600. *650. *700. » 

1750. *800. *850. *900. *950. *1000. *1050. *11 00. *11 50 .*1200.* 
21300. *1400. *1500. *1600./ 

DATA CPNAK/.2166 *.2138*. 2127, .2118*. 21 10*. 2103* 

1.2097*. 20 93, .2 090 *.2 088, .208 8, .2089 *.2 090* .20 93* .2096* 

2. 210 *.211 *.2122, .21 39*. 2159/ 

DATA CPNA/. 320 *.316* .314 *.3 12 *.310, .308 * 

1.306 *.304 *.3032* .3022 *.3014* .3008*. 3004,.3» .3 * 

2. 3 ,.3006, .30 18 *.303 8*. 3064/ 


DATA CPLI /.992 *.992 * .992 » .992 * . 992 , .992 » .99 2 » 

1.992 *.992*. 9 92 *.9 92 *.992*. 992*. 99 2, .99 2*. 992*. 99 2* 
2. 992*. 992*. 992/ 

500 RETURN 

1000 FORMAT GENERATOR 
RESTORE 

PRIMARY TEMPERATURE NON CONVERGENT 
SPACE 5 
K5=- 1 

END OF FORMAT 
END 
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IBFTC PMULT 

FUNCTION 


FUL I ST » SDD 

RETURNS 2-PHASE MULTIPLIER 


FUNCTION PMULT ( X * RHOL ♦ RHOG » V I SL »V ISG»EX ) 
1 A=VI SL /V I SG-1 • 

B=RHOL/RHOG-l. 

C=1 • / ( 1 «+X*A ) 

D=1 .+X*B 

PMULT =C**EX*D 
50 RETURN 
END 
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$ibftc glance fulist,sdd 

SUBROUT INE GLANCE(MG»TF»HG»DG*YG»ZG) 

MG= MODE INDICATOR 
TF= INDEPENDENT VARIABLE 
HG= STEP SIZE 

DG ( I ) = DERIVATIVE OF YG EVALUATED AT TF. 

YG { I ) — DEPENDENT VARIABLE. INTEGRATED VALUE OF DG ( I ) AT T. 

Z G ( I ) = R E F • VALUE OF YG ( I ) FOR ESTIMATE OF TRUNCATION ERROR 
I MAX LESS THAN OR EQUAL TO 30. 

INSTRUCTIONS TO GLANCE ARE MADE AS FOLLOWS 
M= N START INTEGRATION WITH N DEPENDENT VARIABLES 
M-0 CONTINUE, NO CHANGE 
M= 1 HALVE HG REPEAT LAST STEP 
M= 2 DOUBLE H , TAKE NEXT STEP 
M= 4 CONTINUE USING RUNGE KUTTA GILL METHOD 

GLANCE INDICATES THE ACTION TAKEN AS FOLLOWS 
M=0 STEP WAS COMPLETED WITHOUT CHANGE 
M= 1 STEP WAS REPEATED WITH HG HALVED 
M= 2 STEP WAS TAKEN WITH HG DOUBLED 
M- 3 STEP WAS PARTIALLY COMPLETED R-K-G. 

M= 5 ERROR, M WAS NOT POSITIVE ON FIRST CALL 

DIMENSION YG(30) ,ZG( 30) ,DG(30) ,Y( 9*30) ,D(9,30 ) 

DIMENSION A ( 4 ) »B(4) »C(4) »J(7) »B2( 5) »B4 ( 5 ) *DY ( 30 ) »DYG ( 30 ) *0(30) 
M=MG 
T=TF 
H=HG 

IF(M)10,10,5 
.5 I XA = 0 
10 I = I X A+l 

90 GO TO { 100,200 *450 )♦ I 
100 IF(.M)110,110,130 
110 M=-5 

GO TO 760 
130 HALFE=0.5*H 

B2(l ) = . 0234375 
B2(2 >=“.15625 
B2(3 >=.703125 
B2(4)=. 46875 
B2 < 5 > =-. 0390625 
B4 { 1 >=“.0390625 
B4(2 > = .21875 
B(4> =“.546875 
B4(4)=l. 09375 
B4-I 5 > = .2734375 
A ( 1 > =0.5 
A(2)=0. 29289322 
A(3)=i. 7071068 
A ( 4 ) =0.1666667 
B ( 1 ) =1 • 0 
B ( 2 ) =A ( 2 ) 

B ( 3 ) =A ( 3 > 

B(4) =0.33333333 
C(1 )=A(1 ) 

C ( 2 ) =A ( 2 ) 

C ( 3 ) =A ( 3 ) 

C ( 4 ) =A ( 1 ) 

J ( 1 ) =0 
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J ( 2 ) =1 
J ( 3 ) =-l 
J ( 4 ) =0 
J ( 5 ) =-l 
J ( 6 ) =— 1 
J ( 7 ) =0 

150 I MAX = M 
M=-3 

DO 151 1=1 ♦ I MAX 

Y ( 1 * I ) =YG ( I ) 

D ( 1 » I ) =DG ( I ) 

151 Q ( I ) =0# 0 
T=T HALFE 
N = 2 

NL = 1 
L = 1 
I XA= 1 

160 LL=N-NL 

DO 161 I = 1 ♦ I MAX 
DY ( I )=D(LL» I )*H 
R=A( L )*DY ( I ) -B ( L ) *Q ( I ) 

Y ( N * I ) =Y ( LL » I ) +R 
Q ( I ) =0 ( I ) +3 • *R-C ( L ) *DY ( I ) 

161 CONTINUE 

170 I F ( M ) 750 * 1 80 *750 

180 DO 181 I = 1 » I MAX 

181 ZG { I ) = Y ( N * I ) 

GO TO 750 

200 DO 201 I = 1 * I MAX 

201 D ( N » I ) =DG ( I ) 

205 IF(L-3)210, 240*260 
210 IF(L-2)230, 220*230 
220 T=T HALFE 
230 M=-3 
NL=0 

GO TO 250 
240 M=0 
250 L=L 1 

GO TO 160 

260 I F ( M ) 270 » 270 * 1 50 
270 IF (N-5 >280*290 *280 
280 N=N 1 

GO TO 320 

290 IF(M+4 >350*300*350 
300 DO 310 1 = 1*1 MAX 
DO 310 M= 1 * 4 
Y ( M , I ) =Y { M+ 1*1 ) 

D ( M » I ) =D ( M+l » I ) 

310 CONTINUE 
320 M=-3 

T=T HALFE 
NL = 1 
L = 1 

GO TO 160 
350 M=0 
I XA = 2 
T=T H 

ASSIGN 430 TO IXC 
GO TO 360 

400 DO 410 I = 1 * I MAX 

Y ( N ♦ I ) = • 2 0* ( Y ( N-l , I )+Y( N-2 * I )+Y (N-3 * I ) +Y ( N- 4 * I ) +Y ( N-5 * I )+Vl*D(N-l * 
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1 I ) +V2*D ( N-2 » I )+V3*D ( N-3 * I )+V4*D (N-4 » I ) +V5*D ( N-5 ♦ I ) ) 

410 CONTINUE 

GO TO I XC *(430 *750) 

430 DO 440 I = 1 » I MAX 
ZG( I ) =Y ( 6 * I ) 

440 CONTINUE 

ASSIGN 750 TO IXC 
GO TO 750 

450 DO 451 1=1* IMAX 

451 D ( N ♦ I )=DG ( I ) 

K AY=K+4 

455 IF(M) 46 0*550*150 
460 IF(M+2)550*470 *490 
470 I F ( K )480»650»480 
480 IF(J(KAY) )550»550*650 
490 I F ( K >500*700*500 
500 I F ( J(KAY) >700*550*550 
550 IF(K 1560*580*570 
560 K=K 1 

GO TO 580 
570 K=K 1 
580 M=0 
T = T H 

585 IF(N-9)590*600*600 
590 N=N 1 

GO TO 620 

600 DO 610 I = 1 » I MAX 
DO 610 L= 1 » 8 
Y ( L * I ) =Y ( L+l * I ) 

D ( L * I ) =D ( L + l * I ) 

610 CONTINUE 

620 DO 621 I = 1 ♦ I MAX 

621 ZG(I)=.20*(Y(N-2»I)+Y(N 3*1) +Y( N-4* I )+Y ( N-5 » I )+Y ( N-6 » I )+Wl*D(N-2 
1*1) W2*D ( N-3 * I ) +W3*D ( N-4* I ) +W4*D( N-5 » I ) +W5*D( N-6 * I ) ) 

GO TO IXD* (400 *660 ) 

650 H=H H 
T = T H 
N = 1 1 

ASSIGN 660 TO IXD 
GO TO 620 

660 DO 670 I = 1 » I MAX 
Y ( 2 * I ) -Y ( 3 » I ) 

D ( 2 * I ) =D ( 3 * I ) 

Y(3* I )=Y( 5* I ) 

D(3* I )=D( 5*1 ) 

Y ( 4 * I ) =Y ( 7 , 1 ) 

D { 4* I ) =D ( 7 » I ) 

Y ( 5 * I ) = Y ( 9 * I ) 

D ( 5 » I ) =D ( 9 » I ) 

670 CONTINUE 
360 N=6 
K = 3 

ASSIGN 400 TO IXB 
ASSIGN 400 TO IXD 
370 V1=14.520833*H 
V2= 14.291667*H 
V3=19.50*H 
V4= 6.7083333*H 
V5=l *9791667 #H 
Wl=61 *263889 *H 
W2= 128*1388 9*H 
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W3 = 142»66667 
W4= 70 • 97 22 2 2*H 
W5=15.180555*H 
GO TO IXB* (400,620 ) 
700 H= • 5*H 
T = T H 
K=-3 

DO 720 I = 1 » I MAX 

W1=0.0 

W2 = 0 • 0 

W3 = 0 • 0 

V1=0.0 

V2 =0 . 0 

V3=0 .0 

DO 710 L = 1 » 5 
K<=6-L 
M = N KK 

VI =V1+B2 ( L ) *D( M» I ) 
W1=W1+R2 ( L ) *Y( M, I ) 
V2=V2+B4(L)*D(M,I ) 
W2 =W2 + B4 ( L ) *Y ( M * I ) 
VO=V3+B2 ( KK ) *D (M * I ) 
. W3=W3+R2 ( KK ) *Y (M » I ) 
710 CONTINUE 

G2 = Y ( N-4 » I ) 

G4=Y ( N-3 ♦ I ) 

G6=Y ( N-2 » I ) 

Y ( 8 , I ) =Y C N— 1 > I ) 

Y ( 6 » I ) =G6 
Y ( 4 » I ) =G4 
Y ( 2 » I ) =G2 
G2=D ( N — 4 » I ) 

G4=D ( N-3 * I ) 

G6 = D ( N-2 > I ) 

D ( 8 ♦ I ) =D ( N-l ♦ I ) 

D ( 6 » I ) =G6 
D ( 4 » I ) =G4 
D ( 2 ♦ I ) =G2 
Y ( 3 » I ) =W3 
D ( 3 » I ) =V3 
Y ( 5 ♦ I ) = W 1 
D ( 5 » I ) =V1 
Y ( 7 » I ) =W2 
D ( 7 » I ) =V2 
720 CONTINUE 
M=-l 
N = 9 

ASSIGN 620 TO IXB 
GO TO 370 

750 DO 751 1 = 1 * IMAX 

751 YG { I ) =Y ( N ♦ I ) 

T F = T 

HG = H 
760 MG=M 
770 RETURN 
END 





SUBROUTINE I NTERP ( T ABX , TABY , TABZ ♦ X ♦ Y *Z ♦ NX »N Y * ER ) 

TABX* ( X- TABLE ( NX ) ) 

TABY* ( Y-TABLE! NY ) ) 

TABZ ( X-ROW» Y-COLUMN ) *2-DIM. ARRAY 
X*X VALUE 
Y*Y VALUE 

Z*INTERPOLATED VALUE 

DIMENSION T ABX (NX ) » T ABY ( NY ) » TABZ ( NX >NY ) 

DO 1 J=1 *NY 
IF(TABY(J)-Y)1»3*2 

1 CONTINUE 
GO TO 100 

2 XX=1. 

CALL LIN I NT ( TABX ( 1 ) *TABZ( 1 * J-l ) *X *ZZl »OUT *NX) 

I F ( OUT ) 8 * 100 *8 

8 IF(OUT-100. )4, 100,4 

3 XX=0 • 

4 CALL LININT ( TABX ( 1 ) ,TABZ ( 1 , J ) ,X ,ZZ2 ,OUT ,NX ) 

I F ( OUT ) 9 » 100 *9 

9 IF(OUT-100. ) 5* 100 *5 

5 I F ( XX ) 6 ♦ 99* 6 

6 Z = ZZ1- ( TABY ( J 1 )-Y )*(ZZ1-ZZ2 )/ (TABY (J-l )-TABY (J) ) 
ER=0 • 

GO TO 101 

100 ER = 1. 

GO TO 101 
99 ER*0 • 

1=112 

101 RETURN 
END 



n ri 


SIBFTC LININT FULIST*SDD 

SUBROUTINE LININT (Z»Y*ARGZ,ARGY *AL* N 1 
SUBROUTINE FOR LINEAR INTERPOLATION 

DIMENSION Z ( 100 ) »Y( 100) 

6 NN = N 

7 M = N 1 
5 DO 111=1 *M 

8 LL= I +1 

9 L = I 

10 A= ( Z ( I ) -ARGZ ) * ( Z ( I+ll-ARGZ) 

12 IF(A>14,14,15 
IS AL=0.0 

11 CONTINUE 

20 I F ( ARGZ-Z ( 1 ) 125,14,26 

25 I F ( Z ( 1 ) -Z ( 2 ) )27, 14,28 

26 I F ( Z ( 1 ) -Z ( 2 ) )28,14,27 

2 7 ARGY =Y ( 2 1 — ( ( Z ( 2 1 -ARGZ ) / ( Z ( 2 ) -Z ( 1 ) 1 ) * ( Y ( 2 ) -Y ( 1 1 1 
GO TO 22 

2 8 ARGY = ( ( ARGZ-Z ( M ) ) / ( Z ( NN ) — Z ( M ) ) )#(Y(NN)-Y(M) ) + Y(M) 
AL= 1 00 • 0 
GO TO 22 

14 RAT= ( ARGZ-Z ( L) )/( Z ( LL )-Z ( L 1 1 
2? ARGY=Y (L 1 +RAT* ( Y ( LL ) -Y ( L ) ) 

24 AL=L 
22 M=N 1 
RETURN 
END 
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sibftc therma fuList,$dd 

SUBROUTINE THERMA ( PR *T »VPSI A »HLL ♦ MEN ,SL »SEN »VOLLL > VENN »CPV»CPL ♦ 
1POLLYV, AV »CHOI CE ♦ SUPER »HGG» SG»W» CP P » POLLY , A ) 

SUBROUTINE RETURNS THERMODYNAMIC PROPERTIES 
OF SATURATED AND/OR SUPERHEATED POT ASSTUM. 

<NRL PROPERTIES). 

$UPER=-1.» SP. VOL » ENTHALPY .AND ENTROPY ONLY. 

SUPER=0. » SATURATED PROPERTIES ONLY. 

SUPER=1.» SUPERHEATED PROPERTIES. 

CHOICE=l . * PRESSURE INPUT ( SATURATED ) 

CHOICE*— 1 . * TEMPERATURE INPUT ( SATURATED ) 

CHOICE=ABS GREATER THAN 6736.0, SP VOL 
ITERATIONS PRINTED. 

CHOICE=ABS GREATER THAN 1362.0* VIRIAL COEF . 

AND THEIR DERIVATIVES PRINTED. 

DIMENSION PLUSH ( 1 8 ) * TUSH (12) *VUME ( 300 ) * FF ( 1 00 ) 

DATA PLUSH/10. *2. »1.». 3333, .2,. 166667, .l25».l, 

10.083333333 *0.07 142857 *0.0625 *0.055 555556 *0.05 *0.03846 1538* 
20.33333 *0.025* 0.0 16667 *0.01 25/ 

DATA TUSH/8. 33333E-4,7.66923077E-4, 7.142857E-4, 6. 666667E-4* 
16.25E-4*5.8823529E-4,5.5555555E-4*5.263157E-4»5.0E-4» 

1 4. 761904E— 4 * 4. 5454E-4»4. 3478 E— 4 / 

DATA ( VUME ( I ) ,1=1 ,208 ) /3254.8 ,6 18. 1 , 285 . 2 ,266. *247. ,228. ,209. ,90. ♦ 

1 171.0. 152.0 *133.0, 114.0 *95.0 *76.0 *5 7.0 *38.0 *19.0* 1.0 *3545. 5* 

2691. 4,333. 6, 91.16 *38. 94 *24.7 1,5. 18, 20. 0,1 8. 0*16. 0*14. 0*12. 0*10.0* 
48. 0*6. 0*4. 0,2. 0,1 0.0,3828. 1*75 5. 13 *370.62 *112.84, 59.89 *46.22* 
428.48*17.142*9.005,2.696,2. *1.8, 1.6 *1.4 *1.2*1. *.8 *.6* 
54106.9*814.8,403.1*128.0,72.37,58.28,40.40*29.37,21.77,16.13*12. ♦ 
68.100,5.10*3.10, 1.60,1.30,1.00* .7,4383.9,872.4,433.3*140.4* 
781.47, 66.64,48.02 *36. 7, 29.03 *23. 45, 19.18*15 .78*12.99,6.86, 

83. 908 .2. 0. 1. 0*0. 5, 4659. 9 *928. 9 *462, 5, 15 1.4, 89. 02 *73. 38 ,53. 76, 

941. 92 .33. 96. 28. 23. 2 3. 88. 20. 46. 17. 68. 11. 75. 8. 99. 4. 18. 2. 0. 1.0* 

94935.2.984.8.491.0. 161.7.97.76.79.25.58.57.46.12*37.79* 
131.8,27.29*23.76,20.91,14.91 ,12.17*7.53*2.39,1.0*5210.3,1040.4, 
2519.16,171.61,102.05 ,84. 65 *62, 87, 49.78, 41.03 ,34. 76 *30. 04 *26. 36, 
323.41*17.21,14.41 ,9 . 75 ,4.49 , 19. 99 , 5485 . 1 * 1095 . 8 * 547 • 09 » 181 .28 , 
4108.09*89.78 ,66. 88 ,53. 22 *43. 94 *37. 3 7, 32. 43, 28. 58, 25. 49 *19. 026, 

516.152. 11.383. 6.419. 1.977.5759. 7. 11 51. 0. 574.9 *190. 8,1 14.0, 
694.74,70.71*56.28,46.65,39.77*34.6, 30.57,27.34,20.605,17.611, 

7 12. 68 ,7.62*4.96, 6034.3 *1206. 1*602. 6 *200.2, 119.7 *99. 58, 74.42* 

8 59.31,49.24,42.03 *36.63*32.42*29.05 ,22.02*18.9,13.78,8.54,5.88, 
96308.8 *1261. 1,630.2 ,209. 5, 125. 4, 104. 4, 78. 05 *62.27*51.74,44.22/ 
DATA ( VUME ( I ) »I =209,216) /38. 5 7, 34.17 ,30.66,23.33*20.08, 14.76,9.35* 
16.61/ 

200 PNRL=PR/14. 69594 

230 IF (PNRL-0. 007)231, 232*232 

231 PNRL=0.007 

232 IF (PNRL-21.0J233, 233*234 
234 PNRL=21.0 

233 AJ=778.26 
E=67 .0 
F=68.0 
G=69.0 
DECKA=0 
DECKB=0 
DECKC=0 
VALA=0 
VALB=0 
VALC=0 
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CHECK=CHO ICE 
J=3 
J J = 1 0 

PRENT=ABS(CHECK) 

300 TNRL=T 

330 IF (TNRL-1250. 0)331 *332*332 

331 TNRL=1250.0 

332 IF (TNRL-2700. 0)350*350*333 

333 TNRL=2700.0 
350 KAM=0 

500 TLAG= ( ALOG { TNRL ) )*. 4342945 

600 IF(CHOICE)800*800,820 
820 PSAT AT =PNRL 

GO TO 1000 

800 PVLAG= (6. 12758-81 28. 77/TNRL-O . 5 3299*TL AG ) *2. 302 5 8 5 
900 PSATAT=EXP(PVLAG ) 

1000 PSAT =P SAT AT *14*69594 
1100 PSATSF=PSATAT*2116.2 

1200 BANLUG= ( 4 544. 3 /TNRL -3. 7641 +T LAG )*2. 302585 
1300 BAN=EXP(BANLUG ) 

1400 B=BAN* (-1*0) 

1500 CLUG=( 0.6988+5 72 5.0 /TNRL) *2. 302 585 
1600 C=EXP ( CLUG ) 

1700 DANLUG=(1. 9274+5671. 4/TNRL)*2. 302585 
1800 DAN=EXP(DANLUG ) 

1900 D=DAN* (-1.0) 

2 000 DBDT = ( B/TNRL )* ( 1 .0-1 046 3. 637 /TNRL ) 

2100 TSQ=TNRL*TNRL 

2200 DCDT =13182.2 *C *(-1.0) /TSQ 

2300 DDDT = 13058. 8 *D *(-1.0 )/TSQ 

2400 T CUBE=TSQ*TNRL 

2500 ELLA=39375.0/TNRL 

2600 ELLEN=EXP (ELLA ) 

2700 H2ERA=998.95+. 127*TNRL+24836. O/ELLEN 
2800 SARAH=31126.0/TNRL 
2900 SALLY=EXP (SARAH) 

3000 SZERA=0. 18075 0 . 1 27*TLAG*2 . 302 585+0 . 76 1 7/SALL Y 
3100 RR=1 545.43 
3200 WAYT = 39.10 
3300 RRM=39. 52506 

3400 PARLUG= (408 0.0 /TNRL -1.803) *2. 302585 
3500 PAR=EXP ( PARLUG ) 

3600 VAL1=(RRM*TNRL /PSATSF-PAR ) *39. 1 
MARK=1 

8200 TFAR=TNRL-459. 69 
8300 TFAR2=TFAR*TFAR 
8400 TFAR3=TFAR*TFAR2 

8500 DUNS L= 5 2 .7274 7 . 3 539E~3*TFAR-0 . 5601 3E-6*TFA R2 +0 . 03 1 58 
1 E-9*TFAR3 
8600 VLAQ=1 .O/DUNSL 
8700 VOLLL=VLAQ 
9500 T4=TSQ*TSQ 
9700 D2B=1.094877E8*B/T4 

9800 D2C=26364.4*C/TCUBE-13182.2*DCDT/TSQ 
9900 D2D=26117.6*D/TCUBE-13058.8*DDDT/TSQ 

601 EMMA=28070. O/TNRL 
701 EMMY=EXP( EMMA) 

801 CZERA=0. 127+2. 888/EMMY 

CPLAG=0 • 22712 0 . 64848E-4*TNRL+0 . 23 1 78E-7*TS0 
CPL=CPLAQ 

IF (PRENT-1 362. 0)3700,102 *102 


- 79 - 



10? WRITE •. 6 , i 0 i , B,C,D,DBDT,DCDT,DDDT,D2B, 
1D2C,D2D»PSAT ,TNRL 
3700 VAL12=VAL1*VAL1 
3800 VAL1 3- VAL 1*VAL 12 

3 900 SLEFT = PSATSF*VAL1/ f 1545.43*TNRL » 

4000 RIGHT-1.0+B/VAL1+C/VAL12+D/VAL13 
4100 7ALLY=0.001*$LEFT 

I F (PRENT-6736. 0)4200, 103 » 103 
103 WRITE (6,104) VAL1 ,SLEFT , RIGHT »DECKA,DECKB , 

1DFCKC » VALA , VALB » VALC »CZERA »HZERA*SZ ERA »PSAT »TNRL 
4200 DECK=SLEFT-RIGHI 
4250 ADECK-ARS (DECK ) 

4300 I F i T ALL Y- ADECK (4400,4500 ,4500 
4400 I F > MARK™ 1 0 14600,4 700,4700 
4700 VAL1- (RRM*TNRL/PSATSF-PAR)*39* 1 
GO TO 4500 

4600 I F f MARK -2 5 4800,5200,5500 
4800 VALA- VAL1 
4900 DFCKA-DECK 
5000 VAL1 -VALA*1 *03 
5100 MARK -MARK+1 
GO TO 3700 
5200 VALB-VAL 1 
5300 DECK6 "-DECK 
5400 VAL1 -VALA*0« 97 
MARK "MARK-*- 1 
5450 GO TO 3700 
5500 VALC VAl 1 
5600 DECKC DECK 
5700 DECK A2 - DECK A*DECKA 
5800 DECKB2 "DECK BEDECK B 
5900 DECKC2-DF CKC*DECKC 
6000 Z 1 Z 3 -DECK A * DEC KC 
6100 till - DEC K A -D EC KB 
6200 W1 W3 -DECKA2 - D E C K C 2 
6300 W1W2 ■• DECKA2 - DECKB2 
6400 Y1Y3 V A L A V A L C 
6500 Y1Y2"VALA“VALB 
6600 ANUM--'ZlZ3*YiY2“ZlZ2*YlY3 
6700 ADEN - Wi W2*Z 1Z3-W1W3*Z1Z2 
6800 AAA-ANUM/ADEN 
6900 BBB- ( Y1Y2-AAA*W1W2 ) /Z1Z2 
7000 V A L 1 - V A L A -- A A A * D E C K A 2 - B B B * D F C K A 
7100 MARK --MARK -*-1 
7400 DECK A-DFCKB 
7500 DECKB-DFCKC 
7600 VALA-VALR 
7700 VAlB-VALC 
GO TO 3700 

4500 VENNA-VAL 1 /39. 1 
7800 VAL2 -VAL 1 *VAL1 
7900 V A L 3 - V A L 1 * V A L 2 

7950 VAL4-VAL2*VAL2 

7951 V A L 5 = V A L 2 * V A L 3 

8000 HG1 HZERA4 :RRM*TNRL/AJ)*( f B-TNRl*DBDT ) /VAL1 +( C-TNRL 
l*DCDT/2 .0 5 / VAL2+1 D-TNRL*DDDT/3. 0) /VAL3 ) 

8100 HFNNA-HG1 

8150 HGSA! -HENNA 

9200 VRT- VAL 1 / (RR*TNRL ) 

VRTAR --ARS ( VRT : 

V R i - : V R t A P 


- 80 - 


9300 VRTLUG=ALOG ( VRT ) 

9400 SGSAT=SZERA-(RRM/AJ )*( ( B+TNRL*DBDT ) /VAL1-VRTLUG 

1+ (C TNRL*DCDT) / ( 2 .0*VAL2 ) + ( D+TNRL*DDDT ) / ( 3 . 0*VAL 3 ) -7 .65738 ) 
9600 SENNA=SGSAT 
VALSAT = VE NNA 

8800 HATTY=(PSATSF/AJ ) * ( 1 871 7 • 2/TNRL-O • 53299 ) * ( VAL SAT-VLAQ ) 

89C0 HLNRLsHGSAT-HATTY 
SATTY=:HATTY/TNRL 
SLNRL=SENNA-SATTY 
IF(SUPER) 2001.201 .201 

201 BRAKET=1 • 0+ ( B TNRL*DBDT ) /VAL 1+ ( C+TNRL*DCDT ) /VAL2 
1+ ( D TNRL*DDDT) /VAL3 
301 CGNUM=BRAKET *BRAKET 

401 CGDEN=1.0+2.0*B/VAL1+3.0*C/VAL2+4.0*D/VAL3 
501 PAR2 =TNRL*D2B 2.0*DBDT+0.5*( TNRL*D2C+2 .0*DCDT ) /VALl 
1+ ( TNRL*D2D+2 .0*DDDT ) / ( 3 .0*VAL2 ) 

901 CPG=CZERA+ ( RRM/AJ )* ( CGNUM/CGDEN-1 .0-( TNRL/VAL 1 ) 

1*PAR2 ) 

1050 BRUCK=1 .0/VAL2+2 . 0*B /VAL3+3 • 0*C/V AL 4+4 . 0*D/ VAL 5 
1101 DRHDT=RR*TNRL*(-1»0)*8RUCK 
1201 DVDP*1.0/(DRHDT*WAYT ) 

1301 DVDT= (PSATSF / ( RR*TSQ ) +DBDT/VAL2+DCDT/VAL3+DDDT/VAL4 ) 

1/ (BRUCK*WAYT ) 

1401 PALE=1 .0/ ( (PSATSF/VENNA)*(DVDP+(TNRL/(AJ*CPG) )*DVDT*DVDT) ) 
1501 PALM =P ALE * (-1.0) 

1601 ALF=PALM*32.174*RRM*TNRL 
1701 ABALF=ABS (ALF ) 

1801 ALE=SQRT ( ABALF ) 

1901 IF(KAM)2101 .2001 .2101 

2001 VENN=VENNA 
2201 HFN=HFNNA 
2301 SEN=SENNA 
SL=SLNRL 
9000 HLL=HLNRL 
9100 VPSI A=PSAT 

2351 IF(SUPER) 2352.2401 ,2401 

2352 CPV=0.0 

2353 P0LLYV=0 • 0 

2354 AV=0 • 0 

GO TO 2801 
2401 CPV=CPG 
2501 POLLYV=PALM 
2601 AV=ALE 

2701 IF(SUPER) 2801.2801 ,3401 
2801 HGG=0 
2901 SG=0 
3001 VV=0 
3101 CPP=0 
3201 POLLY=0 
3301 A=0 
3351 GO TO 4601 
3401 KAM=1 
3501 MARK=1 
3601 PSAT AT =PNRL 
3701 PSAT =PSAT AT *14 .69594 
3801 PSATSF=PSATAT*21 16.2 
MM = 1 8 
NN = 1 2 

RPSAT=1.0/PSAT 
RTNRL=1 .O/TNRL 

3903 CALL DUB I N ( RPSAT ♦ RTNRL , VELD » PLUSH ♦ T USH , VUME »MM ,NN ♦ JAN ♦ JAQ ) 


- 81 - 



MARK=1 


4001 

VAL1 = VELD *3 9 • 1 




GO TO 

3700 



2101 

HGG=HENNA 



4101 

SG=SENNA 

CPP=CPG 



4201 

VV=VALl/39. 1 



4301 

CPP=CPG 



4401 

POLLY= 

PALM 



4501 

A=ALE 




4601 

HARRY= 

6,0 




RETURN 




101 

FORMAT 

GENERATOR 




B 

C 

D 



-1RE3 

-E5 


*E5 


SPACE 

D2B 


D2C 



-E5 

-E5 



SPACE 
END OF 

FORMAT 



104 

FORMAT 

GENERATOR 




SPACE 

VAL1 

LEFT 




-E 5 


-E5 


SPACE 

DEC3 

VOL1 




-E5 


-E5 


SPACE 

HO 

SO 




-E5 


-E5 


SPACE 
END OF 
END 

FORMAT 




DBDT 

-E5 

DCDT 

-E5 

DDDT 

-E5 

-E 5 

D2D 

PRESS 

-E5 

TEMP 

-E5 

RIGHT 

-E5 

DEC1 

-E5 

DEC2 

-E 5 

V0L2 

-E5 

V0L3 

-E5 

CO 

-E5 

PRESS 

-E5 

TEMP 

-E5 
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ana 


fIBFTC DUB I N FUl_IST,SDD 

SUBROUTINE DUB I N ( P * T * V , PLUSH ♦TUSH » VUME * MM »NN » JA * JB ) 

SUBROUTINE FOR DOUBLE LINEAR INTERPOLATION 

DIMENSION PLUSH ( 40 ) »TUSH(40) »VUME(1600) 

100 PP=P 
200 TT=T 
600 M=MM 
700 N=NN 

800 DO 1100 I =1 ♦ M 
850 L = I 

900 C=(PLUSH( I )-P)*(PLUSH( I 1 ) -P ) 

1000 IF(C)2700»2700»1200 
1200 JAM=3 
1100 CONTINUE 

1400 TEST=(PLUSH(1)-P)*(PLUSH(2) —PLUSH ( 1 ) ) 

1500 IF(TEST) 1700*2100*2500 
1700 JAB=100 
1750 MN=M-1 

1800 RAP= { P-PLUSH ( MN ) > / ( PLUSH ( M ) -PLUSH ( MN ) ) 

LL=MN 

1950 GO TO 2800 
2100 RAP=0.0 
2200 LL=1 

JAB=200 
GO TO 2800 

2500 RAP=( P-PLUSH (1 ) )/( PLUSH ( 2 ) -PLUSH ( 1) ) 

2550 JAB^O 
2600 LL=1 
2250 GO TO 2800 
2700 LL=L 

LLL^LL+l 
2950 JAB=LL 

2 900 RAP=( P-PLUSH (LL) )/( PLUSH ( LLL ) -PLUSH ( LL ) ) 

2800 LLL=LL+1 
3000 DO 3500 I=1»N 
3100 K=I 

3200 C = ( TUSH ( I ) -T )* ( TUSH { 1 + 1 )-T) 

3300 IF(C)4500»4500»3400 
3400 JAM= 1 0 
3500 CONTINUE 

3600 TEST = ( TUSH ( 1 ) T ) * ( TUSH ( 2 ) -TUSH 1 1 ) ) 

3700 IF(TEST)3800, 4100, 4300 
3800 N J=N-1 

3900 RAT=(T-TUSH(NJ) )/( TUSH ( N ) -TUSH { NJ ) ) 

KK=N J 

3850 JBB= 100 
4050 GO TO 4600 
4100 RAT=0.0 
4200 KK=1 

JBB=200 
GO TO 4600 

4300 RAT= (T-TUSH ( 1 ) ) / ( TUSH ( 2 ) -TUSH ( 1 ) ) 

4350 JBB=0 
4400 KK=1 
4250 GO TO 4600 
4500 KK=K 

KKK=KK+1 
4675 JBB=KK 

4650 RAT= ( T-TUSH ( KK ) )/( TUSH ( KKK ) -TUSH ( KK ) ) 
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4600 KKK=K+1 
4700 KV1=M*KK+LL-M 
4800 KV2=KV1+1 

4900 V1=VUME (KV1 ) +R AP* ( VUME ( KV2 ) -VUME ( KV 1 ) ) 
kvi=m*kkk+ll-m 
5100 KV2=KV1+1 

5200 V 2= VUME ( K VI ) +RAP* { VUME ( KV2 ) -VUME ( KV 1 ) ) 
5300 VV=V1+RAT*(V2 Vl) 

5400 V=VV 
5500 JA= JAB 
5600 JB= JBB 
RETURN 
END 



5IBFTC CURV FULIST *SDD 
FUNCTION CURV(P) 

17 PYP=P/14. 69594 

18 PUN=(ALOG(PYP) )/2. 3025851 

19 PEEP=(PUN+0« 10705539 ) /2. 1805446 

20 PUP-PEEP*PEEP 

21 PUPPY=PUP*PEEP 

22 PURP=PUP*PUP 

23 PURPLE=PUP*PUPPY 

24 PURS E=PUPPY* PUPPY 

25 PERK=PURP*PUPPY 

26 PERKY=PURP*PURP 

27 PLUP*1. 31 90232*PEEP-.13299096+. 6923 7968#PUP+. 37008034 
1*PUPPY+. 1 98 76127* PURP+* 10754 778 *PURPLE+. 066 299672 
2*PUR$E+ .041 265 195*PERK+ .0134 72973 *PERKY 

28 CURV=1900. 0+700. 0*PLUP 
RETURN 

END 



APPENDIX B 

POTASSIUM BOILER THERMAL DESIGN CODE 
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TABLE IV 


NOMENCLATURE FOR BOILER DESIGN CODE INPUT* 


Symbol 

Quantity 

Units 

TITLE 

Identification block 

- 

NT 

Number of tubes 

- 

I TP 

Number of iterations for primary exit temperature 

- 

ITS 

Number of iterations for tube length in boiling 
region 

— 

ITSUP 

Number of iterations for tube length in 
superheat region 


JHP 

Enter JHP = 0 if primary heat transfer coefficient 
is calculated by program 

Enter JHP = 1 if primary heat transfer coefficient 
is input 


POWER 

Total power 

KW 

TPINS 

Primary fluid inlet temperature 

°F 

TSIN 

Secondary fluid inlet temperature 

°F 

BTSUP 

Degree of superheat 

°F 

PSOUT 

Secondary outlet pressure 

psi 

DT 

Outside tube diameter 

inches 

XW 

Tube wall thickness 

inches 

PODT 

Pitch to diameter ratio of tube bundle 

dimensionless 

PODI 

Pitch to diameter ratio of insert 

dimensionless 

DCB 

Center body diameter 

inches 

HB 

Average boiling heat transfer coefficient 

Btu/ft 2 -hr-°F 

HSUP 

Average superheated vapor heat transfer 
coefficient 

Btu/ft 2 -hr-°F 

HP 

Average primary fluid heat transfer coefficient 

Btu/ft 2 -hr-°F 

HPMARG 

Margin factor for tube bundles with tight pitch 

dimensionless 


♦Input variables are given in order of literal appearance on input cards 
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TABLE IV (Cont'd) 


Symbol 

Quantity 

Units 

A 

GS 

Boiling fluid mass velocity (secondary) 

2 

lb /ft -sec 
m 

WP 

Primary fluid mass flow rate 

lb /sec 
m 

XINST 

Insert vane thickness 

inches 

RADPL 

Length of radial pressure drop region 

inches 

TOLP 

Tolerance on primary exit temperature iterations 

°F 

TOLS 

Tolerance on iteration for tube length in 
boiling region 

Percent 

100 

TOLSUP 

Tolerance on iteration for tube length in 
superheat region 

Percent 

100 

PF 

Primary fluid indicator 

0 = Lithium 

1 = Sodium 


CSLOPE 

Sign of iteration line slope. 
Use 1 

- 

CS 

Sign of iteration line slope. 
Use 1 

— 

XN 

Exponent on friction factor 
model; smooth tube = 0.25 

dimensionless 
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TABLE V 


NOMENCLATURE FOR BOILER DESIGN CODE PRINTOUT 


Symbol 

Quantity 

Units 

ACCL 

Radial acceleration 

ft/sec 2 

ACCL/G 

Ratio of radial acceleration to the local 
gravitational acceleration 

g’ s 

DELTAP 

Primary friction pressure drop or Secondary 
total pressure drop in boiling region 

psia 

DPFRIC 

Secondary friction pressure drop in boiling 
region 

psia 

DPMOM 

Secondary momentum pressure drop in boiling 
region 

psia 

DPRIN 

Primary radial pressure drop at inlet 

psi a 

DPROUT 

Primary radial pressure drop at exit 

psia 

DPSUP 

Secondary total pressure drop in superheat 
region 

psia 

DPSFRIC 

Secondary friction pressure drop in superheat 
region 

psia 

DPS MOM 

Secondary momentum pressure drop in superheat 
region 

psia 

GMAX 

Secondary maximum helical mass 

lb /ft 2 -sec 

GP 

Primary mass velocity 

m 2 

lb /ft -sec 
m 9 

HB 

Boiling heat transfer coefficient 

Btu/ft -hr-°F 

HP 

Primary heat transfer coefficient 

Btu/ft 2 -hr-°F 

HSL 

Secondary heat transfer coefficient in liquid 
region 

Btu/ft 2 -hr-°F 

HSUPS 

Secondary heat transfer coefficient in superheat 
region 

Btu/ft 2 -hr-°F 

HWALL 

Effective wall heat transfer coefficient 

Btu/ft 2 -hr-°F 

LTOT 

Total tube length 

ft 

LB 

Tube length in boiling region 

ft 

LSC 

Tube length in subcooled region 

ft 

LSUP 

Tube length in superheat region 

ft 

NT 

Number of tubes 

dimensionless 
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TABLE V (Cont'd) 


Symbol 

PEL 

PEP 

POWER 

P/D 

QB 

QSC 

QSUP 

QTOT 

Q/A BOILING 

REP 

REL 

TPIN 

TPOUT 

TSIN 

TSAT 

TSOUT 

UIB 

UIL 

VMAX 

VP 

WP 

ws 

xw 


Quantity 


Secondary fluid Peclet number 
Primary fluid Peclet number 
Total power 

Pitch to diameter ratio 

Heat transfer rate in boiling region 

Heat transfer rate in subcooled region 

Heat Transfer rate in superheat region 

Total heat transfer rate 

Heat flux in boiling region 

Primary fluid Reynolds number 

Secondary fluid Reynolds number 

Primary inlet temperature 

Primary exit temperature 

Secondary inlet temperature 

Secondary temperature at inlet to boiling region 

Secondary exit superheated vapor temperature 

Overall heat transfer coefficient in boiling 
regi on 

Overall heat transfer coefficient liquid region 

Maximum secondary helical velocity 

Primary fluid axial velocity 

Primary fluid flow rate 

Secondary fluid flow rate 

Wall thickness 


Units 

dimensionless 

dimensionless 

KW 

dimensionless 

Btu/sec 

Btu/sec 

Btu/sec 

Btu/sec 

2 

Btu/ft -sec 

dimensionless 

dimensionless 

°F 

°F 

°F 

°F 

°F 

2 

Btu/ft -hr-°F 

Btu/ft 2 -hr-°F 

fps 

fps 

lb /sec 
m 

lb /sec 
m 

inches 
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BOILER CODE LISTING 
4IBFTC PBDES FULIST »SDD 

*pbdes 

C 

DIMENSION TC L I ( 20 ) * RHOLI (20) »CPLI (20) *VISLI(20) ♦ 

1 TCNA ( 20 ) » RHONA ( 20 ) » CPNA ( 20 ) ♦ V I SNA ( 20 ) »V I SK ( 20 ) » V I SKG ( 20 ) ♦ 
2TCK (20 ) »TEMP(20 ) »XKW( 20) »TW(20) 

C 

DIMENSION TITLE(8) »BITS( 1) >DUMMY ( 50 ) * REG( 5 ) 

DIMENSION FOFO( 4*4 ) * POD( 4) * RES ( 4 ) »EPSMAX ( 5 ) 

DIMENSION TCKG ( 20) 

C 

NAMELIST/INPUT/TITLE* NEX*NT »ITP» ITS»ITSUP» POWER » 

1TP INS»TSI N»DTSUP »DT »XW »XKW »PODT * POD I »HB »HBTB » 

2HSUP»T0LP »TOLS * TOLSUP * HP , HPM ARG* PF»TW*C SLOPE* 

3CS » XN » FOFO »GS » DC8 » WP »PSOUT >XINST * RADPL » JHP 

1 READ ( 5 » I NPUT ) 

C 

C PROGRAM IS APPLICABLE TO A COUNTERFLOW 

C BOILER WITH SUBCOOLED FEED AT INLET 

C AND 100 PERCENT QUALITY VAPOR AT EXIT. 

C 

2 QTOT=.9478*POWER 
D I T=DT-2.*XW 
DLM=2.*XW/AL0G( DT/DIT I 

C 

C COMPUTE GS OR NT FROM TOTAL POWER 

C AND SECONDARY FLOW RATE 

400 TSINR=TSI N+459 .69 

CALL THERMA (PU * TSI NR » PSI NR »H L I N , HG7 * SL7 *SV7 »VL7 » 
1VG7»CPG7»CPL7*EX*AV»-1.»0. *HS7 *SS7 » VS7 *CPSH7 * EXS7 » AVS7 ) 

C COMPUTE EQUIVALENT SAT. TEMP. AT PSOUT. 

C 

300 TKSATO=CURV (PSOUT ) 

TKSAT 0=TK.SAT0-459.69 
2971 CONTINUE 

TK.SUPO = TKSATO+DTSUP 
T9 9 = TK.SUPG+459 .69 

CALL THERMA (PSOUT ,T99»PV7 ,H L 7 * HG 7 , SL 7 * SV 7 » VL7 » VG7 * 

1CPG7 »CPL7 »EX »AV *0. ♦ 1. *HSUOUT »SS7 * VSUPO* CPSH7 , EXS7 * AVS7 ) 
WS=QTOT/ (HSUOUT-HLI N ) 

I F ( GS )401 *401 *402 

401 AS =3. 1416*DI T**2 /4.*FL0AT ( NT ) / 144 . 

GS =WS/ AS 

GO TO 403 

402 ATUBE = 3.1416*DIT **2 /4./144. 

XNT=WS/GS/ATUBE 
NT=IFIX(XNT ) 

403 CONTINUE 
C 

C DETERMINE PRIMARY EXIT TEMPERATURE AND 
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n n 


PROPERTIES AT AVERAGE PRIMARY TEMPERATURE. 

RADI =D I T/ 2 • -DCB / 2 • 

ABLOCK=RADl*XINST*FLOAT (NT )/144. 

AFAX=3. 1416* (DI T**2-DCB**2 ) / 4 . / 1 44 . *F LOAT ( NT ) -ABLOCK 
GAX=GS*AS/AFAX 

BETA = SQRT ( 1 .+ ( 3 . 1416/ POD I ) **2 ) 

GMAX=BETA*GAX 

DCOD s DCB/DI T 

IF (DC 8)6111*6111*6112 

6112 CONTINUE 

DE NOD= ( 1 .-DCOD**2 ) / ( 1 .+DCOD + ( 1 .-DCOD ) /3 . 1416) 

DEN=D I T *DENOD 
GO TO 6113 
6111 DEN=DI T 

6113 CONTINUE 
IT 1 = 1 

TP = TP I NS 

6 IF (PF ) 3 .3 »4 

3 CALL LI N I NT ( TEMP »CPL I » TP* CPAV *OU T ♦ 20 ) 

GO TO 5 

4 CALL LININT (TEMP »CPNA »TP»CPAV»OUT*20) 

5 TP 1 = TP I NS-QTOT /WP/CPAV 
TZ= (TP1+TPINS ) /2. 

IF (ABS ( TZ-TP ) -TOLP ) 7* 7 *8 

8 IF ( IT1-ITP) 10*11 .11 
11 CALL EXIT 
10 IT1=IT1+1 
TP = TZ 
GO TO 6 

7 TPA=TZ 

TPOUT=TPI NS-tTOT /WP/CPAV 
IF(PF) 15*15*16 

15 CALL LININT ( TEMP.TCLI > TPA* TC P » OUT ♦ 20 ) 

CALL LI NINT ( TEMP *RHOL I » T PA » RHOP ♦ OUT * 20 ) 

CALL LININT (TEMP*VISLI >T PA *VISP» OUT *20) 

CALL LI NINT ( TEMP * RHOL I*TPINS»RHOIN»OUT*20) 

CALL LI NINT (TEMP .RHOL I »TPOUT »RHOUT *0UT»20) 

CALL LI NI NT ( TEMP. V I SLI .TP I NS, VI SIN, OUT *20) 

CALL LININT (TEMP.VISLI .TPOUT .VISOUT .OUT ,20) 

CPP=CPA V 

16 CALL L I NI NT ( TEMP * TCNA * TPA » TC P * OU T » 20 ) 

CALL LININT (TEMP* RHONA » T PA * RHOP * OUT » 20 ) 

CALL LININT ( TEMP * V I SNA » T PA , V I SP , OUT, 20) 

CALL L I N I NT ( TEMP » RHONA * TP IN S » RHO I N » OUT .20) 

CALL LININT ( TEMP .RHONA, TPOUT .RHOUT .OUT ,20) 

CALL LININT ( TEMP » V ISNA » TP I NS * V I S I N * OUT .20 ) 

CALL LI NINT (TEMP.V ISNA .TPOUT .VISOUT .OUT *20) 

CPP=CPAV 
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COMPUTE SHELL SIZE AND PRIMARY PARAMETERS. 

IF (PODT )2322»2322*2321 
322 SHI D= BITS 
SHFA= BITS 
GP=BITS 
VP = B I T S 
RE P=B I TS 
GO TO 2323 

321 SHID=DT*PODT*SQRT(2.*SQRT<3. ) /3 . 1416*FLOAT ( NT ) ) 
SHFA= FLOAT (NT ) *3 • 1416*DT**2 /4»* ( 2.*SGRT (3. )/ 
13.1416*PODT**2-l. ) / 144 • 

NT 1 = 1 
1 = 1 

111 N T 1 =NT 1 +6* ( 1-1 ) 

IF (NT1-NT ) 41 1 2 *4113*4114 

112 1 = 1+1 

GO TO 4111 
114 1=1-1 
113 CONTINUE 

PIT CH = PODT *DT 

RA T 1 0=P I TCH/SQR T (SHID**2-P ITCH** 2 ) 

ANG=ATAN ( RAT 10 ) 

N0UTR=6*( 1-1 ) 

XNROWS=FLOAT (NOUTR) 

XNT1=FL0AT ( NT ) 

WETPER=3.1416*(SHID+XNT1*DT)+XNR0WS* (2.* 

1 P I T CH/ SQR T ( 3 » )-SHID*ANG) 

DEQSH=4.*SHFA*12./WETPER 

GP =WP /SHF A 

VP=GP/RHOP 

REP=DEQSH*3600.*GP/VISP 
PRP=CPP*V ISP/TCP 
PECP= RE P*PRP 

CALCULATE HP USING DWYER TU RELATION. 

HMARG=FUDGE FACTOR FOR CONSERVATISM. 

2323 IF ( JHP) 17»17»18 
17 CONTINUE 

REPL=ALOG10(REP) 

CALL LI NI NT ( REG » EPS MAX » RE PL * EPS > OUT ,5 > 
EPSP=EXP(2.303*EPS) 

PSIBAR=1.-1 .82/PRP/EPSP**l .4 
IF (PS I BAR) 555 5*5553 *55 56 

5555 PNU=.93+10.81*PODT-2.01*PODT **2 
GO TO 5557 

5556 CONTINUE 

PNU=.93+10.81*PODT-2 . 0 1*P0DT ** 2+ . 025 2 *PODT ** . 2 7 3 
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1*( PSI6AR*PECP ) **.8 
5557 CONTINUE 

‘ HP=PNU*TCP/DEQSH*HPMARG 
18 CONTINUE 
C 

IF (PF)280»280»281 

280 CALL LININT ( TEMP»CPLI .TPINS.CPRI MS .OUT .20 ) 

GO TO 282 

281 CALL LININT (TEMP.CPNA.TPINS.CPRIMS.OUT *20) 

282 CONTINUE 
IF(DTSUP)37l2»3712»2972 

3712 TP I N c TP I NS 

T S0UT = T KS ATO 
QSUP = 0 • 

HL SUP = 0 • 

DELPSU = 0 • 

DPFSUP=0. 

DP SMOM= 0 • 

GO TO 377 

ENTER ITERATION LOOP FOR DETERMINATION OF 
LENGTH OF TUBES IN SUPERHEAT REGION. 

2972 CONTINUE 
TKMAX=TKSUPO 
TKMIN=TKSATO 
BALP=0. 

301 CALL D I CHOA( VAR.TKMAX.TKMIN.BALP.TOLSUP.es .ICON. 

1TSOUT.ITSUP) 

2973 IF(IC0N)302, 302.377 

COMPUTE THE SUPERHEATED VAPOR PROPERTIES AT 
TSUPAV AND PSUPAV. 

302 TSUPAV=(T SOUT + T K.SUPO ) / 2 • 

T7 7 = T SUPA V+459 .69 
T88=TSOUT+459.69 

CALL THERMA(PU»T88, PVOUT ,HL7 ,HG7 . SL 7 » SV 7 , VL7 . VGOUT ♦ 

1CPG7 »CPL7 .EX.AV » -1 • »0. .HS7.SS7 »VS7 .CPSH7.EXS7 »AVS7 ) 
PSUPAV=(PSOUT+PVOUT) /2. 

CALL THERMA ( PSUPAV, T77 ,PV7 »HL7 ,HG 7 » SL 7 , SV7 ♦ VL 7 ♦ VG7 ♦ 

1CPG7.CPL7 .EX.AV. 0. , 1 . » HSUP AV , SSU PAV , VSU PAV » CPSUAV » EXSAV ♦ AVS AV ) 
CALL LININT ( TEMP , V I SKG , TSUPAV , V I SUPA .OUT, 20 ) 

CALL LININT ( TEMP , TC KG . TSUPAV , T CKSUP , OUT , 20 ) 

SET TEST3=0 • TWLA3=SUM OF DELTAT S. 

QSUP*WS*CPSUAV*( TKSUPO--TSOUT ) 

TPIN = TPI NS-QSUP / wS/CPRIMS 

TWLA3= ( TPIN+TPINS+TSOUT+TKSUPO)/A. 

TEST3=0 • 

304 CALL LININT (TW.XKW.TWLA3.XKW3.OUT »20) 
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RESW3=XW*DIT/ 12. /XKW3/DLM 
HSUPS=HSUP 
IF(HSUP>309, 310.309 
310 PRG=CPSUAV*VISUPA/TCKSUP 

REGSU=DEN*GMAX*300./VISUPA 
IF (DCB)6211» 6211. 6212 

6212 CONTINUE 

ZNSUP=.39*REGSU**.56*PRG** ( 1./3. ) 

GO TO 6213 

6211 ZNSUP=.023*REGSU**.8*PRG**( 1./3. ) 

6213 CONTINUE 

HSUPS=ZNSUP*TCKSUP*12 ./DEN 
309 USUP=l./( 1./HSUPS+RESW3+DI T/HP/DT ) 

T51 = TPINS-TK.SUP0 
T5 2 = TP I N-T SOUT 

TLMSUP= ( T31-T52 ) /AL0G( T51/T52 ) 
ASUP=3600.*QSUP/USUP/ TLMSUP 
IF ( TEST3 ) 305.305.320 
305 DTSOFM=QSUP*3600./ASUP/HP*DT/DIT 
DTSIFM=QSUP*3600./HSUPS/ASUP 
TW0S= ( TPIN+TPINS ) /2.-DTS0FM 
TWIS= ( T SOUT + T KSUP0J/2.+DTSI FM 
TWLA3=(TW0S+TWIS)/2» 

TEST3=1 * 

GO TO 30A 

320 CONTINUE 

HLS UP =ASUP*12./ 3.1416 /DIT /FLOAT (NT) 

COMPUTE SUPERHEAT LENGTH BASED ON THE 
PRESSURE DROP RELATIONSHIP. 

REGSU=DEN*GMAX*300./VISUPA 

321 CALL FRICtPODI » DCB » D I T ♦ GS ♦ V I SUP A * FR ) 
DELPSU=PVOUT-PSOUT 

DPSMOM=GMAX**2/32.2/ 144.* ( VSUPO-VGOUT ) 

RHOSAV=l. /VSUPAV 

PPS=DELPSU-DPSMOM 

GROUPS=GAX**2*BETA**3/64.4/RHOSAV*12 . /DEN 
PL SUP= 1 44. *PPS/ GROUPS /FR 
DPFSUP=DELPSU-DPSMOM 


COMPUTE VAR AND TO 301 

VAR= (HLSUP-PLSUP ) /HLSUP 
GO TO 301 

377 CONTINUE 

SET UP ITERATION LOOP FOR DETERMINATION 
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OF TUBE LENGTH. 


* JK =0 

TSLAV=(TSIN+TS0UT)/2. 

GO TO 41 

57 T SMAX = TPSA T -2 • 

T SM IN = T SOUT 
BAL=0. 

30 CALL D I CHOA ( V » T SMAX »TSMIN»BAL»TOLS » CSLOPE * KCON » TSAT * I T S ) 
IF (KCON)31»31»100 

COMPUTE RELEVANT PROPERTIES AT T SA T , T SA V , TSOUT . 

31 TSAV= ( TSAT+TSOUT ) /2 . 

TKl=TSAT+459.69 

TK2=TSAV+459.69 

TK3 = TS0UT +45 9.69 

CALL THERMA (PU»TK1*PV1 *HL1 * HG 1 * SL 1 * S V 1 » VL 1 * VG 1 * CPG 1 >CPL1» 
1EX , ALPHA, -1. ,0. ,HS»SS*VSH*CPSH»EXS»AVS) 

RHOLl= 1 • / VL1 

CALL THERMA ( PU , T K2 , PV2 » HL2 > HG2 » S L2 * S V2 » VL2 » VG2 * CPG2 » C PL2 , 
1EX , ALPHA, -1. ,0. »HS»SS»VSH»CPSH»EXS»AVS) 

RHOL2 = 1 • / VL2 
RHOG2 = 1 . / VG2 

CALL LININT (TEMP, VISK,TSAV»VISAV, OUT ,20 ) 

CALL LININT (TEMP, TCK*TSAV,TCAV, OUT, 20 ) 

CALL THERMA (PU , TK3 ,PV3,HL3 ,HG3 ,SL3 ,SV3,VL3 ,VG3 ,CPG3 ,CPL3, 
1 EX , ALPHA , - 1 • ,0. » HS » SS » VSH » C PSH » EXS » A VS ) 

RHOL3= 1 • / VL3 

40 RHOG3= 1 • / VG3 

COMPUTATION OF SUBCOOLED REGION 

41 TK4=TSLAV+459.69 
TK4=TSLAV+459.69 

CALL THERMA ( PU , T K.4 , PV4 , HL4 , HG4 , S L4 , SV4 , VL4 , VG4 ,CPG4 ,CPL4 , 
1EX , ALPHA, -1. ,0. »HS»SS»VSH»CPSH»EXS»AVS) 

©SC=WS*CPL4*( TSAT-TSI N ) 

TP SAT = TPOUT + dSC / WP/C PP 
IF (JK) 55,55,56 

55 JK = 1 

GO TO 57 

56 CONTINUE 

SET TEST1=0. TWLA1=SUM OF DELTAT S 
TWLA1= ( TPOUT+TPSAT+TSI N+TSAT ) /4. 

TEST1=0. 

42 CALL LININT (TW»XKW»TWLA1 »XKW1»0UT ,20) 

CALL LININT ( TEMP , V I SK, TS LAV , V I SKL , OUT ,20) 
REAL=DIT*300.*GS/VISKL 
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CALL LININTf TEMP* TCK*TSLAV,TCKL, OUT ,20) 

PREAL=CPL4*VISKL/TCKL 

PEEAL=REAL*PREAL 

IF ( PEE AL- 1 00 • ) 45 » 45 » 46 

45 ZNUL=3. 

GO TO 47 

46 ZNUL=7 .+.025* ( PEEAL ) **.8 

47 HSL=ZNUL*I2.*TCKL/DIT 
RES W=X W*D IT/ 12. /XKW1/DLM 

U I L= 1 * / ( l./HSL+RESW+DIT/HP/DT ) 

DTLMSC=( TPOUT-TS IN-TPS AT+TSAT ) / A LOG ( ( TPOUT- 
1TSIN) / (TPSAT-TSAT) ) 

ASC=QSC/UIL/DTLMSC*3600. 

IF (TEST1)43*43*44 

43 DTOF = QSC*3600. / HP / ASC*DT /D I T 
DT IF=QSC*3600. /HSL/ASC 

T WO= ( TPOUT + T PSA T J/2.-DTOF 
TWI = ( TSIN + TSAT ) /2. + DT IF 
T W LAI = ( TwO+TWl ) /2. 

TEST1= 1 . 

GO TO 42 

COMPUTE BOILING LENGTH BASED ON HEAT TRANSFER. 
SET TEST2=0. 

44 TWBA= ( TPSAT+TSAT+TPI N+TSOUT )/4. 

TEST2=0 • 

62 CALL LI NINT ( T W » XKLW » T WBA * XKW2 *OUT *20) 
RESW8=Xw*DIT/12. /XKW2/DLM 

U I 6= 1 . / ( 1 ./Nb + RESWB + DI T/HP/DT ) 
DTLMB=(TPSAT-TSAT-TPIN+TSOUT ) /ALOG ( ( TPSAT- 
1TSAT) / ( TP I N-T SOU I ) ) 

§6 = QTO T-QSUP-QSC 
AB=QB*3600. /UI B/DTLMB 
IF(TEST2)63, 63*64 

63 DTOF=QB*3600./HP/AB*DT/DI T 
DTIF=QB*3600./HB/AB 
TWOB=(TPSAT+TPIN )/2 .-DTOF 
T W I B= ( TSAT+TSOUT ) /2.+DTIF 
TrtBA= ( TwOB+TWl 3 ) /2 . 

TEST2 = 1 . 

GO TO 62 

64 HLEN=AB*12./3. 1416/DIT /FLOAT (NT) 

COMPUTE PRESSURE DROPS AND LENGTH BASED 
ON THE PRESSURE DROP 

PRBAV=CPL2*VISAV/TCAV 
RE = DI T*GS*3600. / 12. /VISA V 
PECAV=RE*PRBAV 
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CALL FR I C ( POD I , DCB » DT >GS » V I SA V , F R ) 

CALCULATE TWO PHASE MULTIPLIER FROM 
HOMOGENEOUS MODEL USING CONSTANT Q/A 

CALL LININT ( TE MP » V I SKG » TSAV » V I SG * OUT » 20 ) 
AA=RHOL2/RHOG2-l . 

BB=VISAV/VISG-1. 

R3 =AA/BB**2 

R4=(l.+BB)**<l.-XN)/BB/< l.-XN) 

R5=1./BB/ (l.-XN) 

R6 = l./ (XN-1 . )/ ( l.+BB )**( .XN-1 . ) 

R7=l./ (XN-2. ) / ( l.+BB)**(XN-2. ) 

R8=l./ (XN-2 . ) 

R9 = 1 • / (XN-1 . ) 

PHI=R4-R5+R3*( R6-R7+R8-R9 ) 

DELPT=PV1-PV3 

DPMOM=GMAX**2/3 2.2/144.*( VG3-VL1 ) 

PP = DELP T-DPMOM 

GROUP=GAX**2*BETA**3/64.4/RHOL2* 12. / DEN*PH I 
PLEN=144.*PP/GR0UP/FR 

COMPUTE V AND RETURN TO 30 

V= (HLEN-PLEN) /HLEN 
GO TO 30 


100 IF (PODT 12391 ,2391 ,2392 

2391 DE LPP = B I TS 
DPRAD1 =B I TS 
DPRAD2 =B ITS 
GO TO 2393 

2392 FP=.316/(REP)**.25 

DE LPP= FP*HLEN*GP**2 / 64. 4/DEO SH/R HOP/ 144. 

2393 IF (RADPD2395, 2393, 2394 

2394 RDDP=RADPL 
DB = DT 

WB = WP 
PODB = PODT 
NB=NT 

CALL RDP ( DB,RDDP,WB»RHOIN,VI SIN,PODB»NB , DPR ADI ) 
CALL RDP ( DB»RDDP*WB»RHOUT»VI SOUT , PODB »NB * DPRAD2 ) 

2395 CONTINUE 
DPFRI C=PP 

THETA=ATAN(PODI/3.1416 ) 

VVTMAX=GMAX/RHOG3*COS (THETA ) 

RAD=DI T/24. 

ACCMAX= (VVTMAX) **2/RAD 
AOG= ACCMAX/ 32 . 2 
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HLSC=ASC/3.1416*12./DI T /FLOAT ( NT ) 
HLETOT=HLSC+HLEN+HLSUP 

DELPP=FP*HLET0T*GP**2/64.4/DEQSH/RH0P/ 144. 

HWALL=1./RESWB 

HFLUXB=QB*3600. /AB 


WRITE ( 6 » 2000 ) ( T I T L E ( K ) , K = 1 * 8 ) , NT » PODT , D I T » 

IAS »DT » POD I »XW,SHID*HLETOT »ShFA*HLEN*HLSC» 

2POWER »QTOT ,QB * QSC ,HFLUXB *TPI NS*TSI N * TPOUT * 

3TSAT »WP » TSOUT »GP »WS 

WRITE (6*7000) VP , GMAX ♦ REP » VVTMAX , PECP > ACCMAX »HP ♦ AOG » HWALL ♦ 
1HSL»DELPP»HB»DPRAD1 * RE * DPRAD 2 » PEC AV , U I L * U I B > DE LPT >DPFRIC, 
2DPMOM»QSUP*HSUPS*TKSUPO»HLSUP*DELPSU*DPFSUP*DPSMOM 

DATA BITS/0./ 

DATA TEMP/1000. *1100. * 1200* * 1300. *1350# *1400. *1A£>0. *1300. » 

113 30. » 1600. ,1630. *1700. » 175 0. * 18 00. » 1850. » 1900. *1950. » 2000. ,2100. » 

22200. / 

DA TA VISLI/.81».77*.723».69*.67*.65». 64 ».62*.605» 

1.59. . 585*. 56*. 55*. 53*. 51*. 49*. 47*. 45*. 41*. 39/ 

DATA TC LI/27. 55*27. 70 *27. 86*28. 00*28. 10 *28. 17*28. 25* 

128. 32*28. 40 *28. 48, 28. 5 6, 28. 63, 28. 71, 28. 79*28. 87*28. 95 *29. 03* 
229.11*29.27,29.43/ 

DATA RHOL 1/29. 60 *29. 35*29. 10*28. 85*28. 72*28. 60 *28. 48* 

128. 36*28. 25*28. 12*28. 00*27. 88*27. 78*27. 66*27. 54*27. 42* 
227.30*27.18,26.94,26.72/ 

DATA CPLI/1. 0019*1. 000 3* .9 98 7, .9971, .9963*. 9955, 

1.9 947*. 99 39*. 9 930 *.992 3, .99 15 *.9 9 07*. 98 98*. 989 ».9 882* 

2. 98 74. . 9867. .9839. .9843. .9828/ 

DATA VISNA / .556 » .512 ,..476* .444».43Q» .416* 

1.403* .392 ,.382*. 372* .362*. 3 54, .3 46^.338 ,.330, .323* 

2.316* .309 » .296 , .284/ 

DATA TCNA/37. 60*36. 50*35. 50*34. 60*34. 15, 33. 70* 

133. 30 *32. 95, 32. 60, 32. 25, 31. 95, 31. 65 *31. 35 *31. 10* 30. 80* 
230.55,30.35,30.05,29.50*28.90/ 

DATA RHONA/ 50.95*50.1 *49. 3*48. 6, 48. 25*47. 85*47. 50* 

147. 15*46. 80*46. 45, 46. 15*45. 80*45. 50 *45. 15*44.85*44.50* 

244.20*43. 90*43.30*42.70/ 

DATA CPNA/. 3008 * .300 ». 300* .3006, . 3012 ». 3018 » 

1.3028*. 3038*. 3051*. 3064*. 307 9, .3 097, .3116* .3137 * .3160* 

2.3184, . 3210 *.32 37 ,.3300* .3366/ 

DATA VI SKG/. 0412 * .0427 * . 0442 *.04 56*. 0464, .0471* 

1.0 478*. 04 85, .0492*. 049 9, .0505*. 0512*. 0518*. 0524* 

2. 0529*. 0535, .0540 *.054 5, .0555*. 05 62/ 

DATA VISK/. 414*. 379*. 3 54, .336, .327*. 319, .312* 

1.3 05 *.2 98 *.292*. 287*. 281*. 276, .271, .2 66 *.261*. 257*. 253* 

2. 245*. 238/ 

DATA TCK/ 21. 1*20. 1,18. 9, 17. 9, 17. 3, 16. 8, 16. 3, 

115. 8, 15. 2*14. 7, 14. 15*13. 60 *13. 10 *12.50 *12. 05 *11. 50* 
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211, 00 *10. 45, 9. 40 *8.30/ 

DATA REG/3. ,4. ,5. ,6. *7./ 

' DATA EPSMAX/. 477*1. 2455*2. 072. 2. 9775 *3.929/ 

DATA TCKG/. 007 13 *.0074 *.00 76 5* .0 07925 *. 008025* .0081 5 *.00827 5* 
1.0084 *.00852 5 *.00863 *.008 75* .008 81 *.008 97 *.00908 *.009 175 *.009275* 
2. 00935 *.00945 *.00960 *.00971/ 

5000 GO TO 1 
2Q00 FORMAT GENERATOR 
RESTORE 

POWER BOILER TH ERMAL /H YD RAU L I C DESIGN 


SPACE 

TITLE=-X -A 

SPACE 2 

GEOMETRY 


SPACE 



NT = 

-I 


TUBE 

ID = 

-F4 

INCHES 

TUBE 

OD = 

-F4 

INCHES 

WALL 

XW = 

-F4 

INCHES 

TUBE 

LTOT= 

-F4 

FEET 

TUBE 

LB = 

-F4 

FEET 

TUBE 

LSC = 

-F4 

FEET 


SPACE 2 


TUBE P/D » 

-F5 


TUBE FLOW AREA - 

-F6 

FT 2 

INSERT P/D 

-F2 


SHELL ID = 

-F3 

INCHES 

SHELL FLOW AREA= 

-F 4 

FT 2 


THERMAL PARAMETERS 


SPACE 

POWER QTOT 

KW BTU/SEC 

-F0 -F2 

SPACE 2 

PRIMARY 


QB 

BTU/SEC 

-F2 


QSC 

BTU/SEC 

-F2 


Q/A BOILING 
BTU/FT2 HR 
-F2 


SECONDARY 


SPACE 




TPIN 

= 

-F2 

DEG F 

TPOUT 

= 

-F2 

DEG F 

WP 

= 

-F3 

LB M/ SEC 

GP 

= 

-F3 

LBM/FT2SEC 

END OF 

FORMAT 



7000 FORMAT 

GENERATOR 


VP 

= 

-F3 

FT/SEC 

REP 

= 

-FI 


PEP 

= 

-F2 


HP 

= 

-F2 

B/FT2HRDEG 

HWALL 


-F2 

B/FT2HRDEG 

DELTAP 

= 

-F3 

PS IA 

DPRIN 

= 

-F3 

PSIA 

DPROUT 

= 

-F3 

PS IA 


TSI N 

-F2 

DEG F 

TSAT 

-F2 

DEG F 

TSOUT = 

-F2 

DEG F 

WS = 

-F3 

LBM/SEC 


GMAX 

= 

-F 3 

LBM/FT2 

SEC 

VMAX 

= 

-F2 

FT/SEC 


ACC L 

= 

-F2 

FT/SEC2 


ACCL/G 

= 

-F 2 

GEES 


HSL 

= 

-F2 

B/FT2HR 

DEG 

HB 


-F2 

B/FT2HR 

DEG 

REL 

= 

-FI 



PEL 

= 

-F2 



UIL 

= 

-F 2 

B/FT2HR 

DEG 

UIB 

X 

-F2 

B/FT2HR 

DEG 

DELTAP 

= 

-F 3 

PSIA 


DPFRIC 

= 

-F 3 

PSIA 


DPMOM 

= 

-F 3 

PSIA 
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SUPERHEAT REGION 


SPACE 


QSUP 

HSUPS 

TSUPS 

BTU/SEC 

HR UNITS 

DEG F 

-F 2 

-F2 

-F 2 

END 

OF FORMAT 


END 




TUBE LSUP 

DPSUP 

DPSFRIC 

FT 

PS I 

PSI 

-FA 

-F3 

-F 3 


DPSMOM 

PSI 

-FA 
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SIBFTC FRIC FUL 1ST » SDD 

SUBROUTINE FRI C I POD » DCB .DT » G » V I S * FR 1 

1 BETA=SQRT ( 1 «+3 • 14 16/POD** 2 ) 

? GMAX-G*BETA 
8 A=DCB/DT 

4 B=A**2 

5 DENOD=( l.-B ) /( l.+A+( 1 .-A ) /3 . 1416 ) 

6 DEN=DT*DENOD 

7 REQ=DEN*300.*GMAX/VIS 


) 


8 FR1=.316/REQ**.25 
RETURN 
END 
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$ibftc rdp fulist*sdd 

SUBROUTINE RDP ( DT *RDPL ♦ WP »RHO» V IS »PODT »NT »DPR ) 
SUBROUTINE COMPUTES RADIAL PRESSURE DROP. 


501 F=0 • 

XNT = FLOAT (NT ) 

NT1 = 1 
1=1 

4111 NTl=NTl+6*( 1-1 ) 

I F ( NT1-NT 14112*4113 #4114 

4112 1=1 1 

GO TO 4111 
4114 1=1 1 
411? CONTINUE 
C 

DO 502 J=2*I 
NROW=J-l 

F 1 =FLOAT ( NROW) /XNT 
F = F Fl**l .84 

502 CONTINUE 

TERM1=.25 / ( PODT-1 . 1**1. 84 
TERM2=.118/ ( PODT-1 . )**2.92 
GAM=(TERM1+TERM2 )*F 
TERM=(VIS*RDPL/12./WP/3600. ) **. 16 
TERM3=(WP/RDPL*144. /DT)**2 
DPR=TERM/RHO*. 562 /32 . 2 / 1 44. *GAM*T ER M3 
1067 RETURN 
END 



$ I BF-TC DICHOA FUllST.SDD 

SUBROUTINE D I CHOA ( XVAL , UL I M I T * XL I M I T »R A LPT • TOL , SLOPE »K ♦ ANSWER * NI T ) 

0 

c 

C XVAL-VAR I ABLE TO BE ITERATED OR COMPARED 

C ULIMIT=UPPER LIMIT ON ABClSSA 

C XL IM I T -LOWER LIMIT ON ABClSSA 

C BALPT =BALANCE POINT 

C TOL-TOLERANCE 

C SLOPE=SLOPE OF CURVE 

C FOR NEGATIVE SLOPE=l* 

C FOR POSITIVE SLOPE=0 • 

C K=CONVERGENCE INDICATOR=l FOR CONVERGENCE. 

C ANSWER=ANSWER( ABClSSA ) 

C 

I F ( 1-71 ) 5 4 ♦ 2 » 5 4 
2 IF(AFORK) 10,4, 10 

4 AFORK= 1 . 0 

PEMAX=UL I MI T 
PEMIN=XLIMIT 
6 XMAX I =PEMAX 
XM I N I =PEM 1 N 

8 TDIFF=(XMAXI+XMlNl ) /2.0 
ANSWER=TD IFF 
K-0 

RETURN 

C 

C TEST THE DIFFERENCE. 

C 

10 TDIFF=XVAL-BALPT 

I F ( ABS ( TD I FF ) TOL)12,12.20 
C 

C PASS DIFFERENCE EQUAL OR LESS THAN TOL. 

C 

12 AFORK=0. 

I TCTR = 0 
1=0 
K = 1 

RETURN 

C 

20 I F ( SLOPE ) 26 ♦ 22 »26 

22 IF(TDIFF)24,23»23 
C 

C Y GREATER THAN BALPT 

C 

23 XMAX I =ANS WER 
GO TO 32 

C 

C Y LESS THAN BALPT 

C 

24 XM I N I =ANSWER 
GO TO 32 

26 IF(TDIFF) 30,28»28 
28 XMIN I =ANSWER 
GO TO 32 
30 XMAX I =ANSWER 
32 ITCTR=ITCTR+1 

I F ( ITCTR-NIT)8»42 »8 
42 CFORK=0. 

BFORK=0 • 

AFORK=0. 
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itctr=o 

1=0 

CALL DUMP ( TD I FF » 2 ) 
54 1=71 

AFORK=0 • 

BFORK=0 • 

CFORK=0 • 

GO TO 2 
END 
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SIBFTC" LININT FULIST*SDD 

SUBROUTINE L IN I NT ( Z , Y *ARGZ ♦ ARGY , AL ♦ N ) 

' SUBROUTINE FOR LINEAR INTERPOLATION 

DIMENSION Z( 100) *Y(100) 

6 NN=N 

7 M=N 1 
5 DO 1 1 1 = 1 ♦ M 

8 LL=1 +1 

9 L = I 

10 A=(Z(I)-ARGZ)*(Z( I + D-ARGZ) 

12 IF(A) 14 *14*15 
15 AL=0.0 

11 CONTINUE 

20 I F ( ARGZ-Z (1) )25»14*26 

25 IF(Z(1)“Z(2) )27, 14*28 

26 I F ( Z ( 1 )-Z(2 ))28*14,27 

27 ARGY=Y(2)-( ( Z ( 2 ) -ARGZ ) / { Z ( 2 ) -Z ( 1 ) ) ) *( Y ( 2 ) -Y ( 1 ) ) 

GO TO 22 

28 ARGY= ( { ARGZ-Z ( M ) ) / ( Z ( NN ) -Z ( M ) ) ) *( Y ( NN ) -Y ( M ) )+Y ( M ) 
AL=100.0 
GO TO 22 

14 RAT= ( ARGZ-Z ( L ) ) / ( Z ( LL )-Z ( L ) ) 

23 ARGY=Y(L) +RAT* (Y(LL )-Y(L) ) 

24 AL=L 
22 M=N 1 

RETURN 
END 
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$IBFfC DUB IN FULIST»SDD 

SUBROUTINE DUBIN(P*T 


♦V ♦ PLUSH *TU5H *VUME *MM*NN»JA*JB) 


SUBROUTINE FOR DOUBLE LINEAR INTERPOLATION 


DIMENSION PLUSH (40) * TUSH (40) *VUME(1600) 
100 PP=P 
200 TT=T 
600 M=MM 
700 N=NN 

800 DO 1100 I =1 ♦ M 
850 L=I 

900 C= ( PLUSH ( I ) -P ) * ( PLUSH ( I 1 ) -P ) 

1000 IF(C)2700 *2700*1200 
1200 JAM=3 
1100 CONTINUE 

1400 TEST = (PLUSH (1)-P)*(PLUSH(2)-PLUSH(1 )) 
1500 IF(TEST) 1700*2100*2500 
1700 JAB*100 
1750 MN-M-1 

1800 RAP- (P— PLUSH (MN ) )/( PLUSH ( M ) -PLUSH (MN ) ) 
LL-MN 

1950 GO TO 2800 
2100 RAP*0 *0 
2200 LL=1 

JAB=200 
GO TO 2800 

2500 RAP= ( P-PLUSH ( 1 ) )/( PLUSH ( 2 ) -PLUSH ( 1) ) 

2550 JAB=0 
2600 LL=1 
2250 GO TO 2800 
2700 LL-L 

LLL=LL+1 
2950 JAB*LL 

2900 RAP* ( P-PLUSH ( L L ) )/( PLUSH ( LLL ) -PLUSH ( LL ) ) 
2800 LLL=LL+1 
3000 DO 3500 I=1*N 
3100 K = t 

3200 C*tTUSH( I )-T)*(TUSH( I+1J-T) 

3300 IF (Cl 4500*45 00 *3400 
3400 JAM*10 
3500 CONTINUE 

3600 TEST * ( TUSH ( 1 ) T ) * ( TUSH ( 2 ) -TUSH ( 1 ) ) 

3700 IF (TEST ) 3800 *4100 *4300 
3800 N J=N-1 

3900 RAT = ( T-TUSH ( NJ ) ) /( TUSH ( N ) -TUSH ( NJ ) ) 

KK=N J 

3850 JBB=100 
4050 GO TO 4600 
4100 RAT*0.0 
4200 KK=1 

JBB*200 
GO TO 4600 

4300 RAT = (T-TUSH ( 1 ) ) / ( TUSH ( 2 ) -TUSH ( 1 ) ) 

4350 JBB=0 
4400 KK=1 
4250 GO TO 4600 
4500 KK=K 

KKK=KK+1 
JBB=KK 

RAT= ( T-TUSH { KK ) )/( TUSH ( KKK ) -TUSH ( KK ) ) 


4675 

4650 
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4600-KKK=K+l 
4700 KV1=M*KK+LL-M 
4*800 KV2=KV1+1 

4900 V1=VUME(KV1)+RAP*(VUME(KV2)-VUME(KV1) ) 
KV1=M*KKK+LL-M 
5100 KV2=KV1+1 

5200 V2=VUME(KV1)+RAP*(VUME(KV2)-VUME(KV1) ) 
5300 VV=V1+RAT*(V2 VI) 

5400 V=VV 
5500 JA=JAB 
5600 JB=JBB 
RETURN 
END 
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SIBFTC CURV FuLlST,£DD 
FUNCTION CURV(P) 

17 PYP=P/14. 69594 

18 PUN=(ALOG(PYP))/2. 3025851 

19 PEEP= (PUN+0. 10705539 ) /2. 1805446 

20 PUP =PEEP* PEEP 

21 PUPPY =PUP*PEEP 

22 PURP=PUP*PUP 

23 PURPLE =PUP*PUPPY 

24 PURSE*PUPPY*PUPPY 

25 PERK=PURP*PUPPY 

26 PERKY=PURP*PURP 

27 PLUP=1 *3 1902 32*PEEP-» 13299096+ .6923 7968*PUP+. 37008034 
1*PUPPY+. 1 9876 127*PURP+.10754778*PURPLE+. 066 299672 
2*PURSE+*041265195*PERK+.013472973*PERKY 

28 CURV=1900. 0+700. 0*PLUP 
RETURN 

END 



sibftc thfrma fulist*sdd 

SUBROUT INF THERMA f PR » T * VPS I A *HL L * HFN * SL * SEN * VOLLL * VENN * CPV ♦ CPL ♦ 
IPOLLW 9 AV *CHOICE*SUPER* HGG *SG*VV*CPP*POLLY , A) 

C SUBROUTINE RETURNS THERMODYNAMIC PROPERTIES 

C OF SATURATED AND/OR SUPERHEATED POTASSIUM. 

C (NRL PROPERTIES). 

C SUPER=-1. # sp. VOL»ENTHALPY»AND ENTROPY ONLY. 

C SUPER=0.9 SATURATED PROPERTIES ONLY. 

C SUPER=1.» SUPERHEATED PROPERTIES. 

C 

C CHO I CF = 1 * » PRESSURE INPUT ( SATURATED ) 

C C HO I C F = - 1 • * TEMPERATURE INPUT < SATURATED ) 

C C HO I C F = A B S GREATER THAN 6736.0 9 S p VOL 

C ITERATIONS PRINTED. 

C C HO I C F = A R S GREATER THAN 1362.0 9 VIRIAL COEF . 

C AND THEIR DERIVATIVES PRINTED. 

C 

DIMENSION PLUSH*; 1 8 ) » TUSH ( 12 ) » VUME { 3 00 ) »FF< 1 00 ) 

' DATA PuUSH/10. »2 . 9 1 . 9 .3333 9 . 2 9 . i 66667 9 . 125 » .1 » 

10. 083333333*0. 07 142 857*0. 0625*0. 055555556*0. 05*0. 03846 1538* 

20. 33333*0. 0?5»0. 016667*0.0125/ 

DATA TUSH /8 . 33 33 3 E -4 9 7 . 66923077 F 4*7.142857E~4*6.666667E-4» 

1 6. 2 5 E- 4,6 0 8823 5? 9 E 4 9 5 . 5 5 555 55 E - 4 * 5 . 26 3 1 57 F - 4 9 5 . OF-4 9 
14. 76 1 90-+E - -4 9 <+» 54 5 4 F - 4 9 4 « 3478 E- 4 / 

DAT A f VUME t I 1 9 1 -■ L 9 2 08 / 3 ? 54 . 8 > 6 1 8 . 1 9 2 85 . 2 9 26 6 • 9 2 4 1 • 9 22 8 . 9 2 09 . * 90 . 9 
1 1 7 1 . 0 9 1 5 2 . 0 ■> 1 3 3 . 0 9 1 1 4 . 0 *95 . 0 9 76 .0 9 5 7 . 0 9 38 . 0 9 1 9 . 0 9 1 . 0 9 3 545 . 5 9 
2691. 49333.69 9 I.I 6938 . 94*24. 71 » 5. 18» 20. 0.18. O 9 I 6 .O 9 I 4 .O 9 I 2 .O 9 IO.O* 
48. 096.0*4.0*2. 0*10. 0*3828. 1*75 5. 13*370 .62 *112. 84 *59. 89 *46. 22* 

428. 48*17. 142»9. 005*2.696*2. *1.8*1. 6*1.491. 2 *1»*»8*. 6* 

54106. 9*814. 8*403. 1*128. 0*72. 37*58. 28*40. 40 *29. 37*21. 77*16. 13*12** 
68. 100,5. 10 s 3. 10*1. 60* 1.30*1.00* .7*43 83. 9*872. 4*433. 3*1 40. 4* 

781 .4 7 966®64*48o02 *36.7 929.03 *23.4 59 19.18915 .78912. 99*6.86* 

83. 90 8*2. 0*1. 0*0. 5*4659. 9*928. 9* 46 2. 5*151. 4*89. 02*73. 38,53. 76* 

94 1.9 2 *33. 96 9 2 8*23*23. 88*20. 46*17. 68*1 1.75*8. 99*4. 18*2. 0*1.0* 

9 4 9 3 5 . 2 * 9 8 4 . 8 9 4 9 1 . 0 * 1 6 1 . 7 » 9 7 • 7 6 9 7 9 • 2 5 * 5 8 . 5 7 * 46 . 1 2 * 3 7 . 7 9 * 

131. 8*27. 29*23. 7 6*20. 91*14. 91*12. 17*7. 53*2. 39*1. 0*5210. 3*1040. 4* 

25 19. 16*171.6 I * 102.05* 84.65*62. 87* 49. 78 *41. 03* 34. 7 6*30.04* 26. 36* 
323.4] ,17.21 » 14.41 *9.75 *4.49* 19.99 *5 485.1 * 1 0 95 . 8 * 547 . 09 * 1 8 1 . 2 8 * 

41 08. 09* 89. 78 * 66. 88 * 53. 22 * 4 3. 94 * 37. 3 7*32. 43 * 28. 58 *25. 49 * 19. 02 6 * 
516. i 52 *11. 383*6. 419*1. 977* 5 759. 7 *1151. 0*574. 9*1 90. 8* 114.0* 

694. 74 * 70. 71 * 56 . 28 * 46. 65 *39. 77*34. 6 * 30. 57 *27. 34*20. 605* 17. 611 * 
712.68*7. 62*4. 96* 6034. 3 * 1206. 1,6 02.6 *200.2 * 119. 7 *99. 58*74. 42 ♦ 

859. 31* 49. 24* 42. 03* 3 6 . 63*32. 42*29. 05*22. 02*18. 9* 13. 78*8. 54*5. 88 * 
96308.8*12 61.1*630.2*209.5*125.4*104.4*78.05 *62.27*51.74*44.22/ 
DATA (VUME (I 5* I =209, 216)/ 38. 57*34. 17*30. 66 *23. 33*20. 08* 14. 76 *9. 35* 
16.61/ 

200 PNRL=PR/14. 69594 

230 IF (PNRL-0. 007 )231 *232*232 

231 PNRL=0.007 

232 IF (PNRL-21. 0)233*233*234 

234 PNRL=2 1 . 0 

233 AJ=778.26 
E =67 • 0 

F = 68 . 0 

G=69.0 

DECKA=0 

DECKB=0 

DECKC=0 

VALA=0 

VALB=0 

VALC=0 
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CHECK-CHOICE 

J=3 

JJ=10 

PRENT=ABS (CHECK ) 

300 TNRL=T 

330 IF ( TNRL-1250. 0 ) 331 ♦ 332 »332 

331 TNRL=1250 *0 

332 IF (TNRL-2700. 0)350,350,333 

333 TNRL=2700.0 
350 KAM=0 

500 TLAG=(AL0G(TNRL) )*. 4342945 

600 I F ( CHOICE )800»800»820 
820 PSAT AT =PNRL 

GO TO 1000 

800 PVLAG= (6.12758 -8 128.77 /TNRL-0.5 32 99 *TL AG )*2.302585 
900 PSATAT=EXP ( PVLAG ) 

1000 PSAT=PSATAT*14. 69594 
1100 PSATSF=PSATAT*2116.2 

1200 BANLUG* (4544. 3/TNRL-3. 7641+TLAG )*2. 302585 
1300 BAN=EXP ( BANLUG ) 

1400 B=BAN*(-1 .0 ) 

1500 CLUG=(0. 6988+5725. 0/TNRL)*2. 362585 
1600 C=EXP(CLUG) 

1700 DANLUG=(1. 9274+5671. 4/TNRL)*2. 302585 
1800 DAN=EXP ( DANLUG ) 

1900 D=DAN*(-1 .0 ) 

2000 DBDT=(B/TNRL)*(1. 0-10463. 637/TNRL) 

2100 TSQ=TNRL*TNRL 

2200 DCDT=13182.2*C*(-1.0 )/TSQ 

2300 DDDT =13058.8 *D*(-1«0) /TSQ 

2400 TCUBE=TSQ*TNRL 

2500 ELLA=39375.0/TNRL 

2600 ELLEN=EXP (ELLA) 

2700 HZERA=998.95 + . 127*TNRL + 24836. 0/ELLEN 
2800 SARAH=31126.0/TNRL 
2900 SALLY=EXP (SARAH) 

3000 SZERA=0. 18075 0 . 1 27*TILAG*2 . 302 585+0 .761 7/SALLY 
3100 RR=1545.43 
3200 WAYT = 39.10 
3300 RRM=39. 52506 

3400 PARLUG* ( 408 0 .0 /TNRL-1 . 803 ) *2 . 3025 85 
3500 PAR* EXP ( PAR LUG ) ' 

3600 VAL1* ( RRM*TNRL /PSATSF-PAR ) *39. 1 
MARK=1 

8200 TFAR=TNRL-459. 69 
8300 TFAR2-TFAR*TFAR 

8400 TFAR3=TFAR*TFAR2 * 

8500 DUNSL=52 .7274 7. 3539E-3*TFAR-0 . 5601 3E-6*TFAR2+0. 03 158 
1 E-9*TFAR3 
8600 VLAQ=1 .O/DUNSL 
8700 VOLLL=VLAQ 
9500 T4=TSQ*TSQ 
9700 D2B=1.094877E8*B/T4 

9800 D2C=26364.4*C/TCUBE-13182.2*DCDT/TSQ 
9900 D2D=26117.6*D/TCUBE-13058.8*DDDT/TSQ 

601 EMMA=28070. O/TNRL 
701 EMMY=EXP ( EMMA ) 

801 CZERA=0. 127+2. 888/EMMY 

CPLAQ=0.22712 0*64848E-4*TNRL+0.231 78E-7*TSQ 
CPL=CPLAO 

IF (PRENT-1362. 0)3700, 102 ,102 


- 113 - 



102 WRITE (6*101) B,C*D*DBDT»DCDT»DDDT,D2B, 
1D2C*D2D*PSAT »TNRL 

3700 VAL1 2=VAL 1*VAL 1 

3800 VAL1 3=VAL1*VAL 12 

3900 SLEFT=PSATSF*VA|_1/ ( 1545.43*TNRL ) 

4000 R IGHT=1 *0+B/VALl +C/VAL12+D/VAL 1 3 
4100 TALLY=0.001*SLEFT 

IF (PRENT-6736. 0)4200 *103 *103 

103 WRITE (6*104) VAL1 ♦ SLEFT *R IGHT ♦ DECK A*DECKB ♦ 
1DECKC*VALA,VALB,VALC»CZERA*HZERA*SZERA*PSAT*TNRL 

4200 DECK=SLEFT-R IGHT 
4250 ADECK=ABS(DECK ) 

4300 IF (TALLY-ADECK >4400*4500 *4500 
4400 I F ( MARK-1 0 ) 4600*4700 ,4700 
4700 VAL1=(RRM*TNRL/PSATSF~PAR)*39.1 
GO TO 4500 

4600 I F ( MARK-2 )4800*5200,5500 
4800 VALA=VAL1 
4900 DECKA=DECK 
5000 VAL1 =VALA*1 *03 
5100 MARK=MARK+1 
GO TO 3700 
5200 VALB=VAL1 
5300 DECKB=DECK 
5400 VAL1 =VALA*0* 97 
MARK=MARK+1 
5450 GO TO 3700 
5500 VALC=VAL 1 
5600 DECKODECK 
5700 DECKA2=DECKA*DECKA 
5800 DECKB2=DECKB*DECKB 
5900 DECKC2*DECKC*DECKC 
6000 Z 1Z3=DECKA-DECKC 
6100 Z1Z2-0ECKA-DECKB 
6200 W1W3=DECKA2-DECKC2 
6300 W1W2=DECKA2-DECKB2 
6400 Y1Y3=VALA-VALC 
6500 Y1Y2=VALA-VALB 
6600 ANUM=Z1Z3*Y1Y2-Z1Z2*Y1Y3 
6700 ADEN=W1W2*Z 1Z3-W1 W3*Z 1Z2 
6800 AAA=ANUM/ADEN 
6900 BBB=(Y1Y2-AAA*W1W2 ) /Z1Z2 
7000 VAL1=VALA-AAA*DECKA2-BBB*DECKA 
7100 MARK =MARK+1 
7400 DECKA*DECKB 
7500 DECKB=DECKC 
7600 VALA=VALB 
7700 VALB=VALC 
GO TO 3700 

4500 VENNA=VALl/39. 1 
7800 VAL2=VAL1*VAL1 
7900 VAL3 =VAL 1 *VAL2 

7950 VAL4=VAL2*VAL2 

7951 VAL5=VAL?*VAL3 

8000 HG1=HZERA+(RRM*TNRL/AJ)*( ( B-TNRL*DBDT ) / VAL1 + ( C-TNRL 
l*DCDT/2.0 )/VAL2+(D-TNRlL*DDDT/3.0) /VAL3 ) 

8100 HENNA=HG1 

8150 HGSAT =HENNA 

9200 VRT=VAL1/(RR*TNRL) 

VRTAB=ABS ( VRT ) 

VRT =VRTAB 
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9300 VRTLUG=ALOG( VRT) 

9400 SGSAT*SZERA- ( RRM/A J ) * ( ( B+TNRL*DBDT ) /VAL1-VRTLUG 

1+(C TNRL*DCDT ) / ( 2 • 0*VAL2 ) + ( D+TNRL*DDDT )/( 3. 0*VAL3 ) -7.65738 ) 
9600 SENNA=$GSAT 
VALSAT=VENNA 

8800 HATTY=(PSATSF/AJ )*( 1 8717 • 2/TNRL-0 • 5 3299 ) *( VALSAT-VLAQ ) 

8900 HLNRL=HGSAT-HATTY 

satty=hatty/tnrl 

SLNRL=SENNA-SATTY 
IF(SUPER) 2001.201 .201 

201 BRAKET=1 . 0+ ( B TNRL*DBDT ) /VAL 1+ ( C+TNRL*DCDT ) /VAL2 
1+ ( D TNRL*DDDT ) /VAL3 
301 CGNUM=BRAKET*BRAKET 

401 CGDEN=1.0+2 .0*B/VAL1+3.0*C/VAL2+4.0*D/VAL3 
501 PAR2=TNRL*D2B 2 . 0*DBDT+0 . 5* ( TNRL*D2C+2 . 0*DCDT ) /VAL 1 
1+(TNRL*D2D+2.0*DDDT )/(3.0*VAL2 ) 

901 CPG=CZERA+(RRM/AJ ) * ( CGNUM/CGDEN-1 .0-( TNRL/VAL 1 ) 

1#PAR2 ) 

1050 BRUCK=1.0/VAL2+2 . 0*B /VAL3+3 . 0*C/VAL4+4. 0*D/ VAL5 
HOI DRHDT =RR*TNRL* ( -1 • 0 ) *BRUCK 
1201 DVDP = 1.0/ (DRHDT*WAYT ) 

1301 DVDT* (PSATSF/ ( RR*TSQ ) +DBDT/VAL2+DCDT/VAL3+DDDT/VAL4 ) 

1 / ( BRUCK*WAYT ) 

1401 PALE=1.0/{ (PSATSF/VENNA)*(DVDP+(TNRL/(AJ*CPG) ) *DVDT*DVDT ) ) 
1501 PALM=PALE*(-1.0) 

1601 ALF=PALM*32.174*RRM*TNRL 
1701 ABALF=ABS(ALF) 

1801 ALE-SORT ( ABALF ) 

1901 I F ( KAM ) 21 01 .2001 *2101 
2001 VENN=VENNA 
2201 HEN=HFNNA 
2301 SEN^SENNA 
- SL-SLNRL 
9000 HLL=HLNRL 
9100 VPSl A=PSAT 

2351 IF (SUPER) 2352.2401,2401 

2352 CPV=0 .0 

2353 POLLYV=0.0 

2354 AV=0.0 

GO TO 2801 
2401 CPV=CPG 
2501 POLLYV=PALM 
2601 AV=ALE 

2701 IF(SUPER) 2801, 2801, 3401 
2801 HGG=0 
2901 SG=0 
3001 VV=0 
3101 CPP=0 
3201 POLLY=0 
3301 A=0 
3351 GO TO 4601 
3401 KAM=1 
3501 MARK=1 
3601 PSAT AT =PNRL 
3701 PSAT =PSAT AT *14.69594 
3801 PSATSF=PSATAT*2116.2 
MM=1 8 
NN = 1 2 

RPSAT=1 .0 /PSAT 
RTNRL=1 .O/TNRL 

3903 CALL DUBIN(RPSAT ♦ RT NR L > VELD » PLUSH » TUSH » VUME ,MM,NN* JAN, JAQ) 
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MARK=1 

4001 VAL1=VELD*39.1 
GO TO 3700 
2101 HGG=HENNA 
4101 SG=SENNA 
CPP=CPG 

4201 VV=VAL1/39.1 
4301 CPP=CPG 
4401 POLLY=PALM 
4501 A-ALE 
4601 HARRY=6.0 
C 


RETURN 

101 FORMAT GENERATOR 


J- .1. • ' ' ' ' 

B 

-1PE3 

SPACE 

C 

-E5 

D 

-E5 

DBDT 

-E5 

DCDT 

-E5 

DDDT 

-E5 

TEMP 


D2B 


D2C 


D2D 

PRESS 



-E5 


-E5 


-E 5 


-E5 


-E5 


SPACE 
END OF 

FORMAT 










104 FORMAT 

GENERATOR 









SPACE 

VALl 

E 5 

LEFT 

-E5 

RIGHT 

-E5 

DEC1 

-E5 

DEC2 

-E5 

SPACE 

DEC3 


V0L1 


V0L2 


V0L3 


CO 

-E5 


— 

E 5 


-E5 


-E5 


-E5 


SPACE 

HO 


SO 


PRESS 


TEMP 





— 

■E5 


-E5 


-E5 


-E5 



SPACE 
END OF 
END 

FORMAT 
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ABSTRACT 


Calculation procedures and digital computer codes for thermal design 
of multiple-tube potassium condensers and boilers are presented. The 
computer programs are given, together with typical input and output formats 
as illustrated by representative test cases. The computer codes are useful 
for preliminary thermal design studies of condensers and boilers applicable 
to Rankine cycle space power systems using potassium as the working fluid. 



